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PREFACE 


Work  in  this  paper  is  primarily  aimed  at  the  elementary  development  of  the  theory 
of  time-harmonic  electromagnetic  fields  that  have  the  force-free  form  for  the  electric  and 
magnetic  fields,  as  well  as  the  current.  It  is  hoped  that  enough  detail  is  included  for  the 
paper  to  be  used  as  a  tutorial  for  the  majority  of  people  who  have  never  encountered  force- 
free  fields,  and  wish  to  investigate  what,  at  first  blush,  seems  like  an  impossible  solution 
of  Maxwell’s  equations. 

The  research  was  spurred  by  attempts  to  understand  experimental  results  obtained  in 
the  late  1960’s  and  early  1970’s  in  areas  of  fusion  weapons  and  laser  fusion  in  Los  Alamos. 
At  the  time  there  were  anomalies  such  as  unexplained  symmetries  in  the  fragments  of 
disintegrating  laser-irradiated  pellets,  and  puzzles  about  compression  and  stability  of  small 
spherical  volumes  of  plasma  in  the  presence  of  high-intensity  electromagnetic  fields.  During 
that  time,  and  since,  a  large  number  of  people  investigated  many  of  the  phenomena  of 
interest,  without  fully  bong  able  to  resolve  all  the  challenges  of  plasma  inertial  confinement. 
No  claims  will  be  nude  that  this  work  will  make  any  of  the  inertial  confinement  problems 
easier. 

In  the  period  of  close  to  two  decades  after  performing  the  original  work,  some  of 
the  insights  gained  into  solutions  of  Maxwell’s  equations  in  a  conducting  fluid  were  seen 
to  have  relevance  beyond  the  controlled  fusion  community.  In  particular,  finite-sized, 
non-plane-wave,  force-free  time-harmonic  electromagnetic  excitations,  with  electric  and 
magnetic  field  vectors  parallel  to  each  other  everywhere,  in  the  presence  of  plasma  current, 
seemed  to  be  undiscovered  solutions  that  deserved  exploration.  What  makes  such  solutions 
interesting  is  that  the  Poynting  vector  is  identically  zero  everywhere  inside  a  finite-sized 
force-free  field,  and  it  cannot  dissipate  unless  either  its  boundary  is  radiating  or  internal 
stresses  make  it  break  up. 

Classically,  these  fields  cannot  exist  in  a  vacuum.  In  a  quiesdent  plasma,  internal 
stresses  on  their  accompanying  currents  can  disrupt  them  as  soon  as  they  are  created. 
By  a  remarkable  coinddence,  the  electromechanical  stresses  manifest  themselves  mainly 
as  a  pressure  field,  which  can  be  nearly  exactly  balanced  by  vortical  fluid  motion.  This 
vortical  fluid  field,  otherwise  called  Beltrami  motion,  is  also  self-disruptive,  by  itself,  in  a 
compressible  fluid. 

The  alleged  simultaneous  occurance  of  the  force-free  time-harmonic  electromagnetic 
field  and  the  matching  vortical  fluid  motion,  with  a  well-defined  boundary,  will  be  called 
a  “plasmoid".  Although  the  solutions  to  Maxwell’s  equations,  coupled  to  fluid-dynamic 
modes,  do  not  seem  to  violate  any  physical  prindples,  it  is  not  immediately  obvious  how 
one  would  go  about  producing  them.  Mathematically,  it,  at  first,  seemed  possible  that 
such  excitations  could  exist  at  all  frequencies.  A  more  physical  approach  showed  that 
below  a  critical  frequency,  determined  by  its  size,  a  plasmoid  represents  a  plasma  state  at 
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a  local  energy  ir|WY^TT1nm:  i.e.  it  is  unstable.  In  addition,  electrons  are  unable  to  carry  the 
current  below  this  critical  frequency.  Above  the  critical  frequency,  the  physics  becomes 
more  interesting,  as  the  current  is  carried  by  the  electrons  and  the  plasma  is  at  an  energy 
minimum,  even  though  the  total  energy  is  positive. 

The  author  believes  that  a  number  of  observed  phenomena,  including  ball  lightning, 
“fireballs”  and  “charge  clusters”  that  occur  in  strong  plasma  discharges,  may  be  described, 
at  least  in  part,  by  the  plasmoids  in  tins  theory.  Some  of  the  non-linear  effects  of  powerful 
lasers  impinging  on  small  spheres  of  material  might  also  produce  plasmoids  of  the  type  to 
be  described. 

The  work  in  this  paper  being  only  a  start,  the  theory  is  incomplete.  A  great  deal  of 
work  still  remains  to  be  done  to  determine  whether  real  plasmas  at  elevated  temperatures 
can  sustain  force-free  plasmoids.  Open  questions  of  quantum  and  quantum  electrodynamic 
effects  and  interactions  with  the  electromagnetic  zero-point  background  remain  to  be  in¬ 
vestigated.  Puzzles  concerning  external  fields,  surface  currents  and  plasmoid  interactions 
«lan  need  solving,  and  may  involve  quantum  effects.  The  possibility  that  unexplored  and, 
as  yet,  unknown  principles  are  involved  in  explaining  the  existence  of  compact,  stable,  au¬ 
tonomous,  highly  energetic  concentrations  of  plasma  and/or  electromagnetic  fields  cannot 
be  overlooked. 

Much  of  this  work  can  be  conducted  with  quiet  thought  in  libraries  and  in  front  of 
a  computer.  In  general,  however,  a  complementary  experimental  program  is  absolutely 
necessary  at  this  point  to  focus  the  theoretical  effort. 

The  author  would  like  to  thunk  Dr.  franklin  B.  Mead,  Chief  of  the  Future  Technologies 
Section  in  the  Astronautical  Sciences  Division  at  the  Air  Fbrce’s  Phillips  Laboratory  for 
the  the  interest  he  has  shown  in  this  work.  Without  Dr.  Mead’s  encouragement  and  kind 
hospitality  the  ideas  in  this  work  may  never  have  been  developed  even  this  far.  Research 
in  such  areas  is  very  risky  and  I  acknowledge  Dr.  Mead’s  bravery,  in  an  era  of  rapidly 
diminishing  funds,  for  the  support  of  the  work.  Thanks  are  also  due  to  the  University  of 
Dayton  Research  Institute  and  Dr.  Eugene  Gerber  for  providing  a  pleasant  atmosphere 
for  conducting  research. 
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FORCE-FREE  TIME-HARMONIC  PLASMOIDS 


SUMMARY 

Maxwell’s  equations  allow  time-harmonic  solutions  wherein  the  electric  and  magnetic  fields 
are  parallel  to  each  other  and  both  have  the  force-free  form.  These  fields  are  accompanied 
by  a  time-harmonic  current  field  that  also  has  the  force-free  form.  Such  solutions  are  pos¬ 
tulated  as  a  model  for  long-lived  spherical  plasmoids.  It  is  shown  that  a  critical  frequency 
exits,  which  depends  on  the  size  of  the  plasmoid.  Below  this  frequency,  the  time-harmonic 
current  cannot  be  carried  by  electrons  in.  a  plasma.  Above  the  critical  frequency  local 
perturbations  of  the  current  tend  to  increase  the  total  energy,  lending  stabilization.  The 
electromechanical  stresses,  induced  by  the  fields  on  the  current-carrying  electrons,  can  be 
balanced,  to  a  high  degree  of  accuracy  by  vortical  motion  of  the  supporting  plasma.  The 
velocity  field  of  the  vortex  also  has  the  force-free  (Beltrami)  form.  This  vortical  fluid 
motion  is  stable  and  exactly  compatible  with  the  stresses  in  the  plasmoid.  Although  it  is 
a  theoretical- possibility  for  currents  to  flow  on  the  outside  of  a  plasmoid  to  isolate  it  from 
external  fields,  free-space  time-harmonic  fields  can  be  found  that  have  both  their  electric 
and  magnetic  fields  parallel  to  each  other  over  the  surface  of  a  spherical  plasmoid  having  a 
resonant  size,  making  the  fields  continuous  at  the  boundary.  Surface  currents  can  thereby 
be  eliminated.  Only  at  a  resonant  size  can  the  outward-going  radiation  due  to  the  fields 
at  the  surface  of  the  plasmoid  be  nulled  and  the  external  fields  be  joined  continuously  to 
the  internal  fields.  An  example  is  given  showing  that  a  one-megajoule  force-free  plasmoid 
in  an  STP  atmosphere  could  be  described  with  physically  reasonable  parameters. 


I.  INTRODUCTION 

Bostick  described  a  plasmoid  as  a  structure  made  of  plasma,  “. . .  whose  form  is  de¬ 
termined  by  the  magnetic  field  it  carries  along  with  itself**)  .*  His  definition  attempted 
to  describe  the  coherent  discharges  he  encountered  emanating  from  a  plasma  gun.  These 
plasma  structures  are  evidently  very  stable,  have  long  lifetimes,  and  can  interact  with  each 
other*2*.  Since  the  time  of  Bostick's  original  work,  plasmoid-like  structures  have  been 
reported  by  other  investigators***.  Under  different  conditions  they  may  have  names  like 
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“anomalous  discharges”,  “charge  clusters”,  and  “ball  lightning.” 

The  phenomenona  that  seem  to  be  shared  by  all  the  occurences  of  plasmoids  are  the 
presence  of  strong  electromagnetic  fields  and  a  current  Sow  near  the  time  and  position  of 
plasmoid  creation. 

There  may  be  a  number  of  distinct  species  of  plasmoids,  each  deriving  its  stability  from 
different  physical  principles.  These  principles  may  involve  quantum  as  well  as  quantum- 
elect  rodynamic  effects.  It  is  possible  that  some  undiscovered  physical  laws  govern  plasmoid 
generation.  What  .will  be  done  in  this  work,  however,  is  to  assume  that  at  least  one  family 
of  plasmoids  can  be  described  by  a  classical  model  that  invokes  only  “standard”  physical 
theory.  The  physical  model  to  be  presented  for  plasmoids  may  not  be  complete,  t.e.,  it 
may  have  deficiencies.  It  may  not  describe,  even  approximately,  all  species  of  stable  plasma 
excitations  lumped  under  the  term  “plasmoids.”  There  are,  however,  compelling  reasons 
to  explore  a  simple  theory  based  on  an  elementary  model,  which  appears  to  be  relevant 
to  stable  plasmoid  creation  and  subsequent  evolution.  Pushing  the  theory  to  its  limits,  its 
deficiencies  may  become  evident.  At  that  point,  suggestions  can  be  made  for  improving 
the  model,  and  where  new  physics  has  to  be  introduced. 

A  force-free  time-harmonic  electromagnetic  wave  trapped  in  a  stabilizing  plasma  vor¬ 
tical  mode  is  chosen  as  a  model.  The  functional  forms  of  the  tiine-varying  electric  and 
magnetic  vector  fields  are  identical,  within  a  constant  factor  and  phase,  with  this  form 
bring  identical  to  that  of  the  fluid  velocity  field  of  the  plasma.  The  major  electromagnetic 
stresses,  manifested  by  pressures  proportional  to  the  local  electric  field  energy,  are  bal¬ 
anced  by  reduced  pressures  due  to  fluid  motion.  The  pinch  effect  due  to  the  time-varying 
current  adds  to  the  stability. 

Using  a  standing-wave  model  is  not  new.  Kapitza  conjectured  that  ball  lightning  was 
explained  by  such  a  model*4).  A  model  having  a  vortical  field  identical  to  a  component  of 
the  electromagnetic  field  is  not  new  either.  Wells  explored  such  a  model  a  number  of  years 
ago*5).  What  seems  to  be  new  is  the  combination  of  the  ideas  to  explain  some  aspects  of 
plasmoid  phenomena. 

After  some  mathematical  preliminaries,  it  will  be  shown  that  time-harmonic  solutions 
to  Maxwell’s  equations  can  have  the  electric  and  magnetic  fields,  as  well  as  the  current  in 
the  force-free  form.  It  will  be  shown  that  the  phenomena  of  interest  do  not  occur  below  a 
critical  frequency,  which  is  determined  by  the  size  of  the  plasmoid. 

II.  MATHEMATICAL  PRELIMINARIES 

A  vector  field,  F,  will  be  called  “force-free”  if  it  obeys  the  equation 

VxF  =  01F,  (1) 
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where  a  can  be  a  complex-valued  function  of  position.  Unless  otherwise  noted,  however, 
a  will  be  assumed  to  be  constant  in  time,  spatially  uniform,  and  real,  inside  regions  that 
will  be  called  plasmoids.  Equation  (1)  will  not  necessarily  be  assumed  to  hold  outside  the 
boundary  of  a  plasmoid. 


At  first  sight,  it  may  seem  that  the  only 
solution  to  eq.(l)  is  F  =  0.  The  analogy  usu¬ 
ally  brought  to  mind  is  that  of  the  magnetic 
field  due  to  a  steady  and  uniform  direct  cur¬ 
rent  flowing  in  a  wire.  In  fig.(l)  a  steady  cur¬ 
rent,  J,  with  units  of  amperes  per  square  me¬ 
ter,  is  assumed  flowing  in  an  infinite  cylindri¬ 
cal  wire,  producing  a  magnetic  field  H  inside 
and  outside  the  wire.  From  Maxwell’s  equa¬ 
tions,  VxH  =  J  inside  the  wire.  It  is  quite 
obvious  that  H  and  J  are  not  parallel  to  each 
other  anywhere. 

In  general,  however,  it  is  possible  to  cre¬ 
ate  vector  functions  for  which  eq.(l)  is  true. 
Taking  the  curl  of  eq.(l)  inside  a  plasmoid 
results  in 


Figure  1.  Uniform  Current  in  a  Cylindrical  Wire. 
Inside  the  wire  the  magnetic  field  is  orthogonal  to 
the  current. 


Vx  VxF  =  V(V-F)  —  VJF  =  a2F.  (2) 

From  eq.(l),  F  arises  from  the  curl  of  itself  so  that  V-F  =  0.  Therefore  it  follows  that  F 
satisfies  the  vector  Helmholtz  equation 


(VJ+aJ)F  =  0.  (3) 

Suitable  solutions  to  eq.(3)  in  sphurical  coordinates  can  be  obtained  from  solutions  of 
the  scalar  Helmholtz  equation 

(VJ  +  as)^(r)  =  0,  (4) 

using  the  formulas*8* 


and 


L(r)  =  VV>(r), 

(5) 

M(r)  =  L(r)  x  r, 

(6) 

N(r)  =  ±VxM(r), 

(7) 
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where  r  is  the  radius  vector  from  the  origin  to  the  field  point.  The  solenoidal  vector 
P  =  M  +  N  satisfies  eq.(l). 

Although  this  prescription  will  produce  a  force-free  vector  field  for  each  </>  satisfying 
eq.(4),  it  is  beyond  the  present  scope  of  the  present  work  to  investigate  the  mathematical 
completeness  of  this  procedure. 

A'  concrete  example  of  a  family  of  force-free  vector  fields  is 

Mj(r,  9, 4)  =  (1  + 1)  coe  9  sin1-1 9^^a—r  —  [ii_i(ar)  -  2——]  sin*  99 

otT  or  [oj 

±  ji(ar)  sin1 9$, 

where  2  is  any  integer,  2  >  0,  ji  is  the  familiar  spherical  Bessel  function,  and  the  vectors  r, 
9,  and  ^  are  the  familiar  right-handed  3-D  polar  coordinate  unit  vectors,  with  <j>  the  polar 
angle.  The  vector  function  M|(r,  9,  </>)  is  axially  symmetric,  t.  e.,  there  is  no  dependence  on 
4>  in  any  of  the  vector  components. 

Figures -(2-4)  show  families  of  field  lines  for  the  solution  described  by  eq.(8)  with  2  =  1. 
All  the  points  lying  on  a  given  field  line  define  a  family  whose  locus  is  the  surface  of  a 
torus.  If  any  field  point  approaches  the  radial  distance  where  ji(ar)  =  0,  the  other  points 
in  its  family  lie  on  a  torus  that  approaches  a  spherical  overall  shape,  as  in  fig.(4).  In  the 
limit,  the  field  lines  on  the  spherical  part  of  such  surfaces  approach  the  laminar  velocity 
field  lines  of  an  incompressible  fluid  flowing  over  a  rigid  sphere,  with  the  poles  being  the 
stagnation  points*7* . 

The  2  =  1  case  will  be  very  important  in  the  rest  of  this  work,  since  thi?  case  has  the 
highest  symmetries  and  the  smallest  field-line  curvatures  for  a  given  radius  where  the  first 
zero  of  ji(ar)  occurs.  An  2  =  1  spherical  plasmoid  will  henceforth  be  assumed  to  have  a 
minimum  radius  Rq  =  4.493 . . .  /a.  The  radius  Rq  corresponds  to  the  first  zero  of  j\(ar). 
At  any  radius  corresponding  to  a  zero  of  ji(ar),  the  force-free  vector  field  in  eq.(8)  has  no 
radial  component.  The  importance  of  this  property  will  become  evident  below. 

Analytic  force-free  fields  in  toroidal  coordinates  have  also  been  investigated.*'* 

The  name  “force-free”  came  from  the  observation  that  a  D.C.  magnetic  field,  under 
vacuum  conditions  and  satisfying  eq.(l),  would  have  to  arise  from  a  volumetric  current 
field  that  was  parallel  to  the  magnetic  field  at  every  point.  For  such  a  magnetic  field,  the 
stresses  on  the  currents,  which  come  from  the  Lorentz  force  term  poJxH,  where  J  is  the 
current  density,  H  is  the  magnetic  field,  and  po  is  the  permeability  of  the  vacuum,  would 
be  nil. 
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Figure  2.  Ergodic  flow  line*  for  the  vector  field  in 
eq.(8),  with  asl  and  l  =  1.  The  starting  point 
for  the  flow  lines  is  on  the  equator,  at  r  =  4.49. 
The  figure  traced  out  by  the  lines  has  the  topology 
of  a  torus,  having  a  spherical  “outer”  shape  with 
a  “hole”  along  the  symmetry  axis.  If  the  radius  of 
the  starting  point  was  exactly  at  the  first  sero  of 
;i(or), at  r  =  4.493 .,  the  outer  shape  would 
exactly  be  spherical  and  the  hole  auld  have  sero 
volume.  This  figure  is  viewed  orthogonally  to  ttie 
symmetry  axis. 


Figure  3.  Ergodic  fiow  lines  for  the  vector  field 
in  eq.(8),  with  a  =  l  and  /  =  1.  The  starting 
point  for  the  flow  lines  is  on  the  equator,  at  r  s 
3-37,  about  three-quarters  of  the  distance  from 
the  origin  to  the  first  zero  of  ji{ar).  The  surface 
supporting  the  fiow  lines  is  donut-shaped.  This 
figure  is  viewed  at  an  angle  of  45*  to  the  symmetry 
axis. 


Figure  4.  Ergodic  flow  lines  for  the  vector  field  in 
eq.(8),  with  a  ss  1  and  /  =  1.  The  starting  point 
for  the  fiow  lines  is  on  the  equator,  at  r  =  2.8, 
about  three-quarters  of  the  distance  from  the  ori¬ 
gin  to  the  first  zero  of  ii(or).  The  surface  sup¬ 
porting  the  flow  lines  resembles  a  thin  donut.  This 
figure  is  viewed  at  an  angle  of  45*  to  the  symmetry 
axis. 
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Since  the  stresses  in  such  a  system  would  be  minimized,  it  follows  that  the  energy 
directly  associated  with  the  magnetic  field  must  be  at  an  extremum.  The  attempt  is  made 
to  utilize  this  principle  in  the  design  of  compact  toroids  and  spheromaW8*.  In  these  devices 
an  effort  is  directed  to  produce  force-free  plasma  regions  that  will  have  long  lifetimes  and 
thereby  be  hospitable  for  thermonuclear  reactions. 

Force-free  fields  were  investigated  in  the  nineteenth  century  by  Beltrami^10*,  and  are 
still  under  active  study^11'.  Although  the  interested  reader  will  find  an  ample  literature  on 
the  subject,  the  non-plane-wave,  time-harmonic  case,  including  self-consistent  current, 
has  not  been  published. 

III.  TIME-HARMONIC  FORCE-FREE  ELECTROMAGNETIC  FIELDS 
Assume  the  electric  and  magnetic  fields  in  a  neutral  plasma  can  be  written 

E(r,  t)  =  E(r)  sin  (ut  +  ip),  (9) 

and 


H(r,  t)  =  H(r)  cos  (wt  +  ip),  (10) 

where  w  is  the  angular  frequency  of  the  time-harmonic  electromagnetic  field  and  <p  is  a 
phase  angle.  Note  that  the  fields  are  described  by  real  numbers  in  eqs.(9)  and  (10). 

Maxwell’s  equations  in  a  plasma  can  be  written 

V  xE(r)  =  pwH(r)  =  Z0kH(r),  (11) 

and 

VxH(r)  =  ewE(r)  +  J(r)  =  |-E(r)  +  J(r),  (12) 

where  p  is  the  permeability  of  the  plasma,  e  is  the  permittivity  of  the  plasma  (both  assumed 
spatially  invariant),  k  is  the  wavenumber,  k  =  ,/epw,  and  Zo  is  the  wave  impedance, 
Zo  =  \fpft  «  377.6  ohms  in  vacuo. 
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If  the  magnetic  field,  H(r),  and  the  spatial  part  of  the  current,  J(r)  are  assumed  to  be 
force-free,  as  defined  by  eq.(l),  then  it  is  easy  to  see,  from  eq.(ll)  that  the  electric  field, 
E(r)  will  also  have  this  property.  From  eqs.(ll)  and  (12) 

E(r)  =  Zo-H(r),  (13) 

a 

and 

JW-.(.-§)H(r).£(.-g)«r).  (H) 

If  no  free  charge  exists  associated  with  the  frequency  u>,  then  it  can  be  shown  that  it 
is  merely  sufficient  for  H(r)  to  have  the  force-free  form  for  all  the  other  electromagnetic 
quantities  to  have  this  formal .  The  reason  for  this  is  that  under  such  conditions,  V-E  =  0, 
and  the  electric  field  must  be  derived  entirely  from  the  curl  of  a  vector  field.  By  introducing 
the  vector  potential,  A  (see  eqs.(17-19)),  E(r)  can  be  shown  to  be  entirely  derived  from 
the  curl  of  H(r),  and  must  therefore  have  the  force-free  form  if  H(r)  is  in  this  form.  It 
then  follows  from  eq.(12)  that  J(r)  must  also  have  the  force-free  form. 

Prom  eq.(12)  it  is  easy  to  see  that  the  current  and  the  magnetic  field  have  the  same 
time  dependence  factor,  cos(<ji  4-  <p).  There  is  a  critical  frequency,  u0,  for  a  fixed  value  of 
a,  where 

.  w0  =  -£=  =  ac,  (15) 

and  c  is  the  speed  of  light  in  the  plasma.  At  tins  critical  frequency  the  phase  of  the  current 
changes  sign  with  respect  to  the  electric  and  magnetic  fields. 

The  significance  of  the  last  statement  can  be  appreciated  by  considering  the  energy, 
Q,  associated  with  a  force-free  field,  from  first  principles, 

Q  =  |  J  E-E dv  + 1 J  H-Hdu.  (16) 

Both  E  and  H  can  be  expressed  in  terms  of  the  vector  potential,  A,  where 

A(r,  t)  =  — H(r,  t)  =  A(r)  cos(u)t  +  <p),  (17) 

Qt 
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(18) 


E(r,t)  =  - 


dA(r  ,t) 
dt  ’ 


and 


H(r,t)  =  ivxA(r,t). 


(19) 


Equation  (16)  now  V'omes 

Q  =  |  J  +  \J  HVx  A dv.  (20) 

Prom  the  vector  identity 


V-(AxH)  =  H-VxA  —  A-VxH, 


(21) 


the  fact  that  AxH  =  0,  and  eq.(12), 


Q  =  Y  ain2(wt  +  (p)  J  E(r)- A(r )dv  + 1  cos 2(wt  +  <p)  f  Mr)'  [«“E(r)  +  J(r)]dv.  (22) 
Combining  terms, 

Q  =  y  J  TJ(r)-A(r)du  +  i  coe2(wt  +  <fi)  J  A(r)-J(r)du.  (23) 

The  first  term  in  eq.(23)  is  independent  of  time,  but  the  second  term  is  not.  TUdng 
the  time  average,  <5, 

C  =  y/ E(r)-A(r)du  +  ^J  A(r)-J(r)dv.  (24) 

The  first  term  in  eq.(24)  is  always  non-negative.  The  second  term  is  positive  for 
u  <uo  and  negative  for  u  >wq.  To  see  this,  substituting  for  J(r)  and  A(r), 
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(25) 


/  A(r)-J(r)du  =  (l  -  ^)M/ |H(r)|2du. 

Not*  that  the  integral  is  always  positive,  but  the  factor  multiplying  the  integral  goes 
through  zero  and  changes  sign  at  the  critical  frequency. 

The  second  integral  in  eq.(24)  is  the  por¬ 
tion  of  the  energy  in  the  electromagnetic  field 
arising  from  the  direct  interaction  of  the  cur¬ 
rent  with  the  magnetic  field.  Since  Lorentz 
forces  on  the  current  are  zero,  this  part  of  the 
energy  is  at  an  extremum.  This  extremum 
represents  a  minimum.  It  can  be  shown  that 
the  kinetic  energy  density  of  the  electrons  is 
equal  and  opposite  in  sign  to  this  term,  (see 
Appendix  D.) 

An  interesting  consequence  of  the  phase 
change  above  a-  critical  frequency  is  demon¬ 
strated  using  fig.(5).  figure  (5)  shows  the 
vector  potential.  A,  and  two  currents,  J  and 
J  +  4J,  all  at  the  same  point.  The  current 
J  is  the  current  that  would  naturally  arise 
in  a  self-consistent  force-free  time-harmonic 
field.  Note  that  1  is  exactly  out  of  phase  with 
A.  The  current  J+ S3  represents  the  original 
current  perturbed  in  such  a  way  that  the  two 
magnitudes  are  equal,  i.e., 

|J|  =  |J  +  «J|.  (26) 

A  perturbation  such  as  this  may  arise  from  a  local  statistical  fluctuation  of  the  plasma, 
redirecting  the  current  but  not  changing  its  magnitude.  Since  the  unperturbed  current  and 
the  vector  potential  are  180s  out  of  phase,  and  since  the  vector  potential  is  not  affected, 
to  lowest  order,  by  a  small  local  current  change, 

(J  +  M)-A>J-A,  (27) 

Tins  implies  that  a  small  local  deflection  of  the  current  would  increase  the  energy  of  the 
system,  so  that  it  behaves  like  an  elastic  body.  The  same  logic  predicts  the  system  may 


figure  5.  Perturbation  of  the  J  field.  In  a  force- 
free  plasmoid  above  the  critical  frequency,  a  small 
perturbation,  4J,  of  the  current,  J,  will  add  energy 
to  the  system. 
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be  unstable  below  the  critical  frequency  because  a  small  local  perturbation  of  the  current 
lowers  the  energy. 


As  will  be  seen  in  the  next  sections,  stable  time-harmonic  force-free  plasmoids  may 
not  exist  below  the  critical  frequency,  from  several  considerations.  The  above  arguments 
will  be  revisited  when  fluid  velocities  and  stresses  are  introduced. 

IV.  RADIATION  AND  IDEAL  PLASMOIDS 

Before  discussing  electronic  charge  carriers  and  the  conditions  for  overall  stability,  the 
question  of  radiation  can  be  addressed.  The  assumption  to  be  made  is  that  plasmoids  will 
form  in  such  a  manner  that  the  electric  and  magnetic  fields  will  be  tangent  to  the  surface 
boundary.  The  rationale  for  this  will  become  clearer  when  fluid  velocities  are  introduced.  A 
self-consistent,  stable  plasmoid  cannot  have  a  component  of  velocity  normal  to  its  surface. 
Since  the  velocity  field  to  be  introduced  will  be  identical  in  form  to  the  electric  and  magnetic 
fields,  the  conclusion  is  reached  that  the  electromagnetic  fields  must  be  tangent  to  the 
plasmoid  surface.  Only  spherical  plasmoids  will  be  discussed,  but  the  results  are  easily 
extended  to  cover  possible  non-spherical  cases. 

Instead  of  going  into  extensive  mathematical  proofs,  which  is  done  elsewhere'13* ,  and 
in  Appendix  A  for  the  1  =  1  case,  analogy  will  be  made  to  the  field  inside  a  perfectly 
conducting  spherical  cavity  excited  at  a  transverse-electric  (TE)  resonance  frequency.  In  a 
TE  resonance  the  electric  vector  has  no  radial  component  and  must  therefore  be  zero  at  the 
surface  of  the  cavity.  The  condition  for  resonance  is  the  vanishing  of  the  total  electric  field 
and  the  normal  component  of  the  magnetic  field  at  the  cavity  surface.  These  conditions 
are  met  by  the  wavenumber,  it,  and  the  cavity  radius,  a  satisfying!** 


ji(ka)  =  0, 


(28) 


where  ji  is  the  usual  1-th  order  spherical  Bessel  function. 

The  source  for  the  cavity  field  is  the  surface  current,  nxH,  at  the  boundary,  where 
n  is  the  inward-pointing  normal  and  H  is  the  surface  magnetic  field.  This  surface  current 
produces  no  field  outside  the  cavity. 

The  analogy  to  the  force-free  plasmoid  field  can  now  be  made.  If  #  in  eq.(4)  is  taken 
to  be  other  of  the  functions  (“o*  for  odd  or  “e”  for  even  parity  with  respect  to  the  polar 
angle  <j>) 


=  j«(kr)P^(coetf)  m^, 
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(29) 
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then  eqs.(5),  (6)  and  (7)  give  vector  fields^* 

m;,,m  =  (30) 

and 

njl.m  = 

(31) 

For  a  TE  resonance,  the  electric  field  can  be  written  (“oc”  means  “proportional  to”) 
Eji.m  «  m  ji,mI  (32) 

and  the  magnetic  field  can  be  written 

H;i,m  «  nji,m.  (33) 

The  force-free  field  vector  derived  from  can  be  written 

F:i.m  =  n»:»,  m  +  (34) 

with  a  written  for  k.  If  the  radius  a  of  a  force-free  plasmoid  is  chosen  to  satisfy  eq.(28), 
then  both  the  electric  and  magnetic  fields  at  the  surface  will  have  the  form  given  by  eq.(30), 
because  m«(  m  =Gif  eq.(28)  is  true.  Since  the  field  of  the  resonant  cavity  does  not  radiate, 
neither  of  the  fields  of  the  plasmoid  can  radiate. 

The  implication  is,  then,  that  by  choosing  a  radius,  Ho,  that  simultaneously  satisfies 

ji(oHo)  =  0  and  ;i(kHo)  =  0  (36) 

drastically  reduces  the  radiation  from  the  surface  of  a  plasmoid. 
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The  analogy  between  the  cavity  and  the  plaamoid  is  shown  in  fig.(6).  From  now  on, 
it  will  always  be  assumed  that  the  plasmoids  obey  eqs.(35). 


SPHERICAL  CAVITY 

PIASMOID 


Figure  6.  Plasmoid  ss.  cavity  fields.  At  the  sur¬ 
face  of  a  perfectly  conducting  cavity,  the  electric 
field  of  a  TM  mode  is  tero,  with  the  correspond¬ 
ing  magnetic  field  spatially  identical  in  form  to  the 
plaamoid  surface  fields. 


V.  ELECTRONIC  CURRENTS 

In  a  real  plaamoid  high-frequency  cur¬ 
rents  are  assumed  to  be  carried  by  the  elec¬ 
trons.  In  what  follows,  the  electrons  will  be 
treated  dice  a  gas,  neutralized  by  the  back- 
gound  of  ion3,  as  shown  in  fig.(7).  Difficul¬ 
ties  will  be  ignored  that  might  be  encoun¬ 
tered  when  the  wavelength  of  the  radiation 
becomes  so  small  that  it  approaches  the  size 
of  the  Debye  length.  Likewise,  statistical  ef¬ 
fects  due  to  temperature  will  not  be  consid¬ 
ered.  In  the  absence  of  fluid  motion,  the  over¬ 
whelmingly  dominant  forces  assumed  acting 
on  the  electrons  will  come  from  their  inter¬ 
action  with  the  electromagnetic  field.  Under 
these  conditions  the  ions  will  represent  a  sta- 


Figure  T.  The  seutrsl  plums.  The  electron  den¬ 
sity  in  s  plums  will  be  suumed  high  enough  to  be 
approximated  u  s  gu  in  ter  sc  ting  with  the  neu¬ 
tralising  ions  and  the  neutral  atoms.  In  a  plaa¬ 
moid,  above  the  critical  frequency,  time-harmonic 
current  is  carried  by  the  electrons. 
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tionary  background  and  will  also  be  ignored 
in  this  section.  The  effects  of  fluid  motion 
wilt  be  investigated  in  a  subsequent  section. 

Electronic  motion  is  governed  by  Lorentz 
forces,  demonstrated  in  fig.(8).  In  an  electro¬ 
magnetic  field  the  force  on  an  electron  is 


m,r  =  eE  +  e/iVxH,  (36) 


where  m,  is  the  electronic  mass,  f  is  the  accel¬ 
eration,  v  is  the  electron  velocity,  and  e  is  the 
electronic  charge.  If  electrons  could  exactly 
follow  curved  electric  field  lines  the  second 
term  in  eq.(36)  would  be  zero  in  force-free 
fields. 


figure  S.  Force*  on  »  moving  electron.  Electnc 
forces  act  on  the  electrons  whether  or  not  they  are 
in  motion.  Magnetic  force*  act  only  on  moving 
eiect-ona. 


As  a  first  approximation,  electrons  will  be  assumed  to  execute  very-small-amplitude 
simple  harmonic  motion  tangent  to  the  field  lines,  with  the  v  x  H  term  ignored.  In  that 
case,  mring  eq.(9),  the  first  integral  of  eq.(36)  gives  the  velocity  at  ro. 


v  =  f  =  -  (—)  cos(wt  +  <p).  (37) 

\m,/  u 

This  velocity  can  now  be  used  to  determine  the  free— electron  charge  density,  p,.  Since 
J  =  p,y,  eq.(14)  and  (37)  can  be  combined  to  give 


(38) 


Equation  (38)  says  two  things  about  pt.  It  calculates  the  magnitude  and  the  lign  of  the 
electronic  charge  density.  Since  the  sign  of  the  charge  density  must  be  negative,  because 
the  sign  of  the  electronic  charge,  e,  is  negative,  it  follows  that  electrons  cannot  carry  the 
high-frequency  current  in  force-free  plasm oids  below  the  critical  frequency.  That  is,  if 
k  <  a,  then  the  current  must  be  carried  by  the  ions.  Since  the  ionic  mass  is  so  much  larger 
than  electronic  mass,  this  is  clearly  an  impossibility,  at  least  from  a  classical  viewpoint. 


There  are  several  conclusions  to  be  drawn  from  eq.(38).  The  first  is  that  self-consistent 
force-free  plasmoids  cannot  form  below  the  critical  frequency.  When  they  can  form,  the 
free-elect ron  charge  density  must  be  uniform,  and  increases  quadraticaily  with  frequency. 


Although  the  above  analysis  directly  gave  the  charge  density  by  assuming  the  current 
was  consistent  with  Maxwell’s  equations,  and  electrons  obeyed  Newton’s  laws  of  motion,  no 
mention  is  made  of  stability.  A  full  analysis  of  stability  will  have  to  include  the  effects  of  the 
electromagnetic  stresses  internal  to  the  plasmoid.  ft  will  be  shown  later  that  these  stresses 
can,  for  the  most  part,  be  replaced  by  a  pressure  proportional  to  the  square  of  the  electric 
field,  which  can  be  cancelled  by  vortical  motion,  in  a  non-relativistic  approximation.  The 
simplified  analysis  to  be  given  below  does  not  Include  electromagnetic  stresses,  or  can  be 
thought  of  as  assuming  that  they  have  been  cancelled  by  other  mechanisms. 

Under  this  assumption,  it  will  be  shown  that  stability  is  still  not  assured  unless  the 
frequency  and  charge  density  are  identical  to  those  that  have  already  been  derived  from 
other  considerations.  It  will  also  be  shown  that  small  perturbations  in  the  magnitude  of 
the  free-electron  charge  density  will  cause  the  electron  gas  to  contract  if  the  magnitude 
increases  and  expand  if  the  magnitude  decreases. 

VL  VHUAL  CONSIDERATIONS 

Consider  the  moment  of  inertia,  I,  of  the  electrons  in  a  finite  plasmoid.  and  its  time 
derivatives: 

t=  [  pm|r|*dV,  (39) 

Jv 

~  =  i  =  2jvPmr-TdV,  (40) 

and 


^|  =  /  =  2^pm[r-f+|f|l)dV,  (41) 

where  pm  is  the  (assumed  spatially  uniform)  mass  density  of  the  electrons,  r  is  the  position 
vector  associated  with  an  “infinitesimal”  blob  of  electrons,  and  the  integration  is  over  the 
volume  of  the  plasmoid.  Strictly  speaking,  the  integrations  in  eqs.(39)-{41)  should  be  sums 
over  individual  electrons,  but  their  densities  are  assumed  high  enough  for  the  integral  to  be 
an  appropriate  approximation.  The  quantity  /  is  will  be  called  the  “virial  of  the  electronic 
motion”*14'15'1*'17). 

The  exclusion  of  the  electromagnetic  field  densities  from  the  above  integrals  is  assumed 
to  be  justified,  in  the  non-relativistic  limit,  to  deal  with  lowest-order  effects  arising  from 
variations  of  charge  density  and  current.  Subsequent  analysis  will  deal  with  the  direct 
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effects  of  field  densities  on  the  mechanical  stresses  in  a  plasmoid,  and  how  they  interact 
with  the  fluid  dynamics  of  the  plasma.  It  is  assumed  that,  under  non-relativistic  conditions, 
these  separate  effects  can  be  dealt  with  independently. 

The  quantity  r  =  r(x,y,z,t)  is  the  instantaneous  position  of  an  infinitesimal  volume 
whose  mean  position  is  at  the  fixed  point  ro  =  r(z,  y,z,t  =  to),  and,  to  first  order,  merely 
oscillates  with  the  same  frequency  as  the  applied  time-harmonic  field.  By  definition,  the 
velocity,  v,  is  proportional  to  J: 


v  =  J/pe,  (42) 

where  p,  is  the  (spatially  uniform)  charge  density  of  free  electrons.  The  charge  den¬ 
sity,  p«,  is  simply  related  to  the  mass  density,  pm  fay  />«  =  (e/m)pm,  where  (e/m)  = 
— 1.759xlOn  coul/kg.  Since  all  the  quantities  are  assumed  to  be  harmonic  in  time,  so  that 
J(x,  y,  z,  t)  s  J(z,  v,  z)  cos  (art  +  <p),  then  r  can  be  written,  for  small  oscillations  about  the 
point  ro, 

r  =  r„  +  f  cos (ut  +  <p)dt  =  to  +  sin(u it  +  <p).  (43) 

P«  J  Ptu 

A  necessary  and  sufficient  condition  for  stability  is  that  the  time-integrated  average 
of  /  is  zero.  This  condition  implies  that  the  moment  of  inertia  of  the  electron  cloud  can 
assume  its  original  value  any  number  of  times,  and  cannot  become  uniformly  unbounded  as 
time  increases.  Instead  of  working  directly  with  I,  however,  investigation  will  be  made  into 
conditions  under  which  the  time-integrated  average  of  /  is  zero.  This  condition  is  more 
subtle.  A  necessary,  but  cot  always  sufficient,  condition  for  the  time-integrated  average  of  I 
to  be  zero  is  that  the  time-integrated  average  of  /  be  zero.  What  will  be  shown  is  that  if  the 
uniform  charge  density  has  exactly  the  right  size,  and  the  frequency  is  above  a  critical  value, 
the  time-integrated  average  of  f  is  exactly  zero  if  all  other  stresses  are  balanced.  A  larger 
(smaller)  value  for  the  charge  density  would  drive  the  time-integrated  average  of  I  negative 
(positive),  leading  to  contraction  (expansion)  and  ultimate  instabilities.  This  is  because 
the  exact  value  of  uniform  charge  density  derived  for  a  null  is  the  same  required  from 
classical  mechanics  to  simultaneously  satisfy  Maxwell’s  equations  and  the  laws  of  inertia, 
with  finite  rr.sss  electrons  in  harmonic  motion  along  field  lines.  Since  other  factors  will 
enter  when  considering  the  electromagnetic  field  densities,  the  conditions  to  be  derived  will 
for  the  charge  densities  must  only  be  considered  necessary,  but  not  sufficient  for  stability. 
Gerjouy  and  Stabler*1*!  discuss  a  similar  point. 

Substituting  the  time-dependent  r  from  eq.(43)  into  eq.(41)  and  taking  the  time 
average  automatically  gives  zero,  because  of  the  assumption  that  the  charges  oscillate  with 
fixed  centers  implies  confinement.  It  is  far  more  interesting  to  substitute  for  the  terms 
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in  J,  writing  field  quantities  in  terms  of  the  current,  to  calculate  conditions  for  stability. 
Denote  the  time-integrated  average  of  a  quantity  by  a  bar,  i.e., 


X  =  lim 
r— oo 


i: 


Xdt 


(44) 


To  obtain  new  information  from  oqs.(40),  (41)  note  that  the  term?  in  the  integrand 
containing  ?,  the  acceleration,  are  related  to  the  volume  force  on  the  electrons^1*), 


Pm?  =  PeE(r)  COs(wt  +  if). 
The  term  with  |r|5  can  be  written 


(45) 


(46) 


IVom  eqs.(14),  (40),  (41),  (43),  (44),  (45)  and  (46),  with  a  little  algebra, 

[S+ <«> 


hold: 


In  order  for  the  right  hand  side  of  eq.(47)  to  be  zero  the  following  conditions  must 


k>  a, 


(48) 


and 


(49) 


where  (m/e)  =  -5.69x  10~ukg/coul  is  the  mass-to-charge  ratio  for  the  electron.  Equa¬ 
tion  (49)  can  be  derived  using  a  Lagrangian  formulation.  Prom  variational  considerations 
introduced  in  classical  mechanics,  the  result  in  eq.(49)  can  be  shown  to  arise  from  only 
assuming  that  the  quantities  E,  H,  and  J,  have  the  force-free  form,  and  that  the  free  elec¬ 
trons  carrying  the  current  have  non-zero  mass.  Uniform  charge  density  is  not  an  a-priori 
assumption  in  the  Lagrangian  calculation^15*. 


16 


Besides  providing  the  criterion  necessary 
for  stability,  eq.(47)  predicts  the  consequences 
of  the  magnitude  of  pt  deviating  from  that  in 
eq.(49).  If  \pt\  increases  or  decreases  in  size, 
the  bracketted  factor  in  eq.(47)  will  decrease 
or  increase,  respectively.  If  this  occurs  in  iso¬ 
lated  regions,  the  bracket,  naturally,  cannot 
be  taken  out  of  the  integral.  Physically,  the 
situation  is  diagrammed  in  fig.(9).  A  local 
increase  in  \p,\  means  that  the  contribution 
to  (PI/dt2  will  become  negative  in  the  neigh¬ 
borhood  of  the  increase,  leading  to  the  con¬ 
tribution  to  dll  it  becoming  negative,  lead¬ 
ing  to  a  local  contraction  of  the  electron  gas. 

This  is  an  unstable  situation,  because  it  will 
cause  a  deviation  from  force-free  motion  and 
could  destroy  the  plasmoid.  Similar  reason¬ 
ing  leads  to  the  conclusion  that  a  decrease  in 
charge  density  leads  to  a  net  expansion  of  the 
plasmoid. 

Physically  this  can  be  understood  by  considering  the  energy  of  self-interaction  of  the 
current  in  a  plasmoid.  Two  ride-by-side  line  currents,  of  equal  size,  in  close  proximity, 
attract  each  other  with  a  force  proportional  to  the  square  of  their  current  because  the 
magnetic  field  between  them  has  a  smaller  energy  density  than  in  the  region  outside. 
The  magnetic  “pressure”  arising  from  the  electromagetic  stress  tensor(!S)  pushes  the  wires 
together.  This  is  also  the  origin  of  the  “pinch  effect”  in  plasma  physics.  Increasing  the 
free-elect ron  charge  in  a  plasmoid  will  increase  the  current,  leading,  qualitatively,  to  a 
volume  pinch  effect  tending  to  contract  the  volume  of  the  plasmoid.  Similarly,  decreasing 
the  charge  density  decreases  the  pinch  effect  and  leads  to  expansion. 

This  (increase/ contraction)-^  decrease/expanrion)  seems  to  be  an  inherent  instability 
of  a  plasmoid  where  statistical  fluctuations  of  charge  density  might  occur.  Any  local  in¬ 
crease  in  charge  density,  causing  a  contraction,  would  tend  to  compress  the  region,  making 
the  free-elect  ron  charge  density  increase  even  more.  Similarly  for  the  decrease /expansion 
case.  This  positive  feedback  is  countered  by  the  tendency  of  an  ionized  plasma  in  a  strong 
electromagnetic  field  to  have  its  radiative  recombination  rate  increase  with  density  and  its 
ionization  rate  to  decrease  with  density.  Statistical  fluctuations  in  charge  density  would 
then  be  nullified  to  some  extent  by  ionization/recombination  processes  in  real  gasses.  A 
full  calculation  investigating  this  phenomenon  has  not  been  done  as  yet.  The  sensitivity  of 
the  ionization/recombination  rates  to  gas  density,  as  a  function  of  frequency  and  electric 
field  amplitude  may  be  a  crucial  factor  in  plasmoid  stability. 

As  was  already  mentioned,  the  conditions  in  eqs.(48)  and  (49)  are  necessary  for  stabil- 


Figure  9.  Effects  of  local  charge  perturbations. 
Local  fluctuations  in  electron  charge  density  can 
cause  plasmoids  to  constrict  or  expand. 
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ity,  if  all  other  stresses  are  absent.  The  added  condition  that  makes  (48)  and  (49)  sufficient 
is  that  dl/dt  =  0,  since  equations  (48)  and  (49)  are  the  conditions  also  derived  from  solv¬ 
ing  Newton's  laws  of  motion  for  the  plasmoid  system  when  the  electromagnet  stresses 
are  ignored.  Since  the  fields  are  stationary  in  space  and  time,  except  for  the  harmonically 
oscillatory  part,  it  follows  that  at  t  =  0  the  plasmoid  initially  has  dl/dt  =  0,  where  the  T 
in  eq.(44)  is  time  for  a  single  cycle  of  the  oscillatory  field,  T  =  2 it/u. 

Ignoring  electromagnetic,  thermodynamic  and  fluid  stresses,  eqs.(48)  and  (49),  con¬ 
straining  the  frequency  to  be  above  a  critical  minimum  and  the  free-electron  charge  density 
to  be  uniform  and  a  specific  quadratic  function  of  k  and  a,  represent  necessary  and  suffi¬ 
cient  conditions  for  plasmoid  stability. 

In  the  next  section,  electromagetic  stresses  will  be  discussed  in  more  derail.  Unless 
these  stresses  are  balanced,  they  will  tend  to  make  the  virial  of  the  electronic  motion,  /, 
positive.  In  other  words,  the  electromagnetic  field  trapped  in  a  plasmoid  acts  like  a  positive 
overpressure  tending  to  make  the  plasmoid  expand  and  dissipate.  A  mechanism  will  subse¬ 
quently  be  discussed  that  will  balance  out  all  but  a  very  small  part  of  the  electromagnetic 
stresses. 

VII.  ELECTROMAGNETIC  STRESSES 

The  usual  analysis  involving  the  virial  theorem<14,15,l*,I7)  includes  the  electromagnetic 
energy  density  in  the  integrand  in  eq.(39).  The  assumptions  made  here  are  that  the 
electronic  velocities  are  non-relativistic  and  the  harmonic  fields  E  and  H  are  non-zero 
inside  the  plasmoid  and  do  not  have  the  force-free  form  outside.  It  is  possible  to  have 
an  external  time-harmonic  field  that  obeys  free-space  conditions  with  no  discon. inui ties 
appearing  at  the  plasmoid  surface.  The  external  branch  of  the  electromagnetic  field  has 
no  currents,  does  not  have  the  force-free  form,  and  will  be  assumed  to  have  no  stresses 
associated  with  it. 

It  is  important  now  to  estimate  how  close  a  real  physical  system  can  be  made  to 
approach  the  conditions  for  force-free-field  approximations  on  the  vector  fields  E,  H,  and 
J.  If  electrons  were  massless,  the  problem  would  be  far  simpler,  since  they  would  experience 
no  centripedal  accelerations  and  would  conform  to  electric  field  lines  like  currents  in  wires. 
In  reality,  the  electrons  cannot  confine  their  motion  to  be  strictly  along  the  curved  field 
lines  of  a  force-free  field. 

The  assumption  that  the  electrons  undergo  ample  harmonic  motion  along  electric  field 
lines  becomes  more  accurate  as  the  frequency  increases  and  the  electric  field  decreases.This 
can  be  seen  from  eq.(43),  where  the  amplitude  of  the  harmonic  motion  is  J/(ptu>).  Sub¬ 
stituting  for  J  in  terms  of  E,  a  condition  for  small-amplitude  oscillations  cim  be  written 
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jlp«*o-  (50) 

where  |Elmat  is  the  maximum  absolute  amplitude  of  the  harmonic  electric  field  within  an 
l  =  1  plasmoid  and  Jic  is  the  radius.  Using  R.M.K.S.  units  eq.(5Q)  becomes 

volts/m,  (51) 

where  (  is  an  appropriately  small  factor  that  satisfies  eq.(50),  and  A  is  the  free-space 
wavelength  of  the  harmonic  wave. 

While  eq.(50)  gives  a  criterion  for  small-pmplitude  electronic  motion,  a  more  exact  ex¬ 
pression  can  be  written  for  the  forces  on  individual  electrons  in  a  force-free  time-harmonic 
plasmoid.  When  these  forces  are  known,  trajectories  may  then  be  calculated,  and  an  as¬ 
sessment  can  be  made  of  techniques  for  balancing  the  forces. 

It  will  be. shown  below  that  the  condition  in  eq.(50)  is  related  to  the  problem  of 
electron  “drift”  in  a  plasma.  That  is,  in  this  context,  the  tendency  for  electrons  to  move 
orthogonally  to  electric  and  magnetic  field  lines.  The  approach  will  be  to  expand  the 
acceleration  field  as  a  perturbation  series  about  tbe  harmonic  motion. 

For  purposes  of  perturbative  expansion, 
it  will  be  assumed  that  the  frequencies  to  be 
encountered  will  be  high  enough,  and  the  re¬ 
spective  field  amplitudes  small  enough  to  say 
that  the  major  part  of  the  electronic  velocity 
will  be  harmonic  motion  along  electric  field 
lines  defined  by  solutions  to  eq.(l).  Since  the 
field  lines  have  curvature,  these  oscillations 
will  be  accompanied  by  a  drift  motion  across 
electric  field  lines,  as  in  fig.(10),  unless  the 
effects  of  the  curvature  are  balanced  by  some 
other  mechanism. 

There  axe  several  reasons  why  this  elec¬ 
tronic  drift  is  undesirable:  First,  as  shown  Figure  10.  Drift  due  to  field  -urvsture.  If  the 
above,  the  current  in  an  ideal  force-free  plaa-  field  lines  are  curved,  the  electrons  will  continuslly 
moid  will  not  radiate,  under  the  right  con-  drift, 10  04ci^lte  different  field  lines, 
ditions  of  size  and  frequency.  Drift  currents, 
not  having  the  force-free  form,  axe  capable 
of  radiating,  potentially  unbalancing  the  plas- 


:  =  2.02x107C^ 
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mend  mechanically.  Next,  drift  currents  may  have  a  divergence  that  would  lead  to  a  change 
in  charge  density,  if  the  rest  of  the  plasma  were  not  present.  Since  the  plasma  creates  large 
restoring  forces  to  maintiin  neutrality,  drift  currents  can  be  responsible  for  pressure  gra¬ 
dients  that  could  disrupt  the  symmetry  of  a  plasmoid,  causing  its  self-destruction.  In 
any  case,  drift  currents  may  have  non-linear  effects  that  can  build  in  an  uncompensated 
manner  and  tend  to  shorten  the  lifetime  of  a  plasmoid. 

It  is  important,  then,  to  understand  the  origins  and  dynamics  of  drift  currents.  Toward 
that  end,  the  forces  acting  on  the  free  electrons  will  be  examined  in  detail,  assuming  non- 
relativistic  conditions. 

An  electron  is  acted  upon  by  the  Lorentz  force 


F  =  mer  =  eE  +  erxB,  (52) 

where  B  =  /i oH  and  r  is  the  instantaneous  position  of  the  electron.  Instead  of  an  analysis 
of  eq.(52)  that  will  relate  the  components  of  F  to  elements  of  the  electromagnetic  stress 
tensor120) ,  the  drift  motion  will  be  taken  as  a  perturbation  determined  by  the  oscillatory 
motion  along  electric  field  lines  and  the  curvature  of  the  components  of  the  field  gradients. 
The  connections  to  the  stress  tensor  also  will  be  discussed  below. 

It  will  be  assumed,  now,  that  the  frequency  of  oscillations,  tv,  is  so  high  that  in  the 
time  needed  for  one  full  cycle  of  the  electric  field,  T  =  2ir/u,  the  net  motion  of  the  center 
of  mass  of  an  electron  oscillation  is  very  small  compared  to  the  overall  size  of  the  plasmoid. 
This  slow  net  motion  can  arise  because  of  field  inhomogeneities  in  the  B  and  E  vectors 
as  well  as  the  relative  local  orientation  of  B  and  E.  The  major  electronic  motion,  then, 
over  many  cycles,  will  be  high-frequency  os<  illation  along  local  electric  field  lines.  This 
high-frequency  motion  will  be  superimposed  on  another  component  whose  net  effect  would 
be  to  slowly  transport  electrons  onto  different  electric  field  lines.  The  following  analysis 
closely  follows  the  logic  on  pp.  77-80  of  Gekker’s  book121'. 

Assume  that  the  position  vector  r  may  be  written  as  a  sum  of  an  oscillatory  term, 
ri,  and  and  a  net  transport  term,  r<).  Let  the  oscillatory  term  arise  from  ignoring  the 
magnetic  field’s  contribution  to  eq.(52). 

It  then  follows  from  eqs.  '9)  and  (10)  that 

ri  =  -q«-2  E(r0 ,  t)  cos(wt  +  if),  (53) 

where  jj  =  e/me  =  -1.76xlOucoul/kg,  the  ratio  of  the  electronic  charge  to  its  mass,  and 
the  electric  field  is  calculated  at  position  rg.  The  quantity  ri  takes  into  account  only  the 
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straight-line  motion  tangent  to  a  field  line.  After  the  contribution  from  ro  is  added  on, 
the  trajectory  should  contain  all  the  corrections  for  curvature. 

Equation  (52)  can  be  rewritten: 

fo  +  h  =qE(r0  +  r1,i)  +  »)(fo  +  fj)xB(ro  +  rl)t).  (54) 

Expanding  the  right  hand  side  of  eq.(54)  as  Taylor  series  in  ri,  and  rearranging, 

?o  +  fi  =jj  |E(r0,t)  +  (r1-V)E(r0,t)  +  |(ri-V)*E(r«,f)  + . .  .j  + 

>7*0  x  B(ro,<)  +  (ri-V)B(ro,t)  +  -(ri-V)JB(r0,t)  + . . .  +  (55) 

17*1  x  B(ro,t)  +  (ri-V)B(ro,t)  +  i(ri-V)JB(r0,t)  +  ...  . 

Subtracting  the  second  time  derivative  of  eq.(53),  with  more  rearranging,  the  equation  for 
ro  is 

r0=>7  (n-V)E(ro,t)  +  rixB(ro,t)j  + 

Ji(ri-V)JE(ro,t)  +  ...l  + 

T  J  1  (56) 

>7*0  x  B(r„,  t)  +  (r, ■  V)B(r0, t)  +  - (n •  V)2B(r0, <)  + . . .  + 

>7*i  x  (n  •  V)B(r0,  t)  +  i  (r,  -V) JB(r0,  <)  + . . .  . 

Prom  Maxwell’s  equation  VxE  =  -dB/dt,  eq.(53)  and  the  familiar  vector  identity 
V(A-B)  =  A-VB  +  B-VA  +  AxVB  +  BxVA, 

eq.(56)  becomes 
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iv|aHa1fl|,+ 

2  |  lj  [ 

,[i(r1-V)5E(r0lt)  +  ...]  + 

jjr0x  B(r0,t)  +  (rrV)B(r0,t)  +  j(rrV)?B(r0,t)  +  .. 
ijrix  (r1-V)B(r0,l)  +  i(rt-V)aB(r0,t)  +  ...j. 


(57) 


The  first  term  6n  the  right  in  eq.(57)  can  obviously  lead  to  gross  motion  of  the  center 
of  mass  of  an  electron.  It  says  that  the  electromagnetic  field  can  produce  a  volume  stress 
in  plasmas  that  appears  to  arise  from  a  pressure  proportional  to  the  square  of  the  electric 
field  and  inversely  proportional  to  the  square  of  the  frequency. 

Many  of  the  other  terms  in  eq.(57)  lead  to  oscillatory  motion,  with  no  time-averaged 
velocity  of  the  center  of  mass.  To  lowest  order,  assume  that  the  oscillatory  part  of  r0  is 
negligible.  Then,  since  ri  and  E(ro,  t)  both  have  sin(urf  +  <fi)  dependence,  the  first  higher- 
order  term  to  contribute  a  non-zero  time-averaged  contribution  in  the  first  square  bracket 
is 


(,/6)(rrV)3E(ro  ,t), 


with  sin4(urf +y>)  appearing  as  the  time  dependence.  Similarly,  the  first  higher-order  term 
to  contribute  in  the  last  square  bracket  is 


(r)/2)f1x(rrV)JB(r0,t), 


with  ain2(uit  +  <p)  cosJ(uif  +  ip)  appearing.  To  lowest  order,  none  of  the  terms  in  the 
second  bracket  can  contribute  to  non-osdllatory  motion,  because  B(ro,t)  has  cos(ud  +  y>) 
dependence. 

The  electric  field  must  now  be  restricted  in  size,  not  only  to  make  the  series  in  eq.(57) 
converge,  but  to  prevent  the  higher-order  terms  from  making  a  substantial  contribution 
to  the  motion.  A  physical  restraint  will  be  that  the  ma-rimnm  value  of  ri  be  much  smaller 
than  Eg.  Rom  eq.(53),  define  a  smallness  parameter,  {,  very  similar  to  ( in  eq.(51),  where 


p  _  ’llElnwt 


<1, 


(58) 
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so  that 


|E|mal  =  t~~-  (59) 

Note  that  for  a  fixed  values  of  {  and  (it/a),  u  is  inversely  proportional  to  Ro,  and  the 
pinirimnm  electric  field  intensity  varies  inversely  with  plasmoid  radius.  Maximum  energy 
densities  can  therefore  be  inversely  proportional  to  the  square  of  the  radius  without  vio¬ 
lating  stability  criteria. 

It  can  then  be  shown  that  the  time-averaged  force  on  an  individual  electron  can  be 
written 


m,  ro  =  — -m,V 


qE(ro) 


[l  +  Otf2)]. 


(60) 


In  the  next  section,  it  will  be  shown  that  eq.(60)  produces  an  effective  non-uniform 
fluid  pressure  proportional  to  the  energy  density  in  the  electric  field.  By  introducing 
vortical  motion,  variations  in  the  total  time-averaged  fluid  pressure  governing  the  flow 
within  the  plasmoid  will  be  cancelled  to  order 

VIII.  VORTICAL  PLASMA  STABILIZATION 

The  time-averaged  volume  force  on  the  electron  fluid  at  any  point  is 


(61) 


FYom  eq8.(38)  and  (60), 


Fv  = 


V|E(r)|J 


l  +  0(e) 


(62) 


To  order  this  is  equivalent  to  the  volume  force  that  would  arise  from  the  gradient  of  the 
pressure  field 


p(r)=Ki_^)|E(r)|J+p°’ 


(63) 
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where  Po  is  a  constant  that  may  not  be  zero. 

A  great  deal  of  care  must  be  taken  to  account  for  all  the  stresses  due  to  the  electro¬ 
magnetic  field.  Equation  (57)  is  correct  throughout  a  region  where  all  the  quantities  exist 
and  are  analytic.  It  is  important  to  note  that  the  approximation,  eq.(63),  derived  from 
eq.(57)  is  not  what  can  be  derived  from  a  localized  interpretation  of  the  meaning  of  the 
electromagnetic  stress  tensor  in  a  piasmoid.  This  subtle  point  is  worth  exploring  in  a  little 
mare  depth. 

The  electromagnetic  stress  tensor  in  a  homogeneous,  isotropic  linear  medium,  2S,  is 
physically  defined  by  its  divergence*8’14’20', 


V-(2S)  =  « |(V  xE)  x  E  +  EV-e|  +  ^(VxH)xH  +  HV-H 
which  can  also  be  expressed  in  the  form 


(64) 


V-^SJs^E  +  pJxH,  (65) 

where  pc  is  the  net  free  (unneutralized)  charge  and  J  is  the  current. 

In  an  ideal  force-free  piasmoid,  there  is  no  free  charge,  pc  =  0,  and,  by  definition, 
JxH  =  0.  The  divergence  of  the  electromagnetic  stress  tensor  is  therefore  identically  zero 
there. 

Where  trouble  can  arise  is  in  writing  an  expression  for  the  stress  tensor  having  the 
dyadic  form 

2S  =  eo(EE  -  ~|E|2I)  +  /i0(HH-  |lH|2!),  (66) 


where  I  is  the  unit  dyadic. 

It  is  true  that,  integrating  over  the  volume  and  surface  area  of  a  piasmoid 


J  (peE  +  ^<)JxH)dv  =  £  2S-ndo, 


(67) 


where  h  is  the  outwardly  pointing  normal.  Obviously,  the  integrals  in  eq.(67)  are  equal  to 
zero  for  a  force-free  piasmoid. 
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BVom  eq.(66),  the  integrand  of  the  second  intgrai  in  eq.(67)  am  be  written 


JS-n  =  «  (E-n)E-  i|E|snj  +p  (H-n)H-i|H|2n  .  (68) 


The  “usual”  interpretation  of  the  quantity  in 
eq.(68)  is  that  of  a  local  stress  being  com¬ 
municated  across  the  surface  of  integration, 
see  fig.(ll).  In  this  case  since  both  electric 
and  magnetic  fields  are  have  no  normal  com¬ 
ponents  at  the  surface  of  the  plasmoid,  the 
electromagnetic  stress,  from  eq.(68)  amounts 
to  a  pressure  equal  to  the  total  energy  den¬ 
sity  of  the  electromagnetic  field.  In  reality, 
there  is  no  local  net  farce  in  an  ideal  force- 
free  plasmoid. 

The  analysis  done  to  find  the  pressure 
due  to  nori-zero-mass  electrons  finds  a  quan¬ 
tity,  in  eq.(63),  that  differs  at  the  surface  of 
the  plasmoid  from  assuming  stresses  to  arise 
from  eq.(68).  Even  though,  as  the  frequency 
increases,  with  the  size  of  the  plasmoid  fixed, 
the  differences  become  smaller,  eq.(68)  can¬ 
not  be  used  and  has  no  relevance  to  the  local 
pressure  field  when  the  fields  are  assumed  to  arise  from  the  farce-free  solutions  to  Maxwell’s 
equations.  Relevance  can  be  re-established  by  allowing  plasmoids  to  have  different  fields 
directly-inside  and  outride  their  boundaries,  forcing  currents  to  flow  on  their  surfaces. 

These  currents,  a  single  layer  needed  for  the  jump  in  the  H  field,  and  a  double  layer 
needed  for  a  jump  in  the  E  field,  will  exhibit  a  reaction  to  the  fields  at  a  plasmoid  boundary 
that  will  effectively  act  like  an  outward  load  pressure  proportional  to  the  local  energy 
density  of  the  electromagnetic  field.  Within  a  constant  factor,  this  latter  pressure  has  the 
same  functional  form  as  the  first  term  in  eq.(63). 

As  will  be  shown  below,  the  pressure  represented  by  the  first  term  in  eq.(63)  can  be 
balanced  by  a  vortical  motion  of  the  fluid  in  the  plasmoid.  If  extra  stresses  occur  at  the 
surface,  they  can  be  balanced  by  modifying  the  magnitude  of  the  fluid  velocity  everywhere. 
This  would,  however,  unbalance  the  stresses  everywhere  else  in  the  plasmoid.  In  a  perfectly 
non- viscous  fluid  it  would  also  be  possible  to  modify  the  velocities  found  in  the  outermost 
“toroids”  of  the  vortical  flow  to  compensate  for  the  stresses  surface  currents  would  produce. 
Since  the  flows  in  separate  toroids  of  the  vortical  motion  are  independent,  the  rest  of  the 
plasmoid  would  be  unaffected.  It  is  easy  to  see  that  even  in  the  absence  of  viscosity,  taking 
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Figure  11.  The  strew  tensor.  The  “usual"  inter¬ 
pretation  of  the  stress  tensor  translates  to  a  local 
stress  across  a  surface. 


surface  currents  into  account  complicates  the  theory. 

In  the  next  section  it  is  shown  that  some  or  all  plasmoid  surface  currents  can  be 
eliminated  if  special  external  fields  are  present.  It  will  be  shown  that  these  external  fields 
can  be  described  by  the  resonant  inodes  of  a  perfectly  conducting  spherical  cavity.  In  view 
of  this,  the  theory  will  now  be  developed  ignoring  the  effects  of  surface  currents. 

Even  though  the  gradient  of  the  pressure  field  in  eq.(63)  acts  only  on  the  electron  fluid, 
it  is  transmitted  to  the  ions  via  the  electrostatic  forces  that  enfcrce  quasineutrality  in  the 
plasma.  The  internal  flow  of  the  plasmoid  will  then  be  governed  by  static  and  dynamic 
fluid  pressure  gradients,  in  addition  to  the  forces  in  eq.(62).  In  what  follows,  it  will  be 
assumed  that  the  plasmoid  behaves  fluid— dynamically  like  a  perfect  gas,  with  constant 
thermal  conductivity  and  specific  heats. 

A  flow  field,  with  velocity  distribution  v(r)  will  now  be  introduced  that  counteracts 
the  pressure  in  eq.(63).  Let 


v(r)  =  AM(r), 


(69) 


where  A  is  a  constant  to  be  determined  and  M(r)  obeys  eq.(l).  In  other  words,  the 
streamlines  of  v  are  identical  with  the  harmonic  electric  field  lines.  This  steady  flow  field 
is  potentially  stable,  as  can  be  seen  by  applying  the  equation  for  the  evolution  of  vorticdty, 
<3<7>, 


ss  —  . 

=  Vx(vx<3), 


(70) 


where 


u  =  Vxv  =  av. 


(71) 


In  eq.(70),  the  right-hand  side  is  identically  zero  because  v  has  the  force-face  form, 
with  v  a  multiple  of  Vxv.  Force-free  vortides,  with  vector  field  M(r),  are  not  the  only 
candidates  for  stable  motion.  If  the  vector  function  (vxi!)  is  the  gradient  of  a  scalar 
function,  the  vortex  formed  by  the  velodty  v  is  also  potentially  stable.  Hill’s  spherical 
vortex  is  in  this  latter  class.  Exactly  equivalent  to  eq.(70)  is  the  expression  for  the  rate  of 
change  of  <3  in  an  element  of  fluid  that  is  converted  with  purely  vortical  flow^, 


DCj  _  _ 
—  =w-Vv. 


(72) 


FVom  eq.(71)  this  becomes 
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(73) 


Dy  _ 
—  =  v-W. 


Since  v  has  the  force-free  form,  the  formula  for  the  gradient  of  a  dot  product  of  two  vector 
fields  implies 


Dv 

Dt 


(74) 


Equation  (74)  says  that  the  acceleration  of  an  element  of  fluid  moving  in  a  force-free 
vortex  can  be  derived  from  the  gradient  of  the  square  of  the  velocity.  Multiplying  both 
sides  of  eq.(74)  by  p,  the  fluid  density,  and  assuming  uniform  density, 

<75> 


it  follows  that  the  volume  force  on  an  element  of  fluid  in  a  force-free  vortex  may  be  thought 
to  arise  from  the  pressure  field 


Pv=Po-|p|v|,I  (76) 

where  po  is  a  constant  and  p  is  spatially  invariant.  This  pressure  can  be  adjusted  to  cancel 
the  major  part  of  the  electromagnetic  pressure. 

It  is  immediately  obvious  that  the  value  of  A  in  eq.(69)  should  be 

where  the  spatial  part  of  the  electric  field  is  EoM(r). 

The  above  analysis  leading  to  the  result  that  vortical  motion  can  cancel  the  major 
portion  of  the  electromagnetic  forces  on  the  electronic  fluid,  leaving  only  a  uniform  residual 
pressure,  po,  everywhere  is  so  important  that  it  will  be  derived  in  another  maimer. 

Some  assumptions  will  be  made  first: 

(1)  All  processes  are  adiabatic,  so  that  the  specific  entropy  is  spatially  uniform. 

(2)  The  gas  pressure,  p,  must  equal  the  surrounding  ambient  pressure. 
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(3)  The  gas  density,  p,  must  be  uniform  within  the  plasmoid. 


Any  flows  that  are  incompatible  with  any  of  the  above  assumptions  will  not  be  con¬ 
sidered.  Before  going  ahead  with  the  alternate  derivation,  these  assumptions  should  be 
discussed. 

Assumption  (1)  says  that  all  fluid  motions  were  initiated  without  the  addition  of 
heat.  Also,  no  processes  occured  that  generated  or  transferred  heat.  If  the  plasma  in 
the  plasmoid  was  derived  from  the  surrounding  plasma  then  the  specific  entropy  of  the 
plasmoid  is  equal  to  the  specific  entropy  of  the  surrounding  plasma. 

Assumption  (2)  says  that  the  only  pressure  gradients  that  can  remain  in  the  plasmoid, 
after  the  vortical  motion  cancels  the  electromagnetic  pressure,  must  arise  from  the  ther¬ 
modynamic  gas  pressure.  Asstuning  an  “atmosphere”  surrounding  the  plasmoid  exerts  an 
isotropic  pressure  p„  the  internal  thermodynamic  pressure  must  match  this  for  stability, 
after  surface-current  effects  ate  included.  In  large  systems,  p»  will  probably  arise  from  the 
surrounding  gas  pressure.  In  small  systems,  there  may  be  a  component  of  pressure  from  a 
Casimir  force.  Casimir  forces  will  arise  in  situations  where,  owing  to  the  high  conductivity 
of  the  plasma  due  to  large  numbers  of  free  electrons,  components  of  the  electromagnetic 
zero-point  background  energy  cannot  penetrate  uniformly  into  the  plasma.  This  will  result 
in  a  lower  energy  density  within  the  plasmoid,  and  an  apparent  pressure  gradient  due  to 
the  vacuum.  It  is  beyond  the  scope  of  this  note  to  investigate  zero-point  effects. 

Assumption  (3)  must  hold  in  order  to  allow  the  forces  that  arise  from  vortical  accel¬ 
erations  io  be  derivable  from  a  conservative  field. 

The  key  to  the  alternate  derivation  is  the  understanding  that  the  specific  internal 
energy  of  the  plasma  is  a  sum  of  the  perfect-gas  thermodynamic  specific  internal  energy, 
£,  and  the  electromagnetic  energy  density,  per  unit  mass  of  fluid,  Q.  The  thermodynamic 
quantities  are 


£  =  (c,-c0)T, 


(78) 


where  cf,  c„,  and  T  are  the  specific  heat  at  constant  pressure,  the  specific  heat  at  constant 
volume,  and  the  temperature,  respectively.  fYom  eq.(24)  the  part  of  Q  that  is  affected  by 
the  fluid  flow  is  represented  by  the  term  with  J  in  it,  since  the  current  can  be  transported 
by  the  fluid  flow.  Using  this  logo  an  effective  energy  density  can  be  written 

C.//  =  iXA  =  -^(l-^)|M(r)P,  (79) 


with  A  and  J  defined  in  Sec  I. 
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Bernoulli’s  theorem^7)  now  says  that,  along  every  streamline  the  following  physical 
quantity  is  conserved: 


BB  =  \[v\2  +  Q'ff  +  e+Z. 

Since  Q  is  negative  (it  >  a)  the  logical  choice 


(80) 


v  =  AM(r), 


(69) 


where  A  is  defined  in  eq.(77),  will  exactly  cancel  out  the  electromagnetic  stresses  in  eq.(63), 
a  result  already  derived,  to  within  order  {3.  This  discrepancy  to  order  {3  might  have  beau 
expected,  since  the  higher  order  terms  arising  from  the  curvature  of  the  fields,  in  eq.(57), 
were  ignored. 

In  the  above  proof,  using  Bernoulli’s  theorem,  it  is  assumed  that  the  electronic  current 
is  converted  by  the  fluid  motion.  Since  the  fluid  motion  is  assumed  to  be  a  purely  adenoidal 
mode,  no  compression  takes  place,  and  the  electronic  density  remains  uniform.  Under  these 
conditions  it  could,  alternatively,  be  assumed  that  the  electrons  (and  thereby  the  currents) 
are  not  converted  by  the  fluid,  with  only  the  ions  and  neutral  atoms  participating  in 
the  fluid  flow.  In  this  case  G«//  will  not  appear  in  eq.(80),  but  the  pressure  due  to  the 
electronic  motion  would  still  be  transmitted  to  the  ions  and  would  have  to  be  subtracted 
from  the  pressure  p  in  cq.(80). 

Vi 

It  is  then  easy  to  see  that  exact ly  the  dune  velocity  profile  would  have  to  be  assumed 
to  balance  the  electronic  contribution,  whether  the  electrons  are  converted  with  the  fluid 
or  not.  Rom  simple  energy  considerations,  it  is  more  likely  that  electrons  are  converted 
with  the  fluid  because  there  would  then  be  no  external  static  magnetic  field  due  to  the 
motion  of  the  ions.  (The  Idnetic  energy  of  the  electrons  is  ignored  here.  A  full  analysis 
would  have  to  compare  the  added  kinetic  energy  due  to  electronic  motion  to  the  energy  in 
the  magnetic  field  that  would  arise  if  the  electrons  were  not  entrained.  It  is  assumed  here 
that  the  electronic  kinetic  energy  would  be  smaller.) 

With  the  first  two  terms  in  eq.(80)  cancelling,  the  quantity 


£+Z  =  c,T,  (811 

P 

the  specific  enthalpy,  is  constant  along  any  fiowline.  It  immediately  follows  that  the  tem¬ 
perature  must  be  constant  along  a  flowline.  If  the  motion  of  the  compressible  gas  is  to 
remain  steady,  then  the  pressure  must  be  constant  in  the  plasmoid,  and  assumed  equal  to 
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the  pressure  of  the  “ambient  atmosphere”.  This,  then,  forces  the  mass  density,  p,  to  also 
be  a  constant  within  the  plasmoid.  Since  pressure,  density  and  temperature  are  constant 
throughout  the  plasmoid,  and  may  take  on  the  same  values  as  in  the  ambient  isentropic 
atmosphere,  the  flow  is  isentropic  and,  therefore,  can  be  initiated  by  an  adiabatic  process. 
The  criteria  of  the  flow,  listed  as  assumptions  (1-3),  above,  are  met. 

As  field  strength  increases,  the  fluid  speed,  from  eq.(77),  must  increase  in  the  same 
proportion.  The  question  naturally  arises  as  to  whether  turbulence  will  occur  beyond  a 
threshhold  velocity.  This  cannot  be  answered  in  full,  but  if  the  electrons  are  entrained 
in  the  flow,  there  is  a  mechanism  that  will  resist  small  perturbations  that  may  grow  into 
turbulent  flow. 

Remember  that  the  term  Q  —  (l/2p)J-  A,  which  is  proportional  to  the  time-averaged 
energy  of  interaction  of  the  magnetic  field  and  the  current  in  the  time-harmonic  field,  is 
negative  for  k  >  o^  A  small  perturbation  in  velocity  that  does  not  change  |v|,  but  i  nplies 
a  local  change  of  flow  direction  without  changing  A,  will  decrease  the  magnitude  of  Q. 
The  reason  for  this  is  that  the  entrained  current  will  still  have  the  same  magnitude,  but 
its  component  parallel  to  the  original  direction,  parallel  to  A  must  be  smaller.  A  smaller 
magnitude  of  Q  means  increased  energy,  which  leads  to  a  restoring  force  proportional  to 
|£o|2.  The  fluid  therefore  would  appear  to  become  more  rigid  as  the  field  strength  increases, 
at  least  partially  offsetting  small  statistical  fluctuations  that  could  grow  turbulently.  This 
effect  has  already  been  mentioned  in  association  with  eqs.(26),  (27)  and  fig.(5). 

Note  that  a  force-free  vortex  cannot  occur  by  itself  in  a  compressible  fluid.  The  pres¬ 
sure  gradients  would  immediately  give  rise  to  velocity  components  that  were  not  solenoidal. 
Such  a  situation  would  lead  to  rapid  destruction  of  the  vortical  field.  Similarly,  a  time- 
harmonic  force-free  field  cannot  persist  by  itself  in  a  compressible  fluid  because  of  the 
stresses  due  to  finite-mass  electrons  in  the  presence  of  electromagnetic  fields  with  non-zero 
curvature.  The  remarkable  coincidence  of  force-free  vortices  being  stable  in  the  presence 
of  force-free  electromagnetic  fields  makes  the  entire  theory  viable. 

There  is  a  very  special  case  that  is  worth  mentioning,  even  though  the  possibility  of 
its  physical  existence  may  be  apocryphal.  At  the  critical  frequency,  where  ut  —  ca  the 
current  and  charge  densities  become  zero.  It  follows  that  a  force-free  electromagnetic 
excitation  could  exist  in  otherwise  “empty”  space.  The  energy  would  be  contained  within 
a  spherical  region,  terminated  by  a  set  of  boundary  currents,  and  could  not  radiate  away. 
Classically,  there  is  no  limit  to  the  energy  density  that  might  be  confined  in  such  an 
electromagnetic  mode  if  the  outside  influences  a*e  large  enough  to  contain  the  stresses 
on  any  surface  currents  that  would  be  needed  to  match  boundary  conditions.  If  such  an 
object  could  exist,  the  limit  on  the  field  strengths  would  be  near  the  point  where  the 
vacuum  polarization  due  to  virtual  electron-positron  pair  production  was  affected. 

Classically,  such  a  “vactoid”  (vacuum  plasmoid)  would  be  almost  undetectable.  As¬ 
suming  the  “correct”  terminating  currents,  the  electromagnetic  fields  immediately  outside 
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it*  boundary  are  nil,  and  no  charges  inside  it  are  present  to  radiate.  A  more  modem  view 
says  that  the  trapped  energy  represents  an  E  =  me2  mass,  and  would  have  a  gravitational 
field.  Also,  scattering  from  virtual  electron-positron  pairs  might  eventually  lead  to  decay. 
Unless  there  was  a  direct  contact  with  ionizable  or  charged  particles,  the  presence  of  an 
individual  vactoid  would  go  unnoticed.  The  interaction  of  two  vactiods  or,  for  that  matter, 
two  force-free  plasmoids,  is  not  yet  known,  but  can  be  investigated. 

The  major  argument  against  the  existence  of  vactiods  is  the  apparent  need  for  a  set  of 
surface  currents  to  match  the  inside- to-outside  boundary  conditions,  with  no  mechanism 
to  generate  or  support  them.  The  Lorentz  forces  on  such  currents  would  amount  to  a 
pressure  equal  to  the  energy  density  of  the  electromagnetic  field  at  the  surface.  This 
would  tend  to  force  the  vactoid  to  disperse  unless  there  were  processes  outside  the  vactoid 
that  balanced  the  electromagnetic  pressure  acting  at  the  surface  from  the  inside. 

Even  apart  from  this  vacuum  special  case,  a  reliable  process  for  manufacturing  force- 
free  plasmoids  is  still  to  be  found.  In  the  next  section,  some  of  the  issues  involved  will  be 
discussed.  Appendix  B  discusses  the  case  of  an  energetic  (one  mega  joule),  20  cm.  diameter 
plasmoid,  formed  in  air  at  sea-level  atmospheric  pressure. 

X.  DISCUSSION 

The  analysis  so  far  has  assumed  that  force-free  plasmoids  exist  and  proceeded  to  derive 
some  of  their  properties.  This  was  taken  to  be  a  worthwhile  intellectual  excercise,  since 
the  theory  appears  applicable  to  otherwise  unexplicable  observations.  As  it  stands,  the 
weak  points  in  the  theory  concern  prescriptions  for  the  manufacture  of  force-free  plasmoids 
and  the  nature  of  the  fields  external  to,  yet  compatible  with,  force-free  plasmoids.  These 
last  two  issues  lead  to  the  suspicion  that  describing  autonomous  self-confining  energetic 
plasmas,  such  as  ball  lightning,  with  force-free  fields  is  probably  not  incorrect,  but  may 
be  incomplete.  This  section  aims  at  discussing  the  previous  sections,  highlighting  some  of 
the  shortcomings  of  the  theory. 

Arguments  found  at  the  end  of  Sec.  Ill  accompanying  eqs.(28)-{35)  were  used  to  show 
that  the  fields  at  the  surface  of  an  ideal  plasmoid  do  not  radiate.  No  inferences  were  made 
about  the  transition  in  field  amplitudes  at  the  surface  of  a  plasmoid. 

The  moot  point  a3  to  whether  there  can  be  a  true  discontinuity  in  each,  or  either, 
of  the  electric  and  magnetic  fields,  between  inside  and  outside  the  plasmoid,  will  not  be 
considered  in  great  detail.  If  a  true  discontinuity  can  arise,  the  question  of  how  to  determine 
surface  currents  cannot  be  answered  easily.  The  assumption  will  be  nude  that  if  there  are 
discontinuities  the  transition  to  the  “outside’’,  fields  occurs  continuously. 

That  is,  there  must  be  a  small  region  in  the  neighborhood  of  the  boundary  of  a 
plasmoid  where  the  fields,  along  with  their  fir  i  derivatives,  change  continuously.  If  so,  then 
plasmoids  must  have  currents  in  this  region  to  match  the  boundary  conditions  imposed 
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by  classical  electrodynamics.  (If  not,  the  possibility  that  plasmoids  may  exist  having 
no  surface  currents  must  be  investigated.  Lack  of  surface  currents  leads  to  the  apparently 
illogical  consequence  of  a  standing  wave  having  no  sources.  As  an  approximation,  however, 
this  may  not  be  an  impossible  situation.  There  may  be  waves  whose  space-time  symmetry 
severely  limits  their  ability  to  propagate  in  the  “usual”  manner.  These  waves  do  have 
sources,  but  the  compact  support  of  these  sources  may  be  a  very  complicated  function  of 
the  retarded  fields.  In  the  literature  such  waves  have  been  called  Diffractionless/Phase 
Conjugated  Waves,  Electromagnetic  Missiles,  Bessel  Beams,  etc/22*  A  force-free  plasmoid 
unterminated  by  surface  boundary  currents  may  frill  into  this  category.  The  rate  at  which  a 
plasmoid  disintegrates  if  no  boundary  currents  exist  car.  be  investigated.)  As  has  been  seen, 
these  currents  are  needed  to  satisfy  Maxwell’s  equations,  but  since  they  do  not  radiate, 
they  are  not  needed  to  counteract  a  steady-state  flow  of  energy  from  the  surface  of  the 
inside  region. 

To  make  the  last  point  clearer,  consider  the  case  of  three-dimensional  space  divided 
into  two  regions,  “inside”  and  “outride”  by  a  spherical  boundary.  On  the  inside,  for  r  <  Ho, 
the  electric  and  magnetic  fields  are  Ej  and  Hj.  On  the  outside,  for  r  >  Ho,  the  fields  are 
Eg  and  H<,.  If  n  is  the  outwardly  pointing  normal,  it  is  shown  in  standard  texts'25*  that 
the  inside  and  the  outside  will  evolve  independently  if  the  surface,  at  r  =  Ho  sustains  the 
equivalent  electric  (J)  and  magnetic  (M)  currents 

3  =  -H,),  (82) 

and 


M  =  -nx(E0  —  E,).  (83) 

For  solving  problems  on  computers,  the  carriers  of  these  currents  are  irrelevant.  Phys¬ 
ically,  J  can  be  carried  by  normally  charged  particles  like  electrons.  The  magnetic  current 
M  can  only  be  approximated  by  normally  charged  particles,  since  magnetic  charges  evi¬ 
dently  do  not  exist. 

The  purpose  for  this  discussion  is  that  if  a  real,  steady-state,  force-free  plasmoid  were 
to  exist,  with  non-rero  fields  inside  and  free  fields  outside,  some  mechanism  would  have 
to  be  provided  to  sustain  multiple  layers  of  sheet  currents  on  the  surface  of  the  plasmoid 
if  there  are  discontinuities  in  the  fields.  No  such  mechanism  immediately  cornea  to  mind. 

On  the  other  hand,  assume  a  plasmoid  of  this  type  is  suddenly  created  with  no  surface 
currents.  The  jump  conditions  at  the  surface  can  initiate  a  mode  of  disintegration  that 
involves  transient  (one  suspects  non-harmonic)  phenomena  that  will  generate  waves  whose 
initial  energy  densities  can  be  equal  to  the  electromagnetic  energy  densities  at  the  surface  of 
the  plasmoid.  These  waves  would  linearly  add  to  the  fields  inside  and  outride  the  original 
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plasmoid.  Since  they  are  free-space  waves,  these  transients  would  not,  by  themselves, 
disrupt  the  internal  plasmoid  fields. 

Very  quickly,  then,  the  space  around  the  original,  now  decaying,  plasmoid  would 
fill  with  ultimately  outgoing  energy.  The  exact  form  of  the  outgoing  radiation  then,  in 
principle,  can  be  calculated.  It  may  then,  in  principle,  be  possible,  knowing  the  outgoing 
field  quantities  on  a  fixed  surface,  to  establish  an  externally  generated  electromagnetic  field 
that  can  set  up  a  standing  wave  that  exactly  compensates  for  the  losses  at  the  surface  of 
the  decaying  plasmoid.  A  model  like  this  is  very  close  to  that  described  by  Kapitiza'4). 

Surprisingly,  the  spherically  symmetric 
form  of  all  the  assumed  fields  allows  straight¬ 
forward  analytic  solutions  that  permit  plas¬ 
moid  properties  in  a  spherical  region  0  <  r  < 

Ro,  standing  waves  in  the  “outer”  region  Rq  < 
r  <  oo,  and  do  not  require  any  currents  on 
the  surface  at  r  =  Ro-  This  situation  is  graph¬ 
ically  shown  in  fig.(12). 

In  fig.(12)  the  region  0  <  r  <  Ro  con¬ 
tains  a  force-free  harmonic  plasmoid,  at  fre¬ 
quency  to,  with  fields  and  currents  obeying 
eqs.(9)-(14).  The  spatial  part  of  the  fields  in 
this  region  are  assumed  to  be  proportional  to 
M  in  eq.(8).  Since  r  =  Ro  is  a  zero  of  the 
Bessel  function  ji(ar),  the  fields  on  the  sur¬ 
face  of  the  plasmoid  are  directed  along  the 
polar  (6)  direction. 

In  addition,  r  =  Ro  is  a  zero  of  jt(kr).  As  a  consequence,  the  magnetic  field  of  a  TM 
cavity  mode  at  frequency  w  can  have  its  magnetic  field  equal  to  zero  at  r  =  Ro,  and  its 
electric  field  will  only  have  a  6  component.  The  dual  situation  of  a  TE  cavity  mode  can 
have  its  electric  field  equal  to  zero,  and  its  magnetic  field  will  only  have  a  6  component. 
It  is  a  simple  matter  to  adjust  the  amplitudes  of  a  TM  and  a  TE  mode  so  that  all  the 
fields  are  continuous  at  r  =  .Ro.  There  is,  then,  no  need  for  currents  at  r  =  Ro,  but  a 
need  to  maintain  an  external  electromagnetic  field  with  the  correct  boundary  conditions. 
In  principle,  this  can  be  done  without  violating  any  physical  laws. 

Alternately,  a  plasmoid  could  be  created  with  a  single  layer  of  electric  current  on  its 
surface,  which  would  allow  a  surface  discontinuity  in  the  magnetic  field.  To  obviate  extra 
layers  of  current  the  external  standing-wave  field  can  be  a  pure  TM  cavity  mode  with  zero 
magnetic  field  at  the  surface  of  the  plasmoid.  The  boundary  conditions  for  such  a  situation 
are,  in  principle,  achievable  inside  a  conducting  spherical  cavity,  as  shown  in  fig.(13).  The 
radius  of  the  plasmoid,  Ro  satifies  eq.(35),  while  the  radius  of  the  cavity,  a,  is  chosen  to 


SPHERICAL 

PLASMOID 


Figure  12.  Matching  plaamoid  boundary  condi¬ 
tions.  At  the  surface  of  a  plasmoid  the  E  and  B 
fields  are  parallel  and  tangential  to  the  surface. 
This  can  be  exactly  matched  ly  external  vacuum 
fields. 
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satisfy 


as  can  be  seen  by  examining  eq.(31). 

The  latter  problem  is  to  create  and  main¬ 
tain  a  single  sheet  of  surface  current  to  neu¬ 
tralize  jumps  in  the  surface  magnetic  field 
proportional  to  (a/k)  <  1,  see  eqs.  (11)- 
(13).  This  seems  more  tractable  compared  to 
generating  and  maintaining  a  double  layer  of 
extremely  large  currents  to  match  a  discon¬ 
tinuity  in  the  electric  field.  The  magnitudes 
of  these  large  currents  would  have  to  be  in¬ 
versely  proportional  to  the  “skin  depth”  or 
transition  region  between  the  inside  and  the 
outside  of  the  plasmoid,  and  would  have  to 
be  out  of  phase 'with  each  other,  to  simulate  a 


Figure  13.  Matching  plasmoid  boundary  condi¬ 
tion!  in  a  resonant  cavity.  By  matching  only  the 
electric  field  at  the  surface  of  a  plasmoid,  a  cur¬ 
rent  iheet,  3,  is  needed  to  account  for  the  nominal 
jump  in  the  magnetic  field. 

discontinuity  in  only  the  electric  field.  As  a 


result,  the  two  sheets  of  current  would  strongly  repel  each  other,  tending  toward  instability. 


The  general  question  of  terminating  the  surface  of  a  plasmoid  probably  cannot  be 
answered  easily  without  more  work.  Anecdotal  evidence  seems  to  indicate  that  the  au¬ 
tonomous  plasmoids  responsible  for  ball  lightning  do  not  have  external  supportive  fields 
that  need  to  match  stringent  boundary  conditions*24*,  but  measurements  are  virtually 
non-existent.  On  the  other  hand,  workers  in  Japan  have  produced  plasmoids  with  excep¬ 
tionally  long  lifetimes,  in  cavities,  sustained  by  microwave  generators.  Whether  either  of 
these  examples  can  be  described  by  force-free  time-harmonic  waves  remains  to  be  seen. 


In  Appendix  C  speculative  terms  are  added  to  Maxwell’s  equations  that  would  allow 
for  a  jump  in  either  the  magnetic  or  the  electric  field  at  a  surface  without  the  need  for  a 
current  on  that  surface. 


Apart  from  surface  boundary  considerations,  a  picture  emerges  concerning  the  descrip¬ 
tion  of  a  spherical  plasmoid  with  a  “trapped”  electromagnetic  field  inside  it.  Consistent 
with  a  current  given  by  eq.(14)  and  an  electric  field  given  by  eq.(9)  the  plasmoid  will  have 
the  complex  (capacitive)  conductivity 


Ionization  levels  throughout  the  plasmoid  must  be  such  that  the  free-electron  density  is 
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given  by  eq.(49). 


From  a  practical  point  of  view,  the  piasmoid  can  be  thought  to  be  a  linear  system 
responding  to  Ohm’s  law  in  the  form  J  =  cr(w)E,  where  J  and  E  are  complex.  (In  eqs.(9) 
and  (10)  the  sine  and  cosine  functions  could  have  been  replaced  by  the  complex  exponential 
e~'wi,  making  E  and  H  complex  functions.  Taking  the  real  parts  of  these  functions  would 
then  give  the  same  physically  meaningful  results  derived  using  real  trigonometric  time 
functions^*.) 

In  the  high-frequency  limit,  k  — *  oo,  a  piasmoid  has  the  asymptotic  specific  capac¬ 
itance  equal  to  the  vacuum  permittivity  <o  =  8.89 pf/m.  At  high  frequencies,  this  is 
equivalent  to  a  low  impedance,  meaning  that  high  field  strengths  might  be  required  to 
initiate  small  energetic  plasmoids. 

Alternatively,  the  piasmoid  could  be  thought  of  as  a  lossless  dielectric  material.  From 
eqs.(12)  and  (14),  this  material  would  have  an  effective  permittivity, 

a 2 

Uf/  =  F*’ 

Note  that  it  is  possible  that  Uff  <  «o.  implying  a  superluminal  phase  velocity  within  the 
piasmoid. 

The  piasmoid  region  will  have  internal  vortical  motion  that  will  cancel  the  tendency 
t  j  expand  due  to  electromagnetic  pressure. 

A  remarkable  “accident’’  may  therefore  exist  in  Nature  wherein  a  non-propagating 
electromagnetic  field  and  a  special  vortical  fluid  field  with  the  same  functional  form,  neither 
of  which  could  exist  separately,  can  combine  to  producing  a  long-lived  high-energy-density 
piasmoid.  The  outer  supportive  field,  since  it  can  exist  in  a  vacuum,  would  have  no  need 
for  vortical  material  motion  for  stability. 

No  radiation  can  escape  from  the  surface  of  the  piasmoid,  to  first  order,  at  the  fre¬ 
quency  to.  Radiation  may  occur  at  harmonics  of  u,  due  to  the  oscillatory  motion  described 
by  eq.(57).  This  radiation  may  have  a  detectable,  distinct,  but  perhaps  not  unique,  signa¬ 
ture. 

To  produce  the  plasmoids  described  in  this  note,  the  following  steps  must  be  included: 

(1)  Closed  spherical  vortices  must  be  produced; 

(2)  Desired  levels  of  ionization  must  be  maintained; 

(3)  Electromagnetic  energy,  at  frequency  w,  must  be  produced  within  the  piasmoid. 


35 


An  explanation  is  needed  for  (3),  above.  The  function  being  called  M(r)  in  this  note 
corresponds  to  a  sum  of  the  functions  and  on  pp.  416-417  of  Stratton’s  book  5  . 
On  page  419  of  the  same  source,  a  plane  wave  polarized  in  the  x  direction,  propagating  in 
the  z  direction,  is  expanded  in  the  form 


nsO  v 


(85) 


None  of  the  terms  on  the  right  side  of  eq.(85)  can  contribute  to  the  fields  in  the  plasmoids 
described  here.  That  is,  the  terms  in  eq.(85)  have  the  wrong  spatial  symmetries  to  con- 
tribute  to  the  plasmoid  fields.  The  same  comments  can  be  shown  to  be  true  if  the  field 
due  to  a  point  source  was  considered.  As  a  result,  external  localized  sources,  which  may 
be  expanded  in  terms  of  plane  waves  or  point  sources,  cannot  produce  energy  that  will  be 
trapped  in  a  plasmoid. 

If  however,  the  plasmoid  region  was  traversed  by  a  column  of  current  whose  axis  was 
parallel  to  the  akis  of  the  plasmoid,  the  field  of  the  plasmoid  can  be  pumped.  Another  way 
of  looking  at  this  is  that  a  circular  cylinder  of  current,  with  its  symmetry  axis  coincident 
with  the  plasmoid  axis,  will  emit  radiation  compatible  with  bong  trapped  by  the  plasmoid. 

Rom  what  has  been  derived  in  this  note,  a  possible  scenario  can  be  envisioned  for  the 
creation  of  a  spherical  plasmoid:  A  very  strong  discharge  occurs  in  an  ionizable  gas  in  the 
presence  of  a  high-frequency  electromagnetic  field.  The  pinch  effect  due  to  the  discharge 
initiates  a  compression,  which  turns  into  a  swirling  vortex.  At  the  right  level  of  ionization, 
the  external  time-harmonic  field  is  expelled  from  the  vortical  region  and  radiation  from 
the  column  of  the  discharge  inside  the  vortex  at  the  same  frequency  as  the  external  fields 
cannot  escape.  At  the  end  of  the  discharge,  a  highly  energetic  plasmoid  emerges  that  has 
nearly  as  much  kinetic  energy  as  internal  electromagnetic  energy.  Created  in  a  quiescent 
ambient  gas  with  a  low  kinematic  viscosity,  the  plasmoid  will  slowly  come  up  to  a  speed 
that  is  compatible  with  its  tangential  gas  velocity.  At  early  times,  the  plasmoid  could 
appear  to  move  very  slowly  and  may  respond  to  ambient  gas  flows. 

Catastrophic  demise  can  occur  if  the  radiation  does  not  leak  away  at  the  same  rate 
the  fluid  motion  dissipates  heat  via  viscous  forces.  The  fluid  and  electromagnetic  forces 
will  not  balance  and  the  plasmoid  may  come  apart.  On  the  other  hand,  variations  in 
the  external  fields  could  cause  boundary  mismatches  that  might  also  lead  to  catastrophic 
decay. 

Alteratively,  the  viscous  dissipation  may  closely  match  the  electromagnetic  losses,  and 
if  thermal  stresses  are  not  too  large,  the  plasmoid  may  expire  gradually. 

The  assumption  will  be  that  the  smaller  the  value  of  {,  in  eq.(51),  the  longer  the 
liforimo  of  the  excitation,  all  other  things  being  equal.  At  this  time  there  are  no  estimates 


for  lifetimes,  but  some  reasonable  values  for  the  parameters  might  be  Ro  =  10cm,  A  = 
1.0mm,  and  (  —  .01.  Under  these  conditions  a  10cm  radius  spherical,  excitation  may 
have  a  total  energy  density  on  the  order  of  10'  Joules/m'.  This  value  is  larger  than 
average  reported  energies  of  ball  lightning^24),  and  these  excitations  may  represent  an 
electromagnetically  excited  fluid  state  that  could  serve  as  a  model  for  ball-lightning-lilce 
phenomena. 

If  £  is  kept  constant,  the  total  energy  in  any  spherical  excitation  will  be  proportional 
to  the  factor 


JP=(^)4'iJo-  (86) 

This  implies  that  for  a  fixed  ratio  of  the  radius  to  the  wavelength,  for  the  same  value  of 
the  energy  content  of  a  spherical  excitation  is  linearly  proportional  to  its  radius.  Since 
the  volume  is  proportional  to  the  cube  of  the  radius,  the  smaller  the  sphere,  the  higher 
the  energy  density,  under  the  above  fixed  ratio  of  radius  to  wavelength.  For  the  numerical 
example  already  given,  if  the  radius  were  shrunk  to  one  millimeter,  the  total  energy  content 
would  be  about- 15,000  Joules.  The  wavelength  of  the  harmonic  wave  would  be  10pm;  and 
the  energy  equivalent  to  a  fifty-watt  lightbulb  burning  for  five  minutes. 
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APPENDIX  A 


THE  /  =  1  FAR  FIELD 


la  this  appendix,  the  far-field  radiation 
from  the  surface  of  an  ideal  force-free  plaa- 
moid  will  be  analyzed  for  the  1  =  1  case. 
It  will  be  shown  that,  for  appropriate  val¬ 
ues  of  £o,  the  radius,  the  fields  at  the  sur¬ 
face  of  the  plasmoid  will  not  radiate,  even 
if  the  plasmoid  is  not  terminated  by  surface 
currents.  An  important  conclusion  is  that  un¬ 
der  certain  conditions  surface  currents,  which 
are  needed  to  satisfy  Maxwell’s  equations  if 
field  discontinuities  occur  at  the  boundary, 
are  not  needed  to  prevent  steady-state  radi¬ 
ation  from  an  l  =  1  force-free  plasmoid. 

Let  the  wavenumber  be  k  =  ui/c,  with  c 
being  the  speed  of  light,  and  u  the  angular 
frequency  of  the  fields.  Referring  to  fig.(Al), 
the  Stratton-Chu  formula,  on  page  467  of 
Stratton’s  book(s),  for  the  electric  field  due 
to  the  radiation  from  the  surface  of  a  closed 
region  can  be  written 


Figure  Al.  The  plasmoid  far  field.  Radiation 
from  a  plasmoid  involves  integrating  the  fields  over 
the  surface  with  a  Green’s  function  linking  the  sur¬ 
face  with  the  far-field  point  R. 


E(r)  =  iufio  £  (nxH(r')]s(|r-  r'|)da'  +  Vx^  [nxE(r')]j(|r- r'|)da' 
+  jf  [n-E(r')]s(|r-r'|)da', 


(Al) 


where  n  is  the  outwardly  pointing  normal  to  the  region,  r  is  the  vector  distance  from  a 
fixed  origin  to  the  oberservation  point,  F  is  the  vector  distance  from  a  fixed  origin  to  the 
integration  point  on  the  surface,  the  harmonic  time  dependence  of  the  field  quantities  is 
e_lwl,  and  the  Green’s  function  is 


(A2) 


At  the  surface  of  a  plasmoid,  the  fields  have  no  components  normal  to  the  boundary, 
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so  that  the  last  integral  in  eq.(Al)  is  zero.  Since  E  and  H  rue  identical,  within  a  constant 
factor,  only  the  first  integral  in  eq.(Al)  has  to  be  evaluated.  It  will  be  shown  that  the 
value  of  this  integral  is  also  zero  if  the  radius  is  chosen  correctly.  Showing  that  this  integral 
is  zero  in  the  feu  field  is  equivalent  to  showing  that  the  term  involving  the  second  integral 
is  also  zero.  The  curl  of  a  null  function  is  zero. 

The  strategy  will  be  to  first  show,  rigorously,  that  the  value  of  the  first  integral,  call 
it  Vi,  can  be  zeroed  in  the  far  field,  «.£.,  as  r  — >  oo. 

In  fig.(Al),  as  r  -*  oo,  the  quantity 


|r  —  r'l  ~  ft  —  r'-r, 


(A3) 


where  R  =  |r|,  r  is  the  unit  vector  in  the  direction  of  r,  and  an  angle,  7,  is  defined  by 


r'-r  =  |r'|cos7. 


(A4) 


Now  define  an  equivalent  surface  current,  J,  where 


J  =  nxH. 


(AS) 


The  integral  to  be  evaluated  becomes 

,tkR  t 

V,  =  Js J(r')e",i|r  'C°nda'-  (-6) 

In  eq.(A6)  the  integral  extends  over  the  entire  surface  current  distribution. 

In  the  l  =  1  spherical  plasmoid  case,  this  equivalent  surface  current  at  radius  Ra  can 
be  written  as  a  vector  function  proportional  to 


J.«r/  =;o(aUo)sin0^, 


(A7) 


where  jo(aRo)  =  —.2172 . . .  corresponds  to  the  value  at  the  first  zero  of  ji(ar).  The  first 
zero  was  chosen  arbitrarily,  any  other  value  of  r  for  which  j;  ( ar )  =  0  would  give  the  same 
result. 
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The  value  of  Vi  becomes  proportional  to  the  integral 

I  =  ^sjne'e-,*R»co,W,  (48) 

where  the  distribution  is  centered  at  r1  =  0,  standard  spherical  coordinate  notation  is  used, 
and  the  integration  is  over  the  surface  of  the  sphere  jr*  j  =  Rq. 

The  integral  in  eq.(A8)  can  be  evaluated  in  closed  form  by  expanding  all  quantities 
in  the  integrand  in  terms  of  spherical  harmonics*’9)  Ti,m(n'),  in  the  primed  coordinates, 
and  in  the  unprimed  (fax-field)  coordinates.  First  convert  the  spherical  vector  ft 

into  a  sum  of  unit  Cartesian  vectors  x  and  y  in  the  x  and  y  directions: 


Next  note  that 


ft  =  -ianft  +  ycoeft. 


(49) 


sins'  sin*'  .  ^/f§|Vw(<0  +  ri.-i(O')] . 


(410) 


and 


sin  <f  cos  ft 


r»,,(tf)+n.-i(*0]. 


(411) 


The  exponential  in  the  integral  has  the  expansion*’9) 

«-'*** co,T  =  4x f; i-'j,(kRo)  £  y^(n')y,,m(fi),  (412) 

ImO  mm— l 


where  the  asterisk  (*)  denotes  complex  conjugation.  Using  the  orthogonality  property  of 
the  spherical  harmonics, 


jf  yl;m(fi')yl.,m.(fi')dn' = 


(413) 
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a  little  algebra  shows  that,  finally, 


1=  -AfKj^kRo)^, 


(414) 


where  is  the  unit  spherical  vector  in  the  far-field. 

If  the  product  kR>  is  one  of  the  zeroes  of  ji(kHo)  then  the  amplitude  of  the  outgoing 
electric  field  due  to  the  l  =  1  surface  radiation  goes  to  zero  more  quickly  than  the  inverse 
of  the  distance  in  the  far-field.  The  fields  at  the  surface  of  a  plasmoid  with  this  radius 
will  therefore  not  radiate.  An  exactly  analogous  argument  shows  that  the  l  =  1  magnetic 
surface  field  will  not  radiate. 

There  is  an  interesting  interpretation  to  this  result.  It  was  shown  in  Sec.  IX  that  it  is 
posable  to  have  physically  realizable  fields  outside  a  plasmoid  such  that  no  discontinuities 
occur  at  the  boundary.  In  general,  if  such  a  situation  would  arise,  the  external  field  would 
be  exchanging  energy  with  the  plasmoid,  even  though  the  net  exchange  would  be  zero.  At 
the  radii  that  correspond  to  zeroes  of  the  spherical  Bessel  function,  call  them  resonances, 
this  exchange  of  energy  is  halted.  It  would  appear  that  under  such  conditions  plasmoids 
would  be  more  mechanically  stable.  The  reasoning  for  this  is  that  momentum  transfer  and 
energy  transfer  are  related. 

Away  from  resonance,  the  energy  and  momentum  balances  are  maintained  by  ex¬ 
changes  between  the  inride  and  outride  of  a  plasmoid.  Disruption  of  either  the  flow  from 
the  inside  to  the  outride,  or  the  outride  to  the  inride,  will  produce  a  net  momentum  trans¬ 
fer  rate,  which  is  equivalent  to  a  net  force.  At  resonance,  forces  arising  from  fluctuations 
of  the  fields  would  thereby  be  minimized. 
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APPENDIX  B 
A  NUMERICAL  EXAMPLE 


On  the  basis  of  anecdotal  evidence,  very  energetic  plasmoids,  known  as  “Ball  Light¬ 
ning”  may  have  total  energy  contents  of  one  megajoule  with  a  radius  of  10cm.  With 
no  further  justification,  it  will  be  instructive  to  calculate  a  set  of  parameters  that,  could 
produce  such  a  force-free  plasmoid  in  the  Earth’s  atmosphere  at  sea  level  and  at  0°C. 

For  a  sphere  with  R o  =  10cm  ,  from  eq.(8),  assuming  an  /  =  1  plasmoid,  with  a  radius 
at  the  first  zero  of  ji(ar), 

a  =  44.93.  (B1) 

Assuming  k  will  be  large  enough  to  ignore  a3/k3,  then  the  total  energy  will  be  shared 
nearly  equally  by  the  electromagnetic  field  and  the  fluid  kinetic  energy,  each  having  about 
half  a  megajoule.  Expressed  as  an  integral 

^  |M(r)|2ar  =  .5x10®,  (B2) 

4  Jo 

where  M(r)  is  defined  by  eq.(8),  with  l  =  1.  Direct  integration  shows  that 

f  HMf*-*?.  (**> 

so  that 

Eo  =  3.566x  lO’volts/meter.  (-B4) 

Taking  £  =  .01  in  eq.(58), 

w = y/1-759xlQQ11xx31566-— = 7.920x10”,  m 

is  equivalent  to  a  126GF*  field  with  it  =  2642meters“1  To  satisfy  the  second  part  of  eq.(36), 
this  value  of  k  must  be  changed  slightly  to 


T 


k  =  2686(meters)  *, 


m 


which  is  equivalent  to  a  12 8  GHz  field. 

-  -  FVom  eq.(38)  the  charge  density  is 


pe  —  — 32.63coul/m3,  (B  7) 

amounting  to  8.532xl017  electrons.  Considering  that  there  are  about  1.126xl023  diatomic 
molecules  of  gas  in  the  plasmoid,  on  the  order  of  only  one  to  two  molecules  in  105  needs 
to  contribute  an  electron. 

At  the  surface  of  the  sphere,  from  eqs.(5),  (20)  and  (27)  the  speed  of  the  gas  u,  is  a 
function  of  the  polar  angle  6: 


u,  =  1.431  x  103  sin  0(meters/sec). 


(B8) 


While  u,  represents  a  very  large  speed,  it  is  not  true  that  the  plasmoid  will  be  traveling  at 
955  meters/sec,  the  equivalent  solid-sphere  speed  to  match  the  surface  speed.  The  reason 
is  that  the  plasmoid  was  assumed  to  be  created  in  a  quiescent  atmosphere.  Analysis  of 
the  Navier-Stokes  equation*7^  shows  that  the  boundary  layer  emerging  from  the  surface  of 
the  sphere  as  a  result  of  the  surface  motion  diffuses  like  a  thermal  wave.  To  grow  out  to  a 
distance  of  one  extra  Bo  takes  on  the  order  of 


At  =  R%/v  =  .01/1.32X  IQ’5  =  757sec,  (B9) 


where  v  is  the  kinematic  viscosity  of  air  at  1  Atm.  and  0°C. 

In  the  first  few  seconds,  the  plasmoid  will  therefore  hardly  be  moving  relative  to  the 
quiescent  atmosphere.  On  the  order  of  tens  of  seconds  some  motion  will  be  evident,  but 
it  is  not  clear  that  the  plasmoid  has  a  lifetime  long  enough  to  get  up  to  any  appreciable 
speed.  Its  entire  liftime,  then,  may  be  spent  “floating”,  or  bring  carried  by  sir  currents, 
or  responding  to  electrodynamic  forces. 

It  is  instructive  to  scale  the  radius  down  to  50/im.  The  total  energy,  assuming  the 
a/k  ratio  stays  the  same,  might  naively  be  projected  to  500  joules.  Since  the  charge 
density  will  go  up  inversely  as  the  square  of  the  radius,  the  density  of  the  plasmoid  must 
be  increased.  Since  only  about  one  molecule  in  10s  will  be  ionized  in  the  larger  system, 
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a  density  corresponding  to  about  forty  times  sea-level  atmospheric  air  must  be  assumed, 
which  will  be  totally  ionized.  The  wavelength  of  the  trapped  radiation  will  be  A  =  1.17pm, 
in  the  infrared  range.  The  surface  speed,  eq.(B8),  *— ill  be  u,  —  7.16  x  104  sin  6,  with  a 
maximum  internal  speed  of  2.2xlOsmeters/sec,  still  not  relativistic.  The  number  of  free 
electrons  will  be  4.3X1013 .  It  would  take  on  the  order  of  8ms  for  such  a  plasmoid  to  come  up 
to  speed  at  0°C.  At  different  temperatures,  where  the  kinematic  viscosity  may  be  higher, 
this  time  could  be  much  shorter.  The  assumption  has  been  that  the  ambient  atmosphere 
is  also  at  a  density  forty  times  that  of  sea-level  air.  It  is  not  clear  how  the  stability 
will  be  affected  if  the  plasmoid  travels  into  a  lower-density  and/or  pressure  region.  The 
electromagnetic  pressure,  balanced  by  the  vortical  field,  is  almost  1010  atmospheres.  If 
even  one  joule  is  trapped  the  pressure  is  on  the  order  of  107  atmospheres.  It  would  seem 
that  a  perturbation  of  a  few  atmospheres  of  pressure  may  not  be  immediately  destabilizing. 
More  work  must  be  done  in  this  area. 
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APPENDIX  C 


AD  HOC  FIELD  TERMS 


In  this  appendix,  terms  will  be  added  to  Maxwell’s  equations  that  will  allow  true 
discontinuities  to  arise  in  the  tangential  components  of  electromagnetic  fields  at  surfaces, 
generating  accompanying  surface  conditions  that  cannot  be  derived  from  the  ‘‘usual”  set  of 
Maxwell’s  equations.  It  should  be  understood  that  there  may  be  no  physical  justification 
for  these  terms.  They  serve  as  an  illustration  that  it  is  possible  to  assume  new  operators 
in  Maxwell’s  equations  that  will  not  affect  any  solutions  that  are  continuous,  yet  can  give 
a  contribution  at  a  field  discontinuity. 

The  terms  to  be  eventually  introduced  are  classically  unobservable,  because  the  elec¬ 
tromagnetic  field  is  classically  assumed  to  be  analytic  everywhere.  An  objection  may  be 
raised  that  the  introduced  terms  are  not  Lorentz  invariant.  In  fact,  they  are.  Lorentz 
invariance  will  be  investigated  separately.  The  important  point  is  that  classically  unob¬ 
servable  terms  can  be  added  to  Maxwell’s  equations,  causing  observable  deviations  from 
the  usual  formulation  at  true  field  discontinuities. 

Consider  the  operator,  O,  acting  on  the  vector  function  M(r),  at  a  point  r', 


OM(r)  =  J  <fflrxr-VM(r) 


(Cl) 


where  r  is  the  unit  radial  vector  centered  at  r',  and  the  integration  is  over  the  full  4 jr 
stear adians.  From  the  properties  of  the  gradient  operator,  this  can  be  written, 


OM(r)  =  J  dflr  x  ~M(r) 


(C2) 


Integration  in  eq.(C2)  is  over  the  angular  parts  of  the  integrand.  The  operator  d/dr 
does  not  affect  the  angular  parts  of  M(r).  In  general,  any  wavefunction  obeying  Maxwell’s 
equations  can  be  expanded  as  a  sum  of  multipole  fields  in  the  form^5) 


M(r)=  +  a£iz!+i(r)Tj,i+i,M  +  a'j'iz,ii(r)Tj1,_1,M  J 


(C  3) 


where  Tjt is  a  vector  spherical  harmonic  (J5*,  a'j^'a  are  complex  constants,  and 


s/(r)’s  are  (unctions  of  the  radial  distance  from  the  center  of  the  expansion.  The  descrip¬ 
tion  of  a  wavefunction  using  eq.(C3)  is  also  known  as  a  Helmholtz  decomposition. 

Vector  spherical  harmonics  are  orthogonal  on  the  surface  of  the  unit  sphere,  ».e., 

J  =  Sj'j'St.t'&M.M',  (C4) 

where  the  asterisk  denotes  complex  conjugation. 

The  unit  radial  vector  can  be  expressed  as 

r=-±=T0,i,0(n).  (C5a) 
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Similarly,  the  unit  azimuthal  vector  can  be  written 

(C54) 

Prom  the  rules  for  coupling  representations  of  the  rotation  group*23),  the  only  terms 
in  the  multipole  expansion  of  M(r)  that  can  contribute  to  the  integral  in  eq.(C2)  are 
those  associated  with  Ti,i,mi  — 1  <  m  <  1.  When  expanded  about  different  origins,  the 
vector  function  M(r)  will  appear  to  have  differing  amplitudes  associated  with  the  the 
vector  spherical  harmonics  The  amount  of  at  a  point  of  expansion  of  M(r) 

is  proportional  to  V  X  M(r)  at  that  point.  The  integral  in  eq.(C2)  will  generally  define 
a  different  vector  at  each  point  r1,  depending  on  the  local  expansion  of  M(r)  about  the 
point  r*.  If  this  local  expansion  does  not  have  harmonic  terms  of  the  form  there 

will  be  no  contribution  at  that  point. 

It  is  now  clear  what  form  an  operator  can  take  to  give  no  contribution  when  applied 
to  analytic  vector  fields  arising  horn  solutions  of  the  classical  Maxwell  equations: 

m(r)=  idnix~jM(T)-  £  <Tu,m,M(r))TuJ  ,  (C6) 

*  m=!-l  irmr1 

where 


46 


(C7) 


In  eq.(C6),  O  is  a  linear  operator  that  gives  a  null  result  when  applied  to  an  analytic 
multipole  field.  At  a  discontinuity  in  the  field  M(r)  the  radial  derivative  term,  in  eq.(C2), 
will  produce  a  singularity  in  the  neighborhood  of  the  discontinuity.  It  should  also  be  noted 
that  removing  harmonics  of  the  form  Ti(iifn  does  not  affect  the  vector  field  M(r;  at  |r|  =  0 
because 


lim  zf  (r)  =  0;  1^0. 


(C  8) 


Therefore 


lim 
|r|  *0 


M(r)-  £  <T 


=  M(r) 


(C  9) 


which  means  that  if  the  function  M(r)  has  a  discontinuity  at  r  =  r'  the  integral  in  eq.(C2) 
will  define  a  generating  function  for  jumps  in  the  field  M(r),  in  the  sense  of  the  theory  of 
distributions(J4).  Less  rigorously,  the  integral  in  eq.(C2)  behaves  like  a  delta  function  at 
discontinuities  in  M(r),  is  finite  at  discontinuities  of  the  gradient  of  M(r),  and  is  otherwise 
zero. 


Maxwell’s  equations  can  be  rewritten  for  any  given  time-harmonic  frequency  ui,  in 
vacuo,  incorporating  the  operator  O, 


VxH(r)  +  GhOH(t)  =  -iwe0E(r)  +  J(r),  (CIO; 


xnd 


^ +  GeOE(t)  =  twfj|jH(r),  (Cll) 


where  the  time-harmonic  factor  is  exp(  -iwt),  and  the  numerical  factors  Gh  and  Ge  will 
be  determined  below. 

The  constants  Gh  and  Gb  can  be  determined  by  assuming  that  no  current  will  flow 
on  the  surface  defining  a  discontinuity  in  either  E  or  H.  The  usual  method  for  deriving 
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the  current  at  a  boundary  can  be  modified  to  find  Gu  and  6’g.  Enclose  a  small  portion  of 
the  boundary  between  two  regions  where  either  H  or  E  has  a  jump  discontinuity  with  a 
closed  curve  in  a  plane  that  cuts  through  the  boundary.  In  the  usual  proof,  both  sides  of 
eq.(ClO),  for  example,  are  integrated  over  the  area  of  the  closed  curve,  Stoke’s  theorem  is 
applied,  and  the  current  flowing  through  the  area  is  found  to  be  equal  to  the  line  integral  of 
the  component  of  II  tangent  to  the  closed  curved.  At  a  discontinuity,  the  current  flowing 
through  the  closed  curve  is  proportional  to  the  inference  in  tangential  components  of  H 
on  both  sides  of  the  boundary  of  the  discontinuity,  i.e., 


Jj»r/  ~  nx(Hj  —  Hj), 


(C12) 


where  n  is  the  unit  normal  to  the  boundary,  pointing  from  the  region  with  the  field  H* 
into  the  region  with  the  field  Hj. 

At  a  field  discontinuity,  the  integral  of  OH(r)  over  the  surface  of  the  closed  curve, 
discussed  above,  is  proportional  to  the  right-hand  side  of  eq.(C12).  It  can  be  shown  that 
if 


=  <?£  =  —,  (£713) 

no  current  will  flow  tangential  to  true  field  discontinuities. 

It  can  be  argued  that  Maxwell’s  equations,  including  O  terms,  with  eq.(C13),  are 
Lc-entz  invariant.  The  logic  is  as  follows:  All  observers  moving  relative  to  an  analytic 
c.,ctromagnetic  field  will  conclude  that  the  field  they  see  is  also  analytic,  and  the  O  terms 
will  not  contribute.  On  the  other  hand,  any  true  discontinuity  in  either  the  electric  or 
magnetic  field  will  be  described  differently  by  different  observers,  but  all  observers  will 
agree  that  no  current  would  be  flowing  tangentially  to  the  boundary  between  the  regions 
where  the  discontinuity  occurs. 

The  reader  should  again  be  cautioned  that  the  added  O  terms  in  Maxwell’s  equa¬ 
tions  may  have  no  physical  justification.  Experimentally,  effects  such  as  ball  lightning  and 
anomalous  charge-cluster  formation  occur  in  the  vicinity  of  very  strong  local  discharges 
that  may  cause  field  discontinuities.  The  purpose  of  this  appendix  is  to  show  that  such  ef¬ 
fects  can  arise  from  terms  that  are  classically  unobservable,  causing  no  changes  whatsoever 
in  any  other  solutions  where  true  discontinuities  do  not  arise. 

Applied  to  the  force-free  case,  for  cases  that  functionally  look  like  eq.(8)  of  the  main 
text,  at  a  radius  corresponding  to  a  zero  of  the  function  ji(ar)  the  fields  are  tangential 
to  the  spherical  surface.  A  discontinuity  at  this  surface  equivalent  to  a  drop  to  zero  fields 
outside  the  plasmoid  with  large  fields  inside,  will  not  need  a  set  of  boundary  currents  to 


48 


sustain  the  wavefunctions.  Indeed,  the  remarkable  special  case,  called  a  vactiod  in  the  text, 
wherein  a  spherical  nonpropagating  field  is  “trapped”  in  the  vacuum,  becomes  realizable. 

If  the  above  theory  applied  at  zero  frequency  to  “perfect”  superconductors,  where 
the  transition  from  non-zero  field  to  zero  field,  outside  to  inside,  occurred  as  a  true  dis¬ 
continuity,  there  would  be  no  current  residing  on  the  surface.  Since  the  quantity  nxfi 
appears  in  the  equations  for  the  external  field,  not  the  current,  per  se^,  there  would  be 
no  obvious  contradiction  introduced  into  classical  electromagnetism.  The  logical  extension 
of  such  assumptions  would  be  the  existence  of  zero-frequency  vactoids,  having  a  static 
electric  or  magnetic  field  trapped  in,  say,  a  toroidal  closed  loop  with  no  enclosing  currents. 
Any  further  speculation  concerning  such  pathological  states  of  the  electromagnetic  field  is 
beyond  the  scope  of  this  report. 

At  the  present  state  of  development  of  technology,  it  is  impossible  to  either  accept 
such  extra  terms  or  reject  them  out  of  hand. 
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APPENDIX  D 


TOTAL  PLASMOID  ENERGY 


Equation  (23)  expressed  the  energy  in  the  electromagnetic  field  of  a  plasmoid  as  having 
a  time  dependence.  Since  energy  must  be  classically  be  conserved  instantaneously,  the  time 
variations  must  be  accounted  for  by  the  kinetic  energy  of  the  electrons  oscillating  along 
electric  field  lines.  When  this  energy  is  added  to  Q  in  eq.(23),  the  time  dependence  of  the 
total  energy  of  the  plasmoid  at  any  point  should  be  independent  of  time. 

This  is  easily  proven.  The  kinetic  energy  of  an  individual  electron  at  the  point  r  is 

K.  =  |m,|v|2,  (131) 

where  mt  is  the  mass  of  the  electron  and  v  is  its  (non-relativistic)  velocity.  The  kinetic 
energy  of  the  electron  gas  can  be  written  in  terms  of  the  current  density  at  r, 

H(t)  m 


where  e  =  -1.602x10  19C. 

Using  eq.(37),  and  the  time  dependence  of  the  current, 


T  =  ~  cos  2(wt  +  <p)  J  J(r)-E(r )dv.  (DZ) 

Since  the  spatial  part  of  the  vector  potential  can  be  related  to  spatial  part  of  the  electric 
field, 


A(r)  =  ~E(r), 

Uf 


m 


then 


cosJ(u)t  +  ifi)  J  A(r)-J(r)du. 


m 
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The  quantity  Q  +  T  >0,  which  is  the  total  energy  of  the  plasmoid,  is  now  independent  of 
time. 

Although  the  total  energy  is  positive,  two  parts  of  the  sum,  i.e.,  the  second  integral 
term  in  eq.(23)  and  the  electron  kinetic  energy,  are  equal  and  opposite  in  sign.  Even 
the  local  energy  densities  associated  with  these  terms  are  equal  and  opposite  in  sign,  the 
kinetic  energy  density  always  being  positive.  Only  frequencies  above  the  critical  frequency 
will  be  considered,  since  electrons  then  carry  the  current. 

The  free  electrons  making  up  the  current  can  be  considered  as  moving  about  the 
minimum  of  a  potential  well  whose  average  depth  at  any  point,  r,  is  (l/4)A(Vr)- J(V>). 
From  eq.(25),  the  depth  of  this  equivalent  potential  well  increases  quadratically  as  the 
frequency  increases.  The  argument  made  in  the  main  text  is  that  if  the  electronic  kinetic 
energy  is  not  changed,  a  small  local  geometric  perturbation  will  raise  the  energy  of  the 
plasmoid.  That  argument  was  based  on  the  out-of-phase  relationship  of  J(r)  and  A(r). 

Another  simple  argument  can  be  made  that  geometrically  perturbing  the  plasmoid 
raises  its  energy.  Following  GekW31)  (p.80)  each  oscillating  electron  can  be  thought  of  as 
an  oscillating  electric  dipole,  whose  potential  energy  is  at  a  minimum  when  its  dipole  mo¬ 
ment  is  aligned  with  the  fields.  If  the  perturbation  does  not  leave  the  plasmoid  in  another 
force-free  state,  a  torque  is  needed  to  turn  any  given  dipole  to  a  different  orientation.  This 
torque  will  have  components  arising  from  the  electric  as  well  as  from  the  magnetic  field. 
A  restoring  torque  will  tend  to  line  up  the  dipole  moment  with  the  electric  field  and  force 
the  electronic  motion  to  follow  magnetic  field  lines.  Since  both  the  electric  and  magnetic 
fields  are  parallel,  all  the  contributions  to  the  restoring  torque  in  a  small  perturbation  will 
add  constructively  and  be  opposed  to  the  perturbing  motion. 
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ABSTRACT 

Elementary  theory  of  force-free  fields  is  applied  to  electromagneti- 
cally  excited  charge  densities.  It  is  shown  that  mechanically  stable 
non-radiating  modes  may  be  excited  in  plasmas  by  harmonic  fields, 
even  in  the  absence  of  D.C.  magnetic  fields.  The  free-charge  den- 
-  -  sity  needed  to  sustain  a  stable  spherical  non-radiating  mode  can  be 
shown  to  be  uniform,  and  dependent  only  on  the  electromagnetic 
wavenumber  and  spatial  extent  of  the  mode,  in  a  neutral  plasma. 
The  time-averaged  electric  and  magnetic  fields  in  such  modes  are  all 
zero.  A  virial  theorem  check  verifies  the  consistency  of  an  intuitive 
approach  and  a  mote  rigorous  lagrangian  formulation.  Electromag¬ 
netic  stress  effects,  which  can  be  accurately  balanced  by  vortical  fluid 
motion,  are  ignored. 

L  INTRODUCTION 

If  the  vector  field  H(r)  is  a  force-free  field  then 


VxH  =  oH,  (1) 

where  a  is  a  complex  constant.  It  may  seem  strange  to  imagine  a  vector  field  that  is 
proportional  to  its  own  curl,  without  the  field  being  identically  null  Nonetheless,  such 
fields  can  be  constructed,  at  least  mathematically,  having  non-zero  magnitudes.  The 
properties  of  such  fields  may  be  very  valuable  in  diverse  areas  of  engineering.  For  example, 
if  the  windings  of  a  D.C.  electromagnet  followed  the  streamlines  of  a  feme-free  field  the 
force  on  the  wires,  due  to  the  magnetic  field  they  produced,  would  be  zero.  (Hence  the 
name  for  the  fields.)  In  fluid  dynamics  the  contribution  of  the  v  x  V  x  v  term  in  the 
Navier-Stolces  equation  is  zero  in  a  vortical  fluid  mode  whose  flow  follows  streamlines  of  a 
force-free  field.  (Such  modes  are  known  as  Beltrami  vortices  in  fluid  dynamics/1)) 
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In  electromagnetic  engineering,  electromagnetic  excitations  that  obey  eq.(l)  can  be 
shown  to  exist  in  vacuo  as  well  as  in  the  presence  of  charge  and  currents.  These  modes  are 
stationary,  i.e.,  provided  with  the  proper  boundary  conditions,  they  will  not  propagate. 

In  this  note,  some  of  the  simple  properties  of  force-free  non-radiating  fields  will  be 
derived.  Any  references  to  a  “non-radiating  mode”  in  this  note  will  always  imply  reference 
to  a  “force-free  non-radiating  mode”. 

Applications  will  be  made  to  the  theory  of  stable  electromagnetically  self-confining 
excitations  that  may  have  long  lifetimes  and  very  high  free-charge  and  energy  densities. 
Many  of  the  phenomena  to  be  described  may  never  have  been  observed.  Some  of  the  results 
of  the  simple  theory  may  explain  experiments  performed  by  Bostick  and  his  associates,  as 
well  as  results  reported  by  Puthoff  and  Shoulders*2*. 

The  idea  of  applying  the  force-free  mode  to  confining  plasma  is  not  new*3'4*.  Tech¬ 
niques  employing  names  such  as  the  “Spheromak”  and  “Compact  Toroid”,  which  nomi¬ 
nally  depend  on  force-free  plasma  motion,  have  been  implemented  with  varying  degrees  of 
success*5'8*.  In  these  schemes  the  object  is  to  induce  the  gross  fluid  velocity  of  the  plasma 
to  conform  to  a  Beltrami  vortical  mode.  This  lowers  the  magnetic  interaction  energy  and, 
presumably,  raises  confinement  times. 

A  common  feature  of  the  past  applications  of  force-free  modes  to  plasmas  is  that  the 
fluid,  or  ion,  velocity  is  non-zero  and  directed  along  the  ergptic  streamlines  of  a  Beltrami 
vortex.  The  stationary  (zero-frequency)  magnetic  field  also  lies  along  these  streamlines. 
Since  the  plasma  is  assumed  to  be  neutrally  charged,  Maxwell’s  equations  predict  that  the 
zero-frequency  magnetic  field  is  not  necessarily  accompanied  by  an  electric  field. 

In  what  follows,  the  force-free  idea  is  taken  one  step  farther,  with  the  introduction  of  a 
high-frequency  electromagnetic  field.  A  self-consistent  solution  will  be  shown  to  exist  that 
allows  the  electric  and  magnetic  fields,  as  well  as  the  current  to  have  the  force-free  form. 
Beyond  a  critical  frequency,  the  magnetic  interaction  energy  of  the  oscillating  free-electron 
density  is  not  only  minimized,  it  changes  sign  and  becomes  negative. 

An  initial  physical  model  is  postulated  that  has  the  ions  stationary,  ie.,  there  is  no 
net  fluid  motion,  only  the  free  electrons  oscillate  very  nearly  along  force-free  streamlines. 
Applying  Lagrange’s  equations  to  this  model  forces  the  free  electron  density  to  have  a 
specific  value  with  respect  to  the  size  parameter  a  and  the  frequency,  but  independent  of 
the  amplitudes  of  the  fields.  It  will  turn  out  that  if  a  is  a  constant,  uniform  free-charge  and 
ion  densities  are  consistent  with  force-free  high-frequency  electromagnetic/fluid  modes. 

This  model  is  incomplete,  owing  to  the  f  ict  that  the  free  electron  gas,  with  non-zero 
mass  electrons,  cannot  exactly  oscillate  along  curved  electric  field  lines,  and  will  thereby 
experience  electromagnetic  stresses.  These  stresses  are,  at  high  frequencies,  equivalent 
to  a  positive  pressure  field  proportional  to  the  square  of  the  electric  field*20*.  A  direct 
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application  of  Bernoulli’s  theorem*21!  says  that  introducing  a  fluid  velocity  whose  square 
is  proportional  to  the  square  of  the  electric  field,  i.e.,  a  Beltrami  vortex,  will  counterbalance 
the  major  electromagnetic  stresses.  The  introduction  of  fluid  vortical  motion  for  enhanced 
stability  will  not  be  addressed  in  this  note,  but  will  not  affect  the  results  derived  here. 

The  virial  theorem  predicts  that  no  excitation  can  be  confined  indefinitely  by  only 
classical  electromagnetic  fields,  even  in  the  relativistic  limit*T,a,*,1°!.  The  model  proposed 
here  does  not  violate  this  prediction.  Application  of  the  virial  theorem,  ignoring  radiar 
tion  from  surface  currents  and  non-linear  effects,  yields  conditions  for  long-term  stability 
that  are  identical  to  those  derived  for  self-consistent  force-free  motion  of  non-zero  mass 
electrons  in  an  electromagnetic  field. 

The  major  approximation  made  is  that  the  massive  electrons  exactly  follow  the  curved 
force-free  electric  field  lines.  In  a  sequel,  this  approximation  and  an  analysis  of  the  virial 
will  be  pursued^11!.  It  will  be  shown  that  there  may  be  tradeoffs  between  energy  density, 
frequency  and  confinement  time,  for  a  fixed-sized  excitation. 

H.  KINEMATICS  OF  THE  FIELDS 

Tbking  the  curl  of  eq.(l), 


VxVxH  =  VVH-VJH, 


(2) 


and  realizing  that,  since  H  is  derived  from  a  curl. 


VH  =  0, 


(3) 


then  it  follows  that  H  obeys  the  vector  Helmholz  equation, 


V2H  +  a2H  =  0. 


(4) 


Equation  (4)  is  the  key  to  searching  for  solutions  to  eq.(l).  Not  all  solutions  of  eq.(4) 
will  be  solutions  of  eq.(l),  but  all  solutions  of  eq.(l)  will  be  solutions  of  eq.(4).  The  process 
of  generating  solutions  to  eq.(l)  can  be  approached  in  several  ways.  One  way  is  to  choose 
any  vector  field  that  obeys  eq.(4)  and  successively  take  its  curl,  noting  all  new  terms  that 
arise.  If  the  number  of  independent  terms  stops  increasing,  form  a  linear  combination  of 
the  final  set  of  terms  such  that  the  linear  combination  obeys  eq.(4). 

Starting  from 
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Pi  =  sin'  8ji(cxr)4>, 


(5) 


I? 


the  resultant  force-free  vector  field  is 


Fj.bM,*)  =  (/  +  l)co6flam'-1  -  [j,_1(ar)  -  /^l]  sin' 88 

±jt(ar)tanl  8j>, 


(6) 


where  1  is  any  integer,  l  >  0,  jt  is  the  familiar  spherical  Bessel  function,  and  the  vectors 
r,  8,  and  $  are  the  familiar  right-handed  3-D  polar  coordinate  unit  vectors.  The  vortical 
field  described  by  eq.(6)  has  a  useful  application,  as  will  be  shown  below. 

It  is  somewhat  tedious,  but  straightforward,  to  show  that 


VxF  i,B(r,0,^)  =  ±oF||B(r,0». 


(7) 


/ 


A  somewhat  more  elegant  expression  for  force-free  fields  can  be  derived  from  the 
properties  of  the  rotation  group.  Using  the  notation  in  Rose’s  boolc<12)  for  the  vector 
spherical  harmonics  Tji.jf,  force-free  fields  can  be  immediately  written  that  correspond 
to  the  (21  +  l)-dimensional  representations  of  the  three-dimensional  rotation  group: 


B/,ro  =  ±ji(ar)Ti,i,m  +  ^-1— j[V1  +  lji-i(ar)Ti,i-i,m  -  v/iji+i(ar)Ti,i+i,m],  (8) 

where  B;,ro  is  a  force-free  field  for  1  >  0  and  -I  <  m  <  l. 

It  is  not  clear  at  this  time  that  the  Bi,m  form  a  complete  set  for  the  expansion  of 
general  force-free  fields.  Note  that  the  spherical  Bessel  functions  can  be  replaced  by 
spherical  Hanltel  functions,  and  that  a  does  not  have  to  be  real.  Note  also  that  force-free 
fields  can  be  expanded  in  cylindrical  (p,  j>,i)  coordinates;  and  may  be  synthesized  using 
functions  obeying  the  scalar  Helmholz  equation,  yet  possessing  angularities  appropriate 
to  conical  geometries.  Neither  of  these  latter  solution  sets  will  be  pursued  in  this  note. 

For  constant  values  of  a,  Zaghloul  and  Barajas  employ  Hansen’s  straightforward 
method  to  synthesize  force-free  vector  fields  from  arbitrary  solutions  of  the  Helmholtz 
equation.  The  reader  is  referred  to  their  paper  and  Stratton’s  book  for  details(1,’u>1,l. 
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The  real  fields  F„,a,  from  eq.(6),  can  be  shown  to  be  a  subset  of  the  set  of  fields  B/,m 
from  eq.(8).  The  explicit  forms  displayed  in  eq.(6)  for  the  field  components,  however,  are 
easier  to  deal  with  analytically. 

Consider  a  vortical  velocity  field  with  the  form  F;,b.  Within  a  common  constant 
factor  imparting  the  correct  dimensionality,  the  time  rates  of  change  of  the  coordinates  are 


dr 

dt 


=  (/  -f  1)  cos  0  sin1-1 


(9) 


f  =  sin'tf. 


(10) 


and 


=  ~jl(ar)  sin1-1 6.  (11) 

Starting  with  eqs.  (9)-(ll)  it  is  possible  to  integrate  the  trajectory  of  a  point  having 
the  above  components  of  velocity.  The  time,  f,  can  be  eliminated  from  any  pair  of  the 
expressions.  If  eq.(10)  is  divided  by  eq.(9),  the  resultant  equation 


d9 

dr 


l  +  l 


tanS 


1 

r 


jl-l(ar) 

Jl(ar) 


can  be  integrated  analytically  to  obtain 


(12) 


sinS  =  A[arji(or)]  r^r, 


(13) 


where  A  is  a  constant. 

Equation  (13)  coupled  with  the  non-analytic  equation  for  <f>, 


i  =  ± 


o  A 

1+ 


l r 

Ur. 


dr' 


r°  [ar';j(ar')] 


(14) 


determine  streamlines. 
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A  very  important  property  of  eqs.  (9}-^ll)  is  that  for  at  every  value  of  rm  for  which 
j l(orrm)  =  0,  the  radial  (r)  and  azimuthal  ($)  components  of  the  vector  fields  Fj,s  vanish. 
At  any  such  r  =  rm  the  only  component  of  the  field  is  the  longitudinal  (6)  component. 
Vortical  fields  like  F|,b  have,  then,  a  layered  structure.  It  follows  that  if  the  simplified 
longitudinal  boundary  condition  at  any  of  the  rm  can  be  met  everywhere  on  the  surface 
of  a  sphere  then  the  field  internal  to  the  sphere  will  obey  eqs.  (9)-(ll)  at  all  times  if  they 
obeyed  these  equations  and  eq.(l)  for  r  <  rm,  or  if  force-free  fields  represented  stable 
configurations  for  the  physical  system  under  consideration. 

For  the  special  case  2  =  1,  the  longitudinal  boundary  condition  is  met  by  the  velocity 
components  of  an  incompressible  invisdd  fluid  at  the  surface  of  a  sphere  immersed  in  an 
infinite  uniform  flow,  or  the  magnetic  field  at  the  surface  of  a  spherical  superconductor  in 
an  otherwise  constant  field,  or,  in  short,  by  any  harmonic  field  whose  radial  component  is 
zero  at  the  surface  of  the  sphere  r  =  rm.  The  2  =  1  case  would  be  indistinguishable  from 
a  solid  sphere  with  no  internal  structure,  from  the  outside. 

In  the  next  sections,  arguments  will  be  made  to  show  that  a  conducting  plasma  may 
be  excited  to  sustain  a  non-radiating  electromagnetic  mode  whose  electric  and  magnetic 
fields  obey  the  force-free  conditions.  Arguments  will  also  be  made  to  show  that  these 
modes  can  be  stable  and  self-confining,  when  terminated  with  the  appropriate  boundary 
conditions.  The  current  in  the  plasma  will  be  forced  to  flow  along  the  streamlines  of  a 
force-free  vortex. 

Electromagnetic  volume  stresses  will  be  ignored  in  this  note,  but  are  not  negligible. 
These  stresses  can  be  counteracted  by  fluid  motion.  A  full  discussion  of  the  interaction  of 
fluid  stresses  and  electromagnetic  stresses  will  not  be  needed  to  understand  the  principles 
presented  here. 

In  the  final  section  and  in  the  Appendix,  some  speculations  will  be  made  as  to  how 
this  theory  can  be  applied  to  predict  the  existence  and  properties  of  stable  “particles”  that 
are  described  by  non-radiating  modes  trapped  in  spheres  of  neutral  plasma. 

m.  INTERACTION  WITH  ELECTROMAGNETIC  FIELDS 

Let  H  in  eq.(l)  be  a  magnetic  field.  If  E  is  the  accompanying  electric  field  and  J  is 
the  current,  in  otherwise  free  space  Maxell’s  equations  are 


VxH  =  eo^+J, 


(15) 


and 


(16) 


„  _  3H 

VxE  =  -/i0-^-, 


where  to  and  po  are  the  permittivity  and  permeability  of  the  vacuum. 

Assuming  a  time-harmonic  dependence  e_,"‘,  eqs.  (1),  (15)  and  (16)  imply  that 

E  =  -  Vrl>,  (17) 

where  rj>  is  a  scalar  potential.  The  vector  potential,  A,  is  then 


a 


Prom  the  Lorentz  condition 


(18) 


it  follows  that 


Sib 

V- A  +  «oA*o-^-  =  0, 


(19) 


^  =  -iu^=0,  (20) 

because  V-H  =  0.  The  above  choice  for  the  vector  potential  was  not  unique.  It  was  made 
to  force  the  current  to  also  have  the  form  of  a  force-free  field.  To  see  this,  solve  for  J: 


J  =  a(l  —  ~j)Hi  (21) 

;  ' 

where  k  is  the  free-space  wave  number  w^/iofio.  EVom  eq.(21)  it  follows  that  V-J  =  0, 
and  that  J  is  also  a  force-free  mode  if  H  is  a  force-free  mode.  This  implies,  also,  that 
at  any  frequency  greater  them  zero,  a  force-free-field  current  does  not  produce  a  scalar 
potential.  Therefore,  an  excess  of  charge  of  a  single  sign  entrained  in  a  force-free  mode 
in  a  time-harmonic  electromagnetic  field  will  have  a  non-time-harmonic  electrostatic  field 
associated  with  it.  In  other  words,  the  charge  density  in  a  stationary  non-radiating  mode 
may  not  be  uniform  in  space,  but  it  is  not  a  function  of  time.  It  will  be  shown  below  that 
the  charge  density  must  also  be  spatially  uniform. 


7 


Note  that  if  fc  =  ja|,  with  real  a,  the  current  becomes  zero  and  the  non-radiating 
mode  corresponds  to  a  vacuum  excitation  of  the  electromagnetic  field  having  finite  energy 
density  and  zero  mean  momentum  density.  (Momentum  density  is  proportional  to  E  x  H, 
whose  mean  is  zero  here.)  Note  also  that  the  relative  phase  of  the  current  and  the  magnetic 
field  it  generates  reverses  at  fc  =  |a|.  As  the  frequency,  u,  increases,  the  value  of  fc  and 
subsequently  the  absolute  volut  of  J  compared  to  H  increases,  beyond  the  frequency  when 
fc  =  |o|. 

Equation  (21)  has  some  interesting  consequences  when  the  interaction  energy  of  the 
current  and  the  magnetic  field  is  considered.  Following  the  logic  in  Stratton’s  book*13', 
but  accounting  for  time-harmonic  variation  variations  of  the  field,  an  expression  can  be 
written  for  the  r.m.s.  energy,  fi,  of  the  current,  J,  in  the  presence  of  the  vector  potential, 
A  (the  asterisk  *  means  complex  conjugation): 

G  =  \j3-A.‘dv.  (22) 

From  eqs.  (18)  and  (21),  this  becomes 

3  =  1 (23) 

For  real  values  of  a  the  magnetic  interaction  part  of  the  energy  density  becomes 
negative  at  high  frequencies  (w  >  ca,  c  is  the  speed  of  light).  The  implication  is  that 
force-free-mode  currents  created  in  short-wavelength  electromagnetic  fields,  residing  in 
local  minimum  of  potential  energy,  are  mechanically  stabilized  by  the  electromagnetic 
fields.  IVom  eq.(22),  the  higher  the  frequency  of  the  field  the  more  “bound”  the  current. 
In  order  for  the  above  analysis  to  apply,  assuming  its  interpretation  is  correct,  it  has  to  be 
possible  for  the  field  to  couple  to  the  current  and  the  boundary  conditions  have  to  be  met. 

In  a  way,  the  streamlines  of  a  force-free  vortex  seem  to  act  like  waveguides,  allowing 
modes  to  penetrate  if  they  are  small  enough  to  “fit”.  At  the  critical  frequency,  u>  =  ca, 
an  electromagnetic  free-space  mode  can  exactly  match  the  size  of  the  force-free  mode.  At 
higher  frequencies,  free-space  modes  can  fit  inside  the  force-free  mode  and  stabilize  it. 

IV.  ENERGY,  FORCE  AND  MOMENTUM  CONSIDERATIONS 

The  ultimate  practical  question  is,  “Can  this  theory  be  applied  to  produce  stably 
confined  plasma  and,  if  so,  how?”  First  consider  the  forces  on  individual  “blobs”  of  charge 
that  make  up  a  current.  The  Lorentz  force  on  a  charge  density,  p,  moving  at  a  velocity  v 
in  an  electromagnetic  field  is 
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F  =  p(E  +  /ioVxH). 


(24) 


Since  J  =  pv  this  reads 

F  =  pE  +  PoJxH  =  pE  (25) 

for  force-free  inodes.  In  other  words,  the  only  non-zero-time-averaged  force  acting  on  the 
charge  in  an  isolated  force-free  mode  arises  from  the  zero-frequency  electrostatic  field  due 
to  the  other  charges.  In  a  neutral  plasma  the  electrostatic  forces  are  evidently  zero.  If 
excess  charge  is  present  in  the  current,  an  electrostatic  field  will  perturb  the  trajectories. 
It  will  be  very  instructive  to  determine  the  sum  of  the  magnetic  interaction  energy  and  the 
kinetic  energy  to  determine  t  e  stability  of  the  overall  system.  The  reason  for  ignoring  the 
energy  in  the  oscillating  electric  field  is  that  it  does  not  directly  contribute  to  the  stability 
of  the  current.  I.e.,  an  infinitesimal  change  in  current  will  affect  the  magnetic  interaction 
energy  and  the  kinetic  energy  only. 

The  electric  field  will,  however,  determine  the  kinetic  energy  of  the  equilibrium  system. 
At  high  frequencies  the  motion  of  the  charges,  ignoring  electrostatic  and  outside  mechanical 
forces,  is  small-displacement  oscillation  about  an  equilibrium,  at  the  frequency  of  the  field. 
(It  is  important  to  realize  that  there  is  no  net  motion  of  the  charges.  The  background 
neutralizing  ions  are  also  assumed  to  be  motionless.)  The  motion  of  the  charges,  assumed 
to  be  carried  by  electrons,  is  driven  by  the  electric  field,  which  is  parallel  to  the  current. 
Calculations  will  all  be  non-relativistic. 

In  a  simple-minded  initial  intuitive  approach,  which  admittedly  should  not  be  giving 
the  correct  results,  Newton’s  law  F  =  ma  wili  be  applied.  The  volume  force  on  a  charge 
distribution  with  density  p  is  pE.  The  acceleration  is  the  time  derivative  of  current  divided 
by  the  charge  density,  p: 

pE  =  pm-,  (26) 

where  pm  is  the  mass  density  of  the  charge.  From  eqs.  (17,i,  (20)  and  (21),  putting  all 
quantities  in  terms  of  E  and  rearranging, 


Cancelling  the  iu>  and  assuming  that  all  the  current  will  be  carried  by  electrons,  the 
expression  for  the  electron-cloud  charge  density  is 
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(28) 


I  I  ( "A  (*J  “  0,J)  L. 

h“(t y—^r1'  *»”• 

where  m/e  =  5.69x10“ 12  Kj/cmil  is  the  mass-to-charge  ratio  of  the  electron.  An  important 
note  is  that  eq,(28)  is  only  valid  for  k  >  a.  The  reason  for  this  is  that  the  sign  of  the 
charge-density  distribution  changes  at  k  —  a  and  for  k  <  a  the  sign  is  inconsistent  with 
negative  charge  carriers.  This  is  an  important  factor  in  stability  considerations*11'. 

Substituting  for  quantities  to  obtain  the  r.m.j.  kinetic  energy  density, 

l^=i|(l-§)h|Hi2-  <29) 

Newton’s  law,  in  its  simple-minded  form  applied  above,  is  not  strictly  valid  in  this 
case.  A  constant  charge  and  mass  density  was  implicitly  assumed  from  the  outset  if 
Newton’s  law  is  applied.  The  fact  that  a  spatially  uniform  charge  density  results  shows 
that  the  assumption  leads  to  self-consistent  results. 

In  the  Appendix,  the  lagrangian  formulation  of  the  system  will  be  investigated.  It 
will  be  shown  that  the  result  in  eq.(28)  can  be  derived  without  assuming  uniform  charge 
density.  It  is  valid  under  all  conditions  where  overall  charge  neutrality  is  maintained  and 
gas-dynamic  effects  are  ignored. 

Equations  (28)  and  (29)  are  extremely  important  from  a  practical  point  of  view.  Prom 
eq.(28)  it  is  evident  that  the  amount  of  charge  participating  in  the  motion  is  determined 
only  by  a  and  the  frequency,  furthermore,  as  proved  in  the  Appendix,  uniform  volume 
charge  density  is  a  mandated,  variationally  correct  solution  of  the  lagrangian  equations. 
An  interesting  point  is  that  the  free-charge  density  is  independent  of  the  amplitude  of  any 
electromagnetic  field  quantity.  By  assuming  that  an  Ohm’s  law  can  be  written  in  the  form 
J  =  <r(w)E  for  a  non-radiating  mode,  then  the  conductivity  o(u>)  is  purely  imaginary  and 
can  be  written 


<r(w)  = 


i(k2  -  a3) 

UflQ  ’ 


(30) 


so  that  at  high  frequencies  a  non-radiating  mode  becomes  purely  capacitive,  with  an 
asymptotic  specific  capacitance  equal  to  to  =  8.S5p//m.  The  endeavor  to  drive  modes  at 
higher  and  higher  frequencies  might  encounter  a  difficulty  equivalent  to  trying  to  match 
a  very  small  impedance.  Very  high  field  strengtns  might  be  required  in  small  (order  of 
ten  wavelengths  on-a-side)  volumes.  Laser  radiation  may  be  a  practical  solution  as  a 
high-power  driver. 
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From  eqs.(23)  and  (29)  the  sum  of  the  kinetic  and  potential  energies  relevant  to 
stability  of  the  overall  motion  is  zero,  if  surface  stabilizing  currents  and  vortical  fluid 
motion  are  neglected.  (Note,  however,  that  the  total  energy  of  the  system  is  positive.)  As 
the  frequency  increases  beyond  the  critical  frequency,  the  amount  of  charge  contributing 
to  the  current  increases,  the  kinetic  energy  increases  and  the  potential  energy  decreases 
for  constant  absolute-value  field  strengths.  Beyond  the  critical  wavenumber,  the  kinetic 
energy  increases  asymptotically  like  the  square  of  the  frequency. 

For  energy  and  momentum  balance,  a  critical  amount  of  free  charge  must  be  available. 
If  not  enough,  or  too  much  charge  is  present  the  conditions  of  the  motion  may  not  be  met 
at  the  given  frequency.  In  addition,  the  system  cannot  radiate  from  its  interior  because 
the  E  and  H  vectors  are  parallel  everywhere  internally. 

Evidently,  there  is  no  restriction  on  the  amplitude  of  the  electromagnetic  fields  within 
the  mode.  Since  the  theory  has  been  derived  here  for  non-relativistic  plasmas,  there  may 
be  some  eventual  dependence  on  field  amplitudes  in  ultra  energetic  modes.  For  very  small 
field  amplitudes,  fluid  and  thermal  forces  might  disrupt  a  non-radiating  mode.  When 
real  systems  are  considered,  massive  electrons,  which  carry  the  current,  cannot  constrain 
their  motion  to  the  curved  electric  field  lines.  This  introduces  a  contraint  on  the  field 
amplitudes,  such  that  the  exclusions  of  the  electron  positions  from  the  field  lines  must  be 
“small’’.  ‘  ' 

Field  curvature  leads  to  electromagnetic  stresses  that  induce  electronic  “drift”  across 
field  lines.  The  interaction  of  the  resultant  pressure  tensor  on  the  free-electron  gas  with 
the  overall  fluid  is  very  important.  None  of  these  complications,  however,  will  pursued 
here. 

The  consequences  of  this  approximation  will  be  discussed  in  a  sequel^11*. 

In  what  follows,  it  is  assumed  that  a  finite  spherical  region  of  a  short-wavelength  non¬ 
radiating  mode  is  sustained  by  matching  only  the  electromagnetic  boundary  conditions  at 
a  fixed  radius.  As  was  mentioned  above,  modes  can  be  chosen  such  that  the  boundary 
conditions  reduce  to  only  tangential  conditions  on  the  velocities  and  the  fields.  This  can 
translate  into  a  requirement  for  surface  currents  to  produce  jumps  in  the  electric  and 
magnetic  fields  going  from  the  force-free  fields  inside  the  spherical  region  to  the  free-space 
fields  outside,  and/or  vortical  fluid  motion  that  can  counteract  the  electromagnetic  surface 
and  internal  stresses. 

No  mechanism  will  be  described  to  fully  invoke  the  needed  conditions.  They  may 
be  very  difficult  to  reproduce  physically.  If  such  boundary  conditions  were  not  provided, 
however,  the  non-radiating  mode  would  decay.  The  speed  of  this  decay  has  not  been 
calculated,  but  it  seems  intuitively  reasonable  that  systems  with  large  densities  of  kinetic 
energy  would  take  longer  to  decay  than  those  with  small  densities,  if  fluid-dynamic  effects 
are  ignored.  Surface  currents  and  vortical  motion  serve  the  purpose  of  counteracting  the 
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normal  electromagnetic  stresses,  which  would  otherwise  tend  to  force  the  excitation  to 
expand. 

Zero  total  energy  of  the  free-electron  cloud,  above  the  critical  frequency,  means  that 
the  electronic  motion  is  in  neutral  stability,  but  might  show  some  rigidity  in  response  to 
compressive  forces.  Unless  “other”  mechanical  forces,  such  as  internal  vortical  motion, 
are  involved,  a  non-neutral,  classical  force-free  system  will  self-destruct  due  to  internal 
electrostatic  forces.  An  ionized  fluid  or  neutral  plasma  will  barely  be  stable,  if  mechanisms 
are  not  present  that  can  lower  the  overall  energy.  Surface  currents  contribute  to  stability 
by  a  “pinch”  mechanism,  which  will  effectively  be  a  surface  tension. 

It  must  be  emphasized  that  lowered,  energy  due  to  either  the  pinch  effect  of  currents 
or  vortical  motion  are  absolutely  required  for  stability.  What  makes  the  present  model 
attractive  is  the  minimization  of  the  body  forces,  induced  by  the  interaction  of  the  current 
with  the  force-free  electromagnetic  field.  By  themselves,  these  body  forces  could  tend  to 
severely  limit  the  lifetime  of  an  excitation. 

The  question  of  meeting  the  boundary  conditions  for  the  electromagnetic  field  also  has 
to  be  addressed  carefully.  A  non-radiating  mode  with  a  given  harmonic  internal  current 
will  generate  external  harmonic  electric  and  magnetic  fields  that  do  not  have  the  force-free 
form.  It  is  a  simple  matter  to  say  that  the  “conect”  boundary  conditions  on  the  field,  i.  e., 
those  that  do  not  permit  radiation  to  escape  from  a  truncated  volume,  can  be  reproduced 
with  surface  densities  of  electric  and  magnetic  current*16).  There  are  difficulties  with  such 
a  solution.  The  first  difficulty  is  that  a  surface  density  of  magnetic  current  may  not  be 
achievable.  Such  a  current  may  be  approximated,  however,  by  “eddies”  at  the  surface. 
The  pure  electric-current  boundary  condition,  on  the  other  hand,  can  be  met  by  having 
a  physically  achievable  surface-current  distribution.  In  either  case  large  accumulations  of 
surface  charge  would  be  required  to  produce  boundary  currents.  The  stabilizing  effects 
of  such  currents  can  also  be  achieved  by  setting  the  supporting  fluid  into  a  vortical  flow, 
counteracting  the  electromagnetic  stresses. 

This  point  will  be  diacussed  further  in  a  sequel,  where  it  will  be  shown  that  the  proper 
choice  of  an  excitation  radius  for  a  given  frequency  can  drastically  affect  the  radiation 
boundary  conditions  and  stability. 

Another  way  of  looking  at  the  overall  picture  is  to  realize  that  in  a  harmonic  field, 
the  electric  field  can  have  a  minimum.  J.e.,  the  first  derivatives  of  the  potentials  can  all 
be  zero,  and  their  second  derivatives  can  all  be  positive.  Unlike  the  physical  situation 
giving  rise  to  Earoshaw’s  theorem*1  J\  the  components  of  the  potentials  and  fields  obey 
the  Helmholtz  equation  instead  of  the  Laplace  equation.  This  allows  the  possibility  of 
stabilization  of  charges  and  currents  under  the  influence  of  only  self-consistent  harmonic 
electromagnetic  fields  and,  perhaps,  added  vortical  fluid  motion. 

In  short,  boundary  conditions  for  truncated  force-free  modes  that  are  met  by  surface 
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currents  and/or  internal  fluid  motion  may  lead  to  plasma  excitations  that  are  mechanically 
stable  and  relatively  long  lived. 

Puthoff  has  suggested^17'  that  by  introducing  the  zero-point-field  (ZPF),  stability 
may  be  greatly  enhanced.  Conducting  systems  will  experience  a  net  outside  pressure, 
analogous  to  the  Ca3imir  force,  due  to  the  exclusion  of  ZPF  modes.  The  more  conducting, 
at  higher  frequencies,  the  more  ZPF  modes  will  be  excluded  from  the  interior,  and  the 
higher  the  pressure  differential.  Additionally,  the  energy  to  stabilize  non-radiating  modes 
might  be  available  from  the  higher  frequencies  in  the  ZPF,  once  the  mode  has  been  formed. 
If  the  above  theory  is  to  be  believed,  The  ZPF  may  play  an  important  role  in  stabilizing 
non-radiating  modes. 

Detailed  calculations  would  have  to  be  done  to  determine  critical  densities  of  charge 
versus  size  and  frequency  for  all  the  quantities  to  come  together  exactly  right.  Further 
research  would  investigate  how  to  inject  charge  into  already-formed  modes,  and  have  them 
stabilized  by  higher  frequencies  that  come  from  directed  sources,  or  are  present  in  the  ZPF. 

A  final  point  relates  to  the  initial  construction  of  non-radiating  modes.  If  these 
modes  are  stable,  why  are  they  not  a  common  phenomenon?  Actually,  they  might  be 
quite  common  in  Nature,  but  with  short  lifetimes  and  low  energies.  They  also  might  be 
very  small  and  hard  to  detect.  Considering  that  the  most  logical  place  to  look  for  charged, 
high-frequency  non-radiating  modes  would  be  in  powerful  electrical  discharges,  such  as 
found  in  lightning  strokes,  the  surrounding  phenomena  might  mask  their  occurauce. 

In  the  laboratory,  strong  microwave  or  laser  (perhaps  accompanied  with  D.C.)  dis¬ 
charges,  with  possible  extra-electron  injection,  might  be  conducive  to  non-radiative-mode 
production.  Discharges  in  a  highly-ionizable  gas  such  as  cesium  or  Penning  mixtures,  un¬ 
der  the  above  conditions  may  enhance  the  effect.  In  any  case,  detailed  calculations  and 
experiments  might  have  to  be  done  to  determine  optimal  conditions  such  as  discharge  field 
strengths,  gas  density,  frequency,  electrode  size,  etc. 

The  investigation  of  non-radiative-mode  effects  may  be  very  fruitful  in  the  future. 
For  example:  Rigid  non-radiative  composite  “particles”  containing  large  numbers  of  fus- 
able  ions  may  be  compressible  in  strong  quadrupole  fields  without  undergoing  chaotic 
turbulence-like  instabilities,  and  have  lifetimes  consistent  with  producing  copious  ther¬ 
monuclear  reactions.  Large-diameter  mooes,  with  significant  lifetimes,  might  be  capable 
of  storing  large  amounts  of  energy,  etc.  The  field  seems  to  be  wide  open,  hinging  on  simple 
experiments  that  may  prove  the  basic  principles. 
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APPENDIX 


The  assumption  that  a  self-consistent  charge  distribution  (neutralized  by  charges 
of  the  opposite  polarity)  can  carry  harmonic  current  in  a  non-radiating  mode  leads  to 
equations  for  the  charge  density  at  any  point.  The  lagrangian  equations  for  the  system 
lead  to  a  non-linear  first-order  vector  partial  differential  equation  for  the  charge  density.  It 
trill  be  shown  that  an  exact  solution  exists  wherein  the  charge  density  is  certainly  uniform. 
Equation  (28)  will  be  found  to  be  valid  at  all  field  strengths. 

It  trill  also  be  shown  that  the  occurrence  of  such  modes  is  consistent  with  a  uniform 
mass  density  and  rigid  motion  in  a  plasma  or  gas  whose  polarizability  is  a  function  of  only 
the  gas  density. 

The  r.m.3.  lagrangian  density^18'  for  a  current  of  massive  particles  in  the  presence  of 
a  non-radiating  harmonic  electromagnetic  field  in  a  neutral  plasma  is  proportional  to 

C  =  \Pm  |*|J  +  J-A’  +  i(eo|E|J  -  mo|H|2)  -  p^,  (Al) 

where  pm  is  the  mass  density  of  the  current,  r  is  the  current  field’s  spatial  coordinates  (the 
dot  means  velocities),  pe  is  the  charge  density,  and  the  scalar  electromagnetic  potential, 
0,  from  eq.(20),  is  zero. 

Using  the  relation  J  =  pt r,  with  equations  (1),  (16),  (18)  and  (20),  the  lagrangian 
equations  for  the  motion  can  be  written  in  the  vector  form,  after  a  little  algebra. 


where  m/e  =  5.69  x  IQ~12 Kg /coul  is  the  electron  mass-to-charge  ratio. 
Equation  (A2)  immediately  yields  one  solution  by  inspection: 


(A3a) 
(A3i) 

which  can  be  solved  for  the  same  p,  given  by  eq.(28): 


Vp,  =  0, 


and 


(5K‘-5)+<p-' 
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*>a’  (A4) 

It  is  not  clear  that  eqs.(A3)  represent  a  unique  exact  solution  for  h?  /  a2.  The 
lagrangian  equations  may  be  integrable  for  non-uniform  density  distributions.  Such  a  case 
might  arise  for  a  shape  perturbed  away  from  spherical.  Under  that  assumption,  it  would 
also  be  expected  that  an  entire  mode  might  oscillate  about  an  equilibrium  with  a  frequency 
that  might  be  much  lower  than  the  frequency  of  the  electromagnetic  excitation,  u.  This 
possibility  was  not  taken  into  account  in  setting  up  eq.(A2),  where  it  was  assumed  that 
the  partial  time  derivatives  of  the  mass  and/or  charge  densities  are  zero.  The  total  time 
derivatives  of  these  quantities  that  appear  in  the  lagrangian  equations  arise  only  from  the 
term  r  •  Vpt. 


It  can  now  be  shown  that  the  above  analysis  is  consistent  with  a  null  virial.  The 
analysis  below  will  be  expanded  somewhat  in  a  subsequent  note*11*. 

The  particles  in  a  high  frequency  non-radiating  plasma  mode  will  all  be  assumed  to 
perform  small-amplitude  oscillations  about  equilibria,  so  that  at  any  time  the  position, 
s(r,<),  of  any  particle  can  be  described  by 

«(r,i)  =  r+f  /'j(r)e-^'<fr',  (AS) 

Pt  J  s 

where  J(r)e~'"‘  is  the  time-dependent  current  at  the  space-fixed  point  r.  The  force, 
F(r,t)  acting  on  the  particle  comes  from  the  electric  field, 

F(r,t)  =  *E(r)e-"«.  (A6) 


The  virial  theorem* 1  ')  requires  that 


Jv  t(r,  t)  •  F*(r,  t)dV  =  -PmJv  |v|W. 
By  straightforward  algebra  it  is  evident  that 


(■47) 


/v«(r,0  •  F*(r,  t)dV  =  (i)  jf  J(r)  •  E*(r )dV  =  jf  |J(r)|JdV,  (A8) 
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and 


pm  /  jvfdV  =  f  |J(r)l*dVi  (AO) 

Jv  P*  Jv 

where  the  integrations  are  over  the  volume  of  the  mode.  For  purely  spherical  modes  the 
“natural”  volume  for  integration  might  be  from  the  origin  to  the  first  zero  of  the  appropriate 
Bessel  function  associated  with  the  mode. 

Substitution  of  eqs.(A8)  and  (A9)  into  (A7)  gives: 


f  |J{r)|2dV  =  ^  (  P(r)|adV.  (A10) 

k 2  -  a2  Jv  P •  Jv 

Cancelling  the  integral  from  both  sides  of  eq.(AlO),  and  solving  for  pt,  obtains 


Pe 


m\(lW) 
.  e  /  #*o 


(All) 


which  is  identical  to  the  result  assumed  from  a  lagrangian  formulation,  as  well  as  the  simple 
analysis  based  on  simple  mechanics  preceding  eq.(28). 

Note  that  eq.(All)  needs  the  relationship  between  the  mass  density  of  the  electrons 
and  the  charge  density,  i.e.,  pm/p.  =  -m/e  =  -5.69xlO~uKg/coul,  where  the  minussign 
arises  because  the  charge  density,  p„  is  negative.  The  somewhat  interesting  insight  from 
eq.(All)  is  that  unless  the  frequency  is  high  enough,  k2  >  a2,  the  required  charge  density 
for  the  electrons  must  be  positive,  dearly  an  inconsistency.  Not  only  must  the  frequency 
be  high  enough  to  predict  a  consistent  charge  density  for  the  current  carriers  (the  much 
heavier  ions  are  never  considered  capable  of  carrying  the  high-frequency  current),  a  full 
virial  analysis^11*  shows  that  one  of  the  conditions  for  stability  is  that  k2>a  . 

The  question  can  also  arise  as  to  what  the  allowable  gas  density  can  be  in  a  non¬ 
radiating  mode. 

Without  taking  compressible  fluid  dynamics  and  its  interaction  with  the  electromag¬ 
netic  stresses  into  account,  the  density  cannot  be  fully  analyzed.  Further  work,  to  be 
reported  on  in  more  detail  at  a  later  time,  seems  to  indicate  that  when  vortical  motion  is 
introduced,  a  stable  exatation  can  arise  that  has  uniform  density  and  entropy,  with  a  gas 
pressure  exactly  equal  to  the  gas  pressure  of  the  surrounding  atmosphere  at  the  time  of 
its  creation,  and  is  stable  against  small  perturbations  of  internal  veloaty(11). 
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la  summitry,  the  assumption  that  non-radiating  modes  can  exist  in  a  neutral  plasma 
leads  to  the  conclusion  that  the  tree-charge  density  entrained  by  the  current  can  be  a 
constant  throughout  the  mode  and  dependent  on  the  tree-space  electromagnetic  wave 
number  tor  a  given  modal  a.  It  is  also  consistent  with  assuming  that  the  only  forces  acting 
in  the  volume  arise  from  fluid  motions  and  the  time-varying  electromagnetic  field,  in  a 
neutral  plasma. 
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ABSTRACT 

Conditions  for  a  null  virial  for  electrons  forming  the  current  in  a 
force-free  harmonic  mode  are  investigated.  Above  a  critical  fre¬ 
quency,  this  current  lies  in  a  potential  well  whose  relative  depth  in¬ 
creases  as  the  frequency  increases.  Conditions  for  greatest  stability, 
ignoring  electromagnetic  stresses,  are  investigated.  Fluid-dynamic 
effects  are  also  ignored.  Under  appropriate  approximations,  the  con¬ 
ditions  for  a  null  electron  virial  turn  out  to  be  identical  to  conditions 
for  the  motion  of  electrons  oscillating  along  the  electric  field  lines  in 
the  force-free  electromagnetic  field,  above  the  critical  frequency.  The 
-  free-electron  density  must  be  uniform,  and  depends  only  on  the  size 
of  an  excitation  and  the  frequency  of  the  time-harmonic  electromag¬ 
netic  field.  Assuming  reasonable  limits  on  physical  parameters  leads 
to  upper  bounds  on  the  energy  content  of  excitations  in  the  mega- 
joule  range  for  a  sphere  with  a  10cm  radius.  A  negative  potential 
well  for  electrons,  above  the  critical  frequency,  allows  the  possibil¬ 
ity  of  quantum  effectk,  which  are  not  investigated.  Further  areas  of 
research  are  briefly  discussed. 

L  INTRODUCTION 

If  the  vector  field  H(r)  is  a  Beltrami,  or  force-free  field'1’2*  then 


VxH  =  oH, 


(1) 


where  a  can  be  a  complex  constant.  Let  H  in  eq.(l)  be  a  magnetic  field.  If  E  is  the 
accompanying  electric  field  and  J  is  the  current,  in  otherwise  free  space  Maxwell’s  equations 
are 


SE 

VxH  =  ec^+J, 


(2) 


1 


VxE  =  -^,  (3) 

where  e0  and  /j0  are  the  permittivity  and  permeability  of  the  vacuum. 

Assuming  a  time-harmonic  dependence  e~'“l  for  all  quantities,  eqs.  (1),  (2)  and  (3) 
imply  that 


a 

where  ij>  is  a  scalar  potential.  The  vector  potential,  A,  is  then 


(4) 


a 


(5) 


Prom  the  Lorentz  condition 


dib 

V- A  +  =  0, 


(6) 


it  follows  that 


^  =  -M-  =  0,  (7) 

because  V-H  =  0.  The  above  choice  for  the  vector  potential  was  not  unique.  It  was  made 
to  force  the  current  to  also  have  the  form  of  a  force-free  field.  To  see  this,  solve  for  J: 


J  =  a(l  —  “j)H,  (8) 

where  k  is  the  free-space  wave  number  wv/ij/ IJ.  From  eq.(l)  it  follows  that  V-J  =  0,  and 
that  J  is  also  a  force-free  mode  if  H  is  a  force-free  mode.  EYom  eqs.(4)  and  (8),  E  can 
also  be  written  in  terms  of  J, 


2 


E  = 


imp o  , 
a2  -  k2 


(9) 


It  is  important  to  note  that  eqs.(l-9)  will  not  all  necessarily  hold  if  a  has  any  spatial 
variation.  It  also  can  be  shown  that,  starting  with  a  force-free  harmonic  current  that  fills 
all  of  space,  the  electric  and  magnetic  fields  will  be  harmonic  and  also  have  the  force-free 
form,  if  a  is  constant. 


A  feature  of  the  above  equations  is  that,  unlike  in  a  vacuum  plane-wave,  as  the 
frequency  increases,  the  field  energy  becomes  more  concentrated  in  the  electric  field.  For 
the  direct  contribution  of  the  magnetic  field  energy  to  the  total  energy  becomes 
negligible  compared  to  the  contribution  from  the  electric  field.  Indeed,  for  k>  a,  with 
critical  frequency  u>cr,(  =  ac,  the  interaction  energy  of  the  current  and  the  magnetic  field 
is  negative.  In  other  words,  the  current  finds  itself  in  a  potential  well,  whose  relative  depth 
increases  with  frequency,  all  other  things  Doing  equal,  after  the  point  where  k  —  a.  This 
energy  well  is  intriguing,  considering  the  quantum  mechanical  possibilities  and  stationary 
states  that  may  arise.  It  was  already  shown(2)  that,  classically,  the  sum  of  the  magnetic 
interaction  energy  of  the  free  electrons  and  their  kinetic  <  -ergy  is  zero.  Although  quantum 
and  relativisitc  effects  may  alter  the  picture  somewhat,  :hey  will  not  be  considered  here. 

The  above  equations  represent  a  valid  solution  of  Maxwell’s  equations  that  can  fill 
all  of  space  with  massless  current.  A  more  realistic  physical  situation  realizing  finite- 
size  excitations  and  non-zero-mass  current  carriers  (electrons)  will  now  be  considered. 
Toward  that  end  it  will  be  shown,  using  the  virial,  that  the  stability  of  finite-size  harmonic 
excitations  in  a  neutral  fluid  be  optimized  by  adjusting  the  free-charge  density.  The 
optimum  free-charge  density  will  be  shown  to  depend  on  the  electromagnetic  wavenumber, 
} ,  the  size  parameter  of  the  excitation,  a,  and  the  condition  k>a. 

A  number  of  assumptions  will  be  involved  in  the  course  of  performing  the  analysis. 
The  first  assumption  is  that  the  background  ions  do  not  participate  in  the  motion  of  the 
excitations,  nor  are  they  affected  by  any  of  the  high-frequency  fields.  It  is  next  assumed 
that  all  the  temperatu*  es  are  low,  and  the  gas-dynamic  pressures  inside  and  outside  the 
excitations  are  equal,  so  that  thermodynamic  and  fluid-dynamic  effects  will  be  ignored. 
In  other  words,  to  first  approximation  the  excitations  to  be  studied  will  “float”  in  an 
ambient  fluid  with  no  net  motion.  The  gas  density  and  free-electron  density  will  also  be 
assumed  uniform  in  space  within  the  excitation.  A  major  result  of  this  note  will  be  the 
determination  of  the  optimal  uniform  free-charge  density. 

The  next  assumption  involves  an  approximation.  It  will  be  assumed  that  the  elec¬ 
trons,  carrying  all  the  current,  oscillate  harmonically  along  curved  electric  field  lines.  This 
approximation  would  be  exact  if  the  electrons  had  no  mass  or  if  the  amplitudes  of  os¬ 
cillation  were  truly  infinitesimal  and  the  temperature  was  low  enough.  As  the  frequency 
increases  and  the  field  amplitudes  decrease,  this  approximation  becomes  more  accurate. 
Excursion  of  electronic  motion  from  the  curved  electric  field  lines  forces  the  current  to 
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have  non-force-free  components  and  can  thereby  lead  to  losses  via  synchrotron  radiation, 
so  it  is  important  that  these  excursions  be  minimized.  In  addition,  internal  pressure  gra¬ 
dients  can  arise  if  the  field  lines  me  curved.  These  pressure  gradients,  proportional  to  the 
square  of  the  field  intensities,  may  become  so  large  as  to  dominate  any  stability  analysis. 
Nevertheless,  such  effects  will  be  ignored  in  this  note. 

At  the  outer  surface  of  the  excitation,  it  will  be  assumed  that  the  jumps  in  the 
electric  and  magnetic  harmonic  fields,  from  large  values  inside  to  zero  outside,  sustain 
equivalent  surface  currents  t5’3*.  Since  these  harmonic  surface  currents  may  be  very  large, 
to  accomodate  their  size  in  a  thin  layer  without  large  excursions  of  the  electrons  might 
require  very  large  concentrations  of  (neutralized)  free  charge  near  the  surface,  and  therefore 
large  concentrations  of  mass  at  the  surface.  These  will  be  assumed  to  exist.  Surface 
currents  are  invoked  to  prevent  radiation  losses.  Subsequent  work  will  deal  with  surface 
radiation  in  a  more  realistic  manner.  To  do  so  here  would  not  be  appropriate. 

Conditions  under  which  all  the  above  assumptions  can  be  met  in  a  reproducible  ex¬ 
perimental  situation  will  not  be  handled  in  this  note,  but  none  of  the  assumptions  appear 
to  violate  any  physical  laws  or  introduce  contradictions.  Subsequent  work  is  being  directed 
toward  eliminating  the  need  for  surface  currents,  and  using  vortical  fluid  modes  to  counter 
electromagnetic  stresses. 

II.  VIRIAL  CONSIDERATIONS 

Consider  the  mean-square  moment  of  inertia,  I,  of  the  electrons  in  a  finite-volume 
excitation,  and  its  time  derivatives: 

I  =  ^pmr*.rdV,  (10) 

=  JvPm[r’-r  +  r'-r]dV,  (11a) 

and 

^  =  /=  J%m[r’-r  +  2|if +  r’.r]dV,  (116) 

where  pm  is  the  (spatially  uniform)  mass  density  of  the  electrons,  r  is  the  complex  position 
vector  associated  with  an  “infinitesimal"  blob  of  electrons,  and  the  integration  is  over  the 
volume  of  the  excitation.  Strictly  speaking,  the  integrations  in  eqs.(10)  and  (11)  should 
be  sums  over  individual  electrons,  but  their  densities  are  assumed  high  enough  for  the 
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integral  to  be  an  appropriate  approximation.  The  quantity  I  is  will  be  called  the  “virial  of 
the  electronic  motion”*'* ,5'8'7*.  The  exclusion  of  the  electromagnetic  field  from  the  above 
integrals  will  be  discussed  below. 

The  quantity  r  =  r(x,y,z,t)  a  the  instantaneous  position  of  an  infinitesimal  volume 
whose  mean  position  is  at  the  fixed  point  ro  =  r (x,y,z,t  =  to),  and  oscillates  with  the 
same  frequency  as  the  applied  harmonic  field.  By  definition,  the  velocity,  v,  is  proportional 
to  J: 


v  =  3/Pt, 


(12) 


where  pt  is  the  (spatially  uniform)  charge  density  of  free  electrons.  The  charge  den¬ 
sity,  pe,  is  simply  related  to  the  mass  density,  pm  by  pt  =  (e/m)pm,  where  (e/m)  = 
1.759  x  10’-1  coul/kg.  Since  all  the  quantities  are  assumed  to  be  harmonic  in  time,  so  that 
3(x,y,z,t)  =  3(x,y,z)e~'ut,  then  r  can  be  written 

e-^dt  =  r0  +  (,^-^Je-“('-'o).  (13) 

A  sufficient  condition  for  self-confinement  is  that  the  time-integrated  average  of  I  is 
zero.  This  condition  implies  that  the  moment  of  inertia  of  the  electron  cloud  can  assume 
its  original  value  any  number  of  times,  and  cannot  become  uniformly  unbounded  as  time 
increases.  Instead  of  working  directly  with  I,  however,  investigation  will  be  made  into 
conditions  under  which  the  time-integrated  average  of  I  is  zero.  This  condition  is  more 
subtle.  A  necessary  condition  for  the  time-integrated  average  of  /  to  be  zero  is  that  the 
time-integrated  average  of  I  be  zero.  What  will  be  shown  is  that  if  the  uniform  charge 
density  has  exactly  the  right  size,  and  the  frequency  is  above  a  critical  value,  the  time- 
integrated  average  of  /  is  exactly  zero.  A  larger  value  for  the  charge  density  would  drive  the 
time-integrated  average  of  I  negative,  leading  to  collapse  and,  perhaps,  instabilities.  This 
is  because  the  exact  valu.  of  uniform  change  density  derived  for  a  null  is  the  same  required 
from  classical  mechanics  to  simultaneously  satisfy  Maxwell’s  equations  and  the  laws  of 
inertia,  with  finite-mass  electrons  in  harmonic  motion  along  field  lines.  For  this  reason, 
the  conditions  to  be  derived  will  be  considered  sufficient  for  self  containment.  Gerjouy  and 
Stabler*8*  also  discuss  this  point. 

Substituting  the  time-dependent  r  from  eq.(13)  into  eq.(llb)  and  taking  the  time 
average  automatically  gives  zero,  because  of  the  assumption  that  the  charges  oscillate  with 
fixed  centers  implies  confinement.  It  is  far  more  interesting  to  substitute  for  tne  terms  in 
I,  assuming  harmonic  motion  for  the  charges,  to  calculate  conditions  for  stability.  Denote 
the  time-integrated  average  of  a  quantity  by  a  bar,  i.e., 


r  =  r0  + 


Ll 

PcJ 
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Q  =  lim 

T  — »oo 


(14) 


To  obtain  new  information  from  eqs.(ll),  note  that  the  terms  in  the  integrand  con¬ 
taining  ?,  the  acceleration,  are  related  to  the  volume  force  on  the  electrons4*', 


Pm  r  =  fit  E(r)e-i"*.. 


The  term  with  |r|*  can  be  written 


(15a) 


(156) 


From  eqs.(9),  (11),  (13),  (14),  and  (15),  with  a  little  algebra. 


d‘I 

dt* 


Pm 

A 


+ 


Pt> 

a*  -  ** 


dV. 


(16) 


In  order  for  the  right  hand  side  of  eq.(16)  to  be  zero  the  following  conditions  must 
hold: 


k 


(17q) 


and 


where  (m/e)  =  5.69xl0”l*tj/cou/  is  the  mass-to-charge  ratio  for  the  electron.  Equation 
(17b)  was  derived  previously1"  using  a  lagrangian  formulation.  From  variational  con¬ 
siderations  introduced  in  classical  mechanics4*'8',  the  result  in  eq.(17b)  arises  from  only 
assuming  that  the  quantities  E,  H,  and  J,  have  the  force-free  form,  and  that  the  free  elec¬ 
trons  carrying  the  current  have  non-zero  mass.  Uniform  charge  density  is  not  an  a-priori 
assumption  in  the  lagrangian  calculation. 

A  diversion  is  needed  at  this  point  to  discuss  the  effects  of  the  harmonic  electromag¬ 
netic  field  assumed  in  the  present  model. 
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The  usual  analysis  involving  the  virial  theorem*4’!l’s,7^  includes  the  electromagnetic 
energy  density  in  the  integrand  in  eq.(10).  The  assumptions  made  here  are  that  the 
electronic  velocities  are  non-relativistic  and  the  harmonic  fields  E  and  H  are  non-zero 
inside  the  excitation  end  zero  outside,  which  is  equivalent  to  effective  electric  and  magnetic 
currents  respectively  equal  to 


J£//  =  -4xH,  (18a) 

and 

Me//  =  nxE,  (186) 


where  n  is  the  outward-pointing  normal  to  a  sharply  defined  excitation  boundary. 

The  first  thing  to  realize  is  that  if  such  currents  existed,  they  would  exactly  counteract 
the  electromagnetic  stresses  at  the  excitation  boundary.  The  reasoning  is  that  internally, 
the  electromagnetic  field,  terminated  by  the  currents  in  eqs.(18a-b),  seems  to  stretch  to 
infinity  with  a  zero  Poynting  vector  (ExH  =  0).  No  energy  is  lost  to  the  “outside”  if  the 
currents  are  maintained.  For  the  lowest  1  =  1  spherical  model2!,  aj  a  —  4.493 . . ,,  both 
E  and  H  are  tangential  to  the  boundary  and  therefore  the  stresses  due  to  the  internal 
fields  are  normal  to  the  surface.  The  oecurance  of  the  currents  in  eqs.(18a-b)  can  be 
thought  of  as  representing  pinches,  which  barely  contain  the  normal  stresses  due  to  the 
electromagnetic  field  trapped  inside.  Viewed  in  this  manner,  the  surface  currents  represent 
a  lowering  of  the  energy  of  the  overall  system. 

It  is  understood  that  surface  currents,  while  countering  the  electromagnetic  surface 
stresses,  will  not  null  the  pressure  gradients  that  arise  from  variations  in  the  fields,  and 
will  tend  to  produce  internal  motion.  No  mechanism  can  produce  ideal  surface  currents. 
The  present  work  is  primarily  aimed  at  finding  the  most  stable  configuration  of  the  charge 
distribution  when  electromagnetic  stresses  and  radiation  are  ignored.  Ongoing  work  ap¬ 
pears  to  indicate  that  judicious  choices  of  k  and  a  can  null  the  outgoing  radiation  field 
and  vortical  fluid  modes  can  null  the  electromagnetic  stresses.  Invoking  idealized  surface 
electric  and  magnetic  currents  will  merely  be  an  a  dodge  that  is  temporarily  needed  to 
obtain  an  important  result  concerning  the  charge  density.  These  surface  currents  will  not 
be  discussed  any  further. 

Returning  to  the  results  in  eqs.(17),  it  might  first  appear  that  if  the  free-charge  density 
were  increased  beyond  the  value  prescribed  in  (1Tb),  stability  would  increase  because  /,  in 
eq,(16)  would  decrease.  If  this  were  to  happen,  the  electromagnetic  fields  and  mechanical 
trajectories  would  no  longer  follow  force-free  field  lines  and  would  no  longer  be  consistent 
with  the  model.  Ftam  the  simplest  considerations,  is  was  shown  that  small  perturbations 
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in  fluid  density  did  not  grow*3).  It  will  be  assumed  that  if  the  conditions  already  discussed 
are  met,  with  the  currents  lying  in  a  negative-potential  well,  the  excitations  under  study 
are  stable  and  only  slowly  evanescent. 

It  is  important  now  to  estimate  how  close  a  real  physical  system  can  be  made  to 
approach  the  conditions  for  force-free-field  approximations  on  the  vector  fields  E,  H,  and 
j,  if  electrons  were  magglera,  the  problem  would  be  far  simpler,  since  they  would  experience 
no  centripedal  accelerations  and  would  conform  to  electric  field  lines  like  currents  in  wires. 
In  reality,  the  electrons  cannot  confine  their  motion  to  be  strictly  along  the  curved  field 
lines  of  a  force-free  field. 

The  assumption  that  the  electrons  undergo  simple  harmonic  motion  along  electric  field 
lines  becomes  more  accurate  as  the  frequency  increases  and  the  electric  field  decreases.This 
can  be  seen  from  eq.(13),  where  the  amplitude  of  the  harmonic  motion  is  J/(peu>).  Sub¬ 
stituting  for  J  in  terms  of  E,  a  condition  for  small-amplitude  oscillations  can  be  written 


IE!”**  <  4;4-9?  (19) 


where  |E|max  is  the  Tna-rirmim  absolute  amplitude  of  the  harmonic  electric  field  within  the 
excitation.  Using  R.M.K.S.  units  and  a  radius,  flo  =  4.493/a  eq.(19)  becomes 


|E|m„  =  2.02xl07C(§),  '  (20) 

where  (  is  an  appropriately  small  factor  that  satisfies  eq.(19),  and  A  is  the  free  -space 
wavelength  of  the  harmonic  wave. 

The  assumption  will  be  that  the  smaller  the  value  of  C,  the  longer  the  lifetime  of  the 
excitation,  all  other  things  being  equal.  At  this  time  there  ere  no  estimates  for  liftimes, 
but  some  reasonable  values  for  the  parameters  might  be  Ho  =  10cm,  A  =  1.0mm,  and 
C  =  .01.  Under  these  conditions  a  10cm  radius  spherical  excitation  will  have  a  total  energy 
density  on  the  order  of  10*  Joules.  This  value  is  larger  than  average  reported  energies  of 
ball  lightning^,  and  these  excitations  may  represent  an  electromagnetically  excited  fluid 
state  that  could  serve  as  a  model  for  ball-lightning-like  phenomena. 

The  total  energy  in  any  spherical  excitation  will  be  proportional  to  the  factor 


r-(f)V  Pi) 
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This  implies  that  for  a  fixed  ratio  of  the  radius  to  the  wavelength,  for  the  same  value  of 
(,  the  energy  content  of  a  spherical  excitation  is  linearly  proportional  to  its  radius.  Since 
the  volume  is  linearly  proportional  to  the  cube  of  the  radius,  the  smaller  the  sphere,  the 
higher  the  energy  density,  under  the  above  fixed  ratio  of  radius  to  wavelength.  For  the 
numerical  example  already  given,  if  the  radius  were  shrunk  to  one  millimeter,  the  total 
energy  content  would  be  about  15,000  Joules.  The  wavelength  of  the  harmonic  wave  would 
be  10/im;  and  the  energy  equivalent  to  a  fifty-watt  lightbulb  burning  for  five  minutes. 

m.  CONCLUSIONS 

A  model  for  electromagnetic/ionized-fluid  excitations  predicts  the  self  containment 
of  high-energy  modes  in  small  volumes.  The  virial  analysis  yields  the  same  results  as 
self-consistent  solution  of  Maxwell’s  equations  with  current  solved  with  a  lagrangian  for¬ 
mulation,  allowing  the  current  to  be  carried  by  non-zero-maas  electrons.  The  important 
result  is  that,  excluding  electromagnetic  stess  effects,  the  most  stable  configuration  must 
have  a  uniform  free-charge  density  that  is  determined  only  by  the  size  of  the  excitation 
and  the  frequency  of  the  harmonic  electromagnetic  field. 

The  artifice  of  idealized  surface  currents  was  used  to  isolate  the  internal  fields  from  the 
outside  and  prevent  radiation.  Internal  electromagnetic  volume  stresses  were  ignored  in  the 
analysis.  Ongoing  work  is  indicating  that  if  k  and  a  are  chosen  correctly,  then  the  outgoing 
radiation  field  can  be  nullified.  Indications  are  also  that  the  internal  electromagnetic 
stresses  can  be  cancelled  by  vortical  fluid  modes. 

The  excitations  trill  have  finite  lifetimes,  but  may  persist  for  a  considerable  amount  of 
time  if  losses  are  not  large.  Upper  limits  were  set  on  the  maximum  electric  field  strength 
within  an  excitation,  to  confine  the  electronic  motion  to  electric  field  lines.  It  U  expected 
that  there  may  be  a  tradeoff  between  maximum  field  strength,  i.i.  maximum  energy 
density,  frequency,  and  lifetime.  Larger  field  strengths  and  lower  frequencies  imply  larger 
excursions  of  the  electrons  from  electric  field  lines,  leading  to  volume  losses  as  well  as 
surface  losses,  with  synchrotron  radiation  probably  being  very  important.  Vortical  modes 
may  greatly  affect  lifetimes. 

Using  reasonable  input  numbers,  maximum  field  strengths  and  subsequent  stored 
energies  are  very  high,  making  the  production  of  such  excitations  look  very  attractive.  It 
is  not  clear,  however,  what  a  production  facility  for  harmonic  force-free  excitations  might 
look  like,  since  the  internal  vortical  structure  of  these  excitations  is  still  under  investigation. 
“Fireballs'’  reported  in  the  literature*®*  may  be  candidates  for  investigating  their  internal 
fields  to  determine  whether  they  fit  the  harmonic  force-free  model. 

Theoretical  areas  still  need  to  be  researched.  Radiation  rates  due  to  non-linear  fluid 
effects  to  be  calculated.  These  rates  might  set  an  upper  bound  to  the  lifetime  of  a  clas¬ 
sical  excitation.  Although  a  rough  stability  analysis  was  done  in  a  previous  paper*3*,  a 
more  rigorous  analysis  involving  perturbations  of  the  current  remains  to  be  done.  This 
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would  involve  the  calculation  of  restoring  forces,  if  any,  that  would  arise  from  a  physical 
deflection,  as  well  as  other  mechanisms  that  would  return  the  excitation  to  its  original 
state  after  a  perturbation.  Energy  losses  due  to  perturbations  of  the  conditions  from  ideal 
should  be  calculated.  Much  of  this  depends  on  the  work  being  done  in  connection  with 
electromagnetic  stresses  influencing  fluid  motion. 

The  idea  of  applying  quantum  mechanics  to  these  excitations  can  be  investigated.  As 
the  size  of  the  excitations  shrinks,  quantum  effects  may  become  very  important,  according 
to  Nardi^10*.  If  a  self-consistent  potential  is  formed  by  the  currents  in  an  excitation,  there 
may  be  a  collective  electronic  wavef unction  that  exactly  confines  its  flux-field  vectors 
to  the  streamlines  of  a  force-free  electric  field.  Quantized  in  this  manner,  lifetimes  of 
states  may  be  much  larger  than  would  be  expected  from  purely  classical  considerations. 
Maximum  energy  levels  may  also  be  much  larger.  The  form  of  the  calculations  would 
probably  resemble  Hartree-Fock  atomic  procedures,  and  would  ultimately  have  to  include 
fluid-dynamic  effects. 

FVirther  investigations  seem  to  hold  great  promise. 
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ABSTRACT 

Spherically  symmetric  force-free  harmonic  solutions  to  Maxwell’s 
equations  may  be  truncated  at  a  radius  where  the  fields  are  tan¬ 
gential  to  a  sphere.  Equivalent  currents,  which  satisfy  the  jump 
boundary  conditions  between  the  inner  fields  and  the  vacuum,  main¬ 
tain  the  excitation  as  though  the  inner  solution  extended  to  infinity 
and  prevent  radiation  flow.  Since  the  Poynting  vector  is  identically 
zero  inside  the  spherical  boundary,  no  energy  can  escape  from  the 
inner  volume  even  if  the  surface  currents  were  not  present.  The  sur- 
.  -  face  boundary  fields  can  radiate  away  the  energy,  however,  which  is 
exactly  cancelled  by  the  radiation  from  the  surface  currents.  It  will 
be  shown  that  for  any  given  frequency  an  infinite  number  of  choices 
for  the  force-free  parameter  are  possible  such  that  the  far-field  radi¬ 
ation  from  the  surface  currents  is  zero.  Each  of  these  choices  implies 
a  different  spherical  truncation  radius.  For  large  wavenumbers,  suc¬ 
cessive  possible  choices  of  truncation  radii  differ  by  approximately  a 
half  wavelength.  For  any  of  these  choices,  no  surface  currents  are 
needed  to  prevent  steady-state  radiation. 


I.  INTRODUCTION 

It  has  been  shown  that  Maxwell’s  equations  for  harmonic  electromagnetic  waves,  with 
angular  frequency  u,  admit  solutions  having  the  electric  and  magnetic  fields,  E  and  H,  as 
well  as  the  current,  J,  in  the  force-free  form^.  A  field,  P,  has  the  force-free  form  if 


VxF  =  aF, 


(1) 


where  a,  the  force-free  constant,  is  not  a  function  of  either  space  or  time. 
Consider  the  following  solution  to  eq.(l)  in  spherical  coordinates  (r,$,j) 


1 


(2) 


Fj(r,  <t>)  =  (/  +  1)  costfsin1-1  9^^—^r  —  [ji_j(ar)  —  /M22l2l  ain1  $$ 

±  jl(ar)  sin1  9$, 

where  l  is  any  integer,  l  >  0,  j|  it  the  familiar  spherical  Bessel  function,  and  the  vectors  r, 
9,  and  4>  are  the  familiar  right-handed  3-D  polar  coordinate  unit  vectors.  This  is  not  the 
only  solution  to  eq.(l),  but  it  has  interesting  properties*1'.  To  be  specific,  let  I  as  1  and 
a  be  real.  At  ar  =  4.493 . . .,  ji(ar)  =  0  and  the  field  becomes  tangential  to  a  spherical 
surface,  with  radius 


iZo  =  4.493 . . .  /a.  (3) 

Both  E  and  H  can  be  written  as  vector  fields  proportional  to  F  in  eq.(2),  multiplied  by 
the  time  factor  exp(-iwt).  This  time  factor  will  henceforth  be  suppresed.  If  the  solution 
to  Maxwell’s  equations  is  truncated  at  r  =  iZo,  with  zero  fields  for  r  >  ii,,  the  resultant 
field  will  continue  to  be  a  solution  throughout  all  of  space  if  currents  are  made  to  flow 
on  the  surface  of  the  sphere  at  r  =  IZo.  The  equivalent  of  an  electric  current,  and  a 
magnetic  current,  M,  must  flow  to  satisfy  the  jump  conditions  in  the  fields  at  r  =  fZ». 
These  currents  are  defined  to  be 


J,  =  -axH„ 


(4) 


and 


M,  =  nxE„  (5) 

where  E,  and  H,  are  the  surface  electric  and  magnetic  fields  and  n  is  the  outward  normal 
at  r  =  fZj  .  The  physical  significance  of  these  currents  is  discussed  elsewhere*3'. 

JVom  eqs.(2-5),  the  surface  currents,  J,  and  M(,  will  be  proportional  to 


J.w/  =  jo(afZo)  sin  9j>, 


(6) 


where  jo(alZo)  =  —0.2172 . . . 


Under  ideal  situations  no  radiation  can  escape  from  the  spherical  surface  if  the  currents 
J,  and  M,  are  present.  The  radiation  from  these  currents  exactly  cancels  the  radiation 
from  the  surface  of  the  excitation.  The  reason  is  that  the  Poynting  vector, 
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S  =  ExH, 


(7) 


is  Miration, t  equal  to  zero  everywhere  farr<  R$. 

la  the  next  xetxa,  the  far-field  raSarion  from  the  surface  currents  will  be  analyzed 
fiar  the  f  *  l  cate.  It  will  be  shown  that,  fcr  appropriate  values  of  flo.  the  surface  currents 
will  ant  ratSate.  the  coochaioe  will  be  that  under  certain  conditions,  surface  currents  are 
not  aaW  to  ptswat  steady  state  radiation  from  an  excitation. 

EL  MILL  SURFACE  RADIATION  CONDITIONS 

Dwjala^the  aassniiirchss  4  •  j/c.  with  c  being  the  speed  of  light,  the  vector  potential 
in  the  far  laid  far  a  volume  i  if  trie  current,  J{r),  approaches  the  expression12) 

Mr)  -  “  jf  -W\  (8) 

where  f  it  a  unit  vector  in  the  direction  of  the  point  in  the  fax  field,  R  is  the  distance  of 
the  far-fidd  paint  from  an  arbitrarily  ehrmm  center  of  the  current  distribution,  r*  is  the 
vector  from  the  center  of  the  distribution  to  the  integration  point,  and 


|t'|cos7  =  fr’.  (9) 

In  eq.(8)  the  volume  integral  extends  over  the  entire  current  distribution  volume,  V. 

When  the  current  distribution  collapses  to  a  surface  distribution  like  that  in  eq.(6), 
the  integral  in  eq.(8)  becomes  proportional  to  the  integral 


I  =  jMrin  0'e-,**>co,W,  (10) 

where  the  distribution  is  t  entered  at  r*  =  0,  and  the  integration  is  over  the  surface  of  the 
sphere  (r1]  =  Rq. 

The  integral  in  eq.(8)  can  be  evaluated  in  closed  form  by  expanding  all  quantities  in 
the  integrand  in  terms  of  spherical  harmonics'4*  Yi,m(ft'),  in  the  primed  coordinates,  and 
in  the  unpriced  (far-field)  coordinates.  First  convert  the  spherical  vector  into 
a  sum  of  unit  Cartesian  vectors  x  and  y  in  the  x  and  y  directions; 
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#'  =  -xsin#'  +  ycos  #'. 


(11) 


Next  note  that 

sin  0' sin#'  =  t^/§[Vi.i(ft')  +  Vk.-i(fi')],  (12) 

and 

sin#1  cos#'  =  ^{rYuW  +  Yt^ltt)],  (13) 

The  exponential  in  the  integral  has  the  expansion^ 

e-.*R.co .7  =  4lgr«jl(iflo)  £  y^(ft')KIim( 0),  (14) 

J;*0  m»-I 

where  the  asterisk  (*)  denotes  complex  conjugation.  Using  the  orthogonality  property  of 
the  spherical  harmonics, 

fy, %(0')rJ.  ,m.(o')<«' = (is) 

a  little  algebra  shows  that,  finally, 

I  =  -47r;i(kflo)#,  (16) 

where  #  is  the  unit  spherical  vector  in  the  far-field. 

If  the  product  kRo  is  one  of  the  zeroes  of  ji(kAo)  then  the  vector  potential  for  outgoing 
waves  due  to  the  1  =  1  surface  electric  current  goes  to  zero  more  quickly  than  the  inverse 
of  the  distance  in  the  far-field.  This  current  will  therefore  not  radiate^3*.  An  exactly 
analogous  argument  shows  that  the  l  =  1  magnetic  surface  current  will  not  radiate  either, 
if  the  same  conditions  are  met. 
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In  the  next  section,  proof  will  be  given  for  the  general  result  that  the  outgoing  radiation 
field  bom  l  >  1  cases  can  also  be  nullified. 

III.  GENERAL  CONDITIONS  FOR  NULL  RADIATION 

In  the  1  =  1  case  the  integrations  needed  to  evaluate  the  far  field  simplified.  For 
any  other  allowed  value  of  l,  Z  >  1,  the  integrals  involve  sums  of  terms  of  the  type 
it(kRo)Fi,±i(n)i  with  L  odd.  Only  for  I  =  1  is  there  a  single  term  in  the  sum.  For  higher 
values  of  /,  then,  it  is  not  clear  whether  the  far-field  radiation  may  be  only  be  minimized, 
or  nulled  in  any  one  direction,  but  cannot  nulled  out  completely  in  all  directions.  The 
azimuthally  symmetric  l  =  1  mode  in  eq.(2),  and  modes  derived  from  it  by  rotations  and 
subsequent  linear  superposition,  are  not  the  only  ones  for  which  the  surface  radiation  of 
a  truncated  force-free  harmonic  excitation  can  be  exactly  nullified  by  a  judicious  choice 
of  the  force-free  constant  a  and  the  free-space  wavenumber  k  without  invoking  surface 
currents.  This  will  now  be  shown  by  comparing  the  fields  in  the  above  excitations  with 
the  internal  electromagnetic  oscillations  of  a  resonant  spherical  cavity® 

The  strategy  will  now  be  to  choose  excitations  that  obey  eq.(l),  choosing  Rq  and  k 
such  that 


jl(oRo)  =  0,  (17) 


and 


jl(kRo)  =  0.  (18) 


Define,  now,  the  complementary  vector  fields 


Mi  =  ji(kr)  sin'  0$.  (19) 


and 


N,  =  (i  +  l)cos0sin,-,0^r—  [jili(Jfcr)  - 1^-]  sin' 00.  (20) 

A  solution  of  Maxwell’s  equations  in  a  perfectly  conducting  hollow  sphere  of  radius 
r  =  Ro  is 


5 


(21) 


E  =  Mi,  H  =  -i— N(. 


If  eq.(18)  holds,  then  eqs.(21)  represent  a  resonant  excitation  of  a  spherical  cavity.  The 
fields  represented  by  eq.(21)  are  confined  to  to  the  region  r  <  Ho-  The  electric  field  at 
r  =  Ho  is  zero.  The  magnetic  field  induces  a  current  I  =  —  nxH  at  r  =  Ho,  where  n  is 
the  outward-pointing  normal  to  the  spherical  surface.  This  current,  considered  by  itself 
in  unbounded  free  space,  is  the  only  source  of  the  electromagnetic  field  both  inride  and 
outside  the  spherical  cavity.  As  such,  this  current  cannot  radiate  into  the  vacuum  region 
r  >  Ho,  so  that  the  electromagnetic  field  external  to  I  will  remain  zero  if  it  is  initially  zero. 

Now  note  that  at  r  =  Ho 


F|(Ho,M)  =  N,(Ho,M).  (22) 

Consider  first  the  magnetic  field  in  a  force-free  excitation  where  eqs.(17)  and  (18)  hold. 

from' the  surface  at  r  =  Ho,  due  only  to  the  magnetic  field,  will  entirely  depend 
on  the  contribution  from  the  term  I'(fl,  d)  =  nxH(Ho,«,  *),  where  n  is  the  outward  normal 
to  the  sphere  at  (Ho,M).  Since  n-H(Ho ,9,4)  =  0,  it  does  not  contribute*5'.  Rom  the 
above  argument  involving  I,  if  the  external  radiation  field  due  to  the  surface  contribution 
from  the  magnetic  component  of  a  force-free  excitation  is  initially  zero,  it  will  remain 
zero.  Exactly  analogous  arguments  apply  to  the  electric  field  contributions  to  the  external 
radiation. 

IV.  DISCUSSION 

The  'vyi'-lnoinn  to  be  drawn  is  that  ideal  spherical  electromagnetic  force-free  excita¬ 
tions  can  exist  with  no  e:.  xnal  radiation  field  for  all  values  of  l  in  eq.(2).  Equation  (2)  has 
no  variation.  It  is  not  difficult  to  generalize  Mi  and  Ni  in  eq.s(19)  <md  (20),  and 

thereby  Ft  using  the  formulas  on  page  416  of  Stratton’s  book*5'.  None  of  the  arguments 
change  when  azimuthal  variation  is  introduced,  so  that  all  the  results  derived  for  modes 
with  no  »vi»mithal  dependence  carry  over  to  the  modes  having  any  consistent  dependence 
on  the  angle  <t>. 

For  reasons  of  simplicity,  all  comments  will  henceforth  be  directed  toward  the  /  =  1 
r,M>  having  no  azimuthal  dependence.  This  is  the  mode  with  the  lowest  symmetry  and 
least  complicated  geometry.  The  surface  fields  topologically  are  eqivalent  to  invisrid-fluid 
flow-lines  over  a  rigid  sphere. 

Since  the  mode  was  truncated  at  the  first  zero  of  the  ji(or)  function,  it  is  easy  to  see 
that  a  necessary  condition  for  nullifying  surface  radiation  is 
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k>a. 


(21) 


Below  a  certain  frequency,  for  a  given  a,  the  surface  radiative  losses  cannot  be  nulled 
(unless  the  k  =  0,  or  zero-frequency  case,  is  allowed).  It  has  already  been  shown  that  far 
k>  a  the  interaction  energy  of  the  volume  harmonic  current  and  the  oscillating  magnetic 
field  has  a  negative  value  everywhere  within  the  excitation*3).  In  other  words,  the  current 
lies  in  a  negative  potential  well  whose  relative  depth  increases  with  frequency  for  k  >  a. 
This  potential  is  positive  for  k  <  a.  Although  it  has  not  yet  been  tested,  this  presumably 
means  that  ’.he  k  >  a  case  is  more  stable  than  the  it  <  a  case.  It  is  interesting  that  surface 
radiation  can  only  be  nulled  in  what  is  presumably  the  stable  region. 

Outside  the  excitation  there  are  null  fields. 

When  non-zero-mass  electrons  are  introduced  as  harmonic  current  carriers,  as  the 
actual  physical  situation  demands,  certain  energy-loss  mechanisms  can  be  identified*3). 
It  is  evident  that  if  force-free  harmonic  excitations  can  be  manufactured  in  a  form  that 
resembles  spherical  plasmoids,  it  is  most  likely  that  a  and  k  will  have  to  be  tailored  in  the 
manner  described  in  this  note,  to  prevent  radiation  losses. 
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